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CHAPTER  I.   INTRODUCTION 


FINLIN  (Mnite  element,  X^oparametric ,  N^on-Unear  with 
Interaction  and  N[o-tension)  is  a  finite  element  computer  pro- 
gram for  analysis  of  flexible  pipe  culverts  buried  in  soil. 
Curved  bar  segments  have  been  used  to  simulate  flexible  pipes. 
A  zero  thickness  special  type  of  frictional  element  is  used 
to  simulate  the  occurrence  of  slip  between  pipe  and  soil. 
Isoparametric,  linear-strain  triangular  elements  are  used  to 
represent  the  soil.   Non-linear,  anisotropic,  state  of  stress- 
dependent  soil  properties  expressed  in  terms  of  octahedral 
stresses  are  used.   Several  option  commands  permit  realistic 
analysis  of  a  culvert  problem  in  two  dimensions.   Incremental 
stages  of  construction,  controlled  slip  between  soil  and  pipe 
and  no-tension  analysis  can  be  performed.   The  program  has 
been  written  in  FORTRAN  IV  source  language  using  hard  and 
soft  wire  facilities  available  in  the  CDC  computer  system. 
In  this  report,  the  procedure  for  preparation  of  input  data 
for  a  problem  will  be  discussed;  a  complete  listing  of  the 
source  program  is  included.   Details  of  mathematical 
formulations,  developments  of  different  element  properties, 
and  related  information  can  be  found  in  the  Interim  Report  - 
"Predicting  Performance  of  Pipe  Culverts  Buried  in  Soil",  by 
M.  B.  Roy,  Purdue  University,  May,  1976,  JHRP-76-15. 

In  any  finite  element  analysis  the  first  step  is  to 
bound  the  problem  by  a  set  of  finite  boundaries  with  appropri- 
ate boundary  conditions.   Then  each  region  of  different 
materials  has  to  be  distinguished.   In  the  next  step,  each 
zone   is  subdivided  into  a  number  of  characteristic 
finite  elements  maintaining  continuity  at  the  boundary  be- 
tween two  elements.   Then  nodes  and  elements  are  numbered  in 
sequence.   Required  material  properties  for  every  element 
need  to  be  defined.   In  the  case  of  nonlinear  analysis,  where 
material  properties  change  depending  upon  the  state  of  stress, 
parameters  which  govern  the  variation  of  properties  need  to 
be  defined.   Also,  construction  in  layers,  load  application 
in  increments,  limited  shear  or  slip  in  the  pipe-soil  inter- 


action,  and  no-tension  in  soil  must  be  accommodated. 

Preparation  of  data  for  a  real  problem  should  take  ad- 
vantage of  the  fact  that: 

1.  The  boundaries  need  only  define  the  areas  of 
primary  interest,  and  obvious  features,  like 
symmetry,  should  be  recognized. 

2.  Smaller  sized  elements  are  needed  in  zones  of 
maximum  interest  and/or  high  stress  gradients; 
larger  elements  may  be  used  elsewhere. 

3.  Each  region  is  divided  into  an  appropriate 
number  and  type  of  elements  whose  node  points  are 
numbered  sequentially  .   The  cost  of  a  solution 
depends  heavily  on  the  numbering  sequence  of  nodes 
even  if  the  total  number  of  nodes  and  elements 
remain  the  same.   A  useful  general  rule  is  to 
minimize  the  maximum  difference  between  the  highest 
and  lowest  node  numbers  in  an  element. 

The  next  step  is  selection  of  the  type  of  analysis,  such 
as  number  of  construction  layers,  number  of  increments  of 
load  in  a  given  layer,  interaction  properties,  no-tension 
in  soil  and  other  similar  decisions.   All  information  has  to 
be  digitized  and  checked  for  correctness. 


CHAPTER  II.   COMPUTER  PROGRAM 

1.   Types  of  Finite  Elements 

Several  types  of  finite  elements  have  been  used  in  the 
program  FINLIN.   Description  of  each  type  follows. 

Type  I,  Curved  Bar  Element 

Segments  of  a  ring  have  been  used  to  represent  flexible 
pipe,  which  has  small  thickness  compared  to  the  radius. 
Figure  1  shows  a  typical  Type  I  element  with  two  nodes, 
each  node  having  three  degrees  of  freedom  -  radial,  tangential 
and  rotational.   The  radius  of  the  pipe,  its  stiffness  EI 
(young's  modulus  times  moment  of  inertia),  and  the  nodal  co- 
ordinates need  to  be  defined.   The  position  of  the  center  of 
curvature  is  also  necessary. 

Type  II,  Interaction  Element 

This  is  a  zero  thickness,  rectangular  element  with  four 
nodes,  each  node  having  two  degrees  of  freedom  in  the  normal 
and  tangential  directions.   Figure  2  shows  a  typical  inter- 
action element  where  the  coordinates  of  nodes  1  and  4  (and 
2  and  3)  are  initially  the  same.   Coordinates  of  all  four 
nodes,  stiffness  values  in  normal  and  tangential  directions, 
E   and  E   respectively,  need  to  be  defined.   The  program  will 
modify  the  values  of  E   and  E   according  to  the  state  of 
stress  in  these  elements. 

For  both  Type  I  and  Type  II  elements,  the  program  will 
perform  the  necessary  coordinate  transformation  depending 
upon  their  position  and  orientation. 


Type  III,   Isoparametric  Triangular  Element 

This  type  of  triangular  element  has  three  corner  nodes 
and  three  intermediate  (midpoint)  nodes,  each  node  having 
two  degrees  of  freedom  (Figure  3).   The  face  1-4-2  can  be  a 


Figure  I,  Type  - 1,  Curved  Bar  Element 
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Figure  2.  Type-H,  Interaction  Element 


Figure  3,  Type -III ,  Isoparametric  Triangular  Element 


curved  boundary.   Nodes  of  the  Type  III  element  must  be  de- 
fined in  the  order  shown  in  Figure  3.   Coordinates  of  only 
corner  nodes  need  to  be  defined.   This  element  is  used  to 
simulate  the  soil  medium  and  for  thick  pipes  where  rotation 
at  pipe  nodes  are  not  significant.   Required  material 
properties  for  this  type  of  element  are  modulus  and  poisson's 
ratio.   Depending  upon  the  nature  of  the  material,  linear  or 
nonlinear  properties  can  be  used. 

In  case  of  analysis  with  no  interaction  element  between 
pipe  and  soil,  placing  Type  III  element  adjacent  to  Type  I 
elements  causes  a  difficulty  in  numerical  procedure  because 
nodes  of  2  and  3  degrees  of  freedom  lie  at  one  point.   To 
eliminate  this  problem  Type  IV  and  Type  V  elements  are  used. 
Type  IV  element  (Figure  4)  has  the  face  1-4-2  adjacent  to  the 
pipe.   For  Type  V  element  (Figure  5)  node  3  touches  the  pipe. 
Improper  representation  of  elements  will  cause  abnormal 
termination  of  the  program. 


Material  Properties  for  Type  III  Elements 

The  program  is  capable  of  using  linear,  nonlinear  and 
anisotropic  material  properties.   Also  different  types  of  soils 
with  distinct  properties  can  be  used.   For  linear  materials  only 
values  of  Young's  modulus   E  and  Poisson's  ratio  v  need  to  be 
defined.   For  anisotropic  material,  the  ratio  of  moduli  in 
vertical  to  horizontal  direction  has  to  be  defined.   For  non- 
linear properties,  experimental  data  are  directly  used.   It  is 
to  be  noted  that,  in  this  program  tangent  modulus  and  tangent 
Poisson's  ratio  values  are  used  for  incremental  analysis.   Also 
octahedral  normal  and  shear  stresses  have  been  used  in  the 
formulation.   Nonlinear  soil  properties  for  any  value  of  stress  level 
and  stress  ratio  are  interpolated  using  cubic  spline  functions. 

It  is  necessary  to  convert  conventional  test  data  (e.g. 
Figure  6)  to  a  form  which  is  acceptable  to  this  program,  as 
fol 1 ows : 


YA 


Side    1-4-2   curved 

Nodes  1-4-2  have 
3 -degrees  of  freedom 


Figure  4»  Type  -  IV,  Triangular  Element 


YA 


Node  3  has    3- 
degrees  of  freedom 


Figure    5,  Type-  V,  Triangular   Element 
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Figure  6,  Plane   Strain  Test  on  Loose   Monterey 

No.O    Sand   (after   Lade.  1972) 
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1.  Select  successive  values  of  {a^-a^)*    e^  and  e^  from 
the  test  data,  at  a  given  o^   starting  from  zero. 

2.  Cubic  splines  are  fitted  to  the  data. 

3.  Small  increments  of  e,  are  chosen,  and  values  of 
a„,  E.  and  v.  are  computed  from  the  generalized 
Hooke's  Law.   With  a 


a      .    and  t  .  and 
oct      oct 


^,  a«,  and  o^   known,  values  of 

"^     ^/ "^     X  r  -1    ai^e  calculated, 
oct  oct  failure 

4.  Step  (3)  is  repeated  for  each  increment  of  e,  up  to 
fai 1 ure. 

5.  Steps  (1)  through  (4)  are  repeated  for  all  curves 
with  different  a. 


3* 
All  values  of  E.  are  plotted  against  corresponding 

a^  .,  and  values  of  stress  ratio 
oct 


( 


QCt 


) 


^oct  fail ure 
are  noted  for  all  points. 

7.  Contour  lines  are  drawn  for  selected  values  of  stress 

ratio.   For  convenience,  values  of  stress-ratio  are 

made  to  range  from  zero  to  unity  with  increments  of 

0.1.   This  will  generate  11  curves  of  E.  vs  a„^.  for 
^  t     oct 

stress  ratio  from  0  to  1  (see  Figure  7). 

8.  Procedures  similar  to  steps  (6)  and  (7)  are  employed 
for  tangent  Poisson's  ratio,  v^,  (see  Figure  8). 

To  aid  this  procedure  of  data  reduction,  a  small  computer 
program  called  "PROPRTY"  can  be  helpful  which  will  cover  steps 
(1)  through  (5)  and  prints  out  data  required  for  steps  (6)  and 
(8). 


Description  of  the  Computer  Program 

The  program  FINLIN  has  been  divided  into  several  primary 
sections  which  are  called   OVERLAYS.    Each  OVERLAY  performs 
specific  computations  and  stores  the  results  for  future  use. 
Figure  9  shows  the  principal  organization  of  the  program  and 
several  OVERLAYS.   Each  OVERLAY  consists  of  one  main  program 
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and  several  subroutines.  OVERLAY  1  and  3  are  the  two  most 

important  sections  of  the  program.   Names  of  subroutines 

associated  with  OVERLAY  1  and  3  are  shown  in  Figure  10.   A 

general  flow-chart  of  the  program  FINLIN  is  shown  in  Figure  11. 

Brief  descriptions  of  the  computations  p3rformed  in  OVERLAY 

(1.0)  and  (3.0)  are  given  below. 

INPUTl  -       Reads  nodal  and  element  data  for  all  elements, 

material  properties,  pipe  geometry  and 

properties . 
INPUT2  -       Reads  modulus  and  Poisson's  ratio  value  for 

linear  elastic  materials. 
PLSTRS  -       Computes  elasticity  matrix  for  plane  stress 

condition. 
PARAMTR  -      Computes  elasticity  matrix  for  plane  strain 

condition. 
AREAS  -        Area  of  triangular  elements,  and  semi-band 

width  for  global  stiffness  matrix  are  evaluated 

in  this  routine. 
SPLINE  -       If  nonlinear  properties  are  specified,  this 

routine  reads  data  for  non-linear  analysis  and 

stores  them  in  a  convenient  form  for  future  use. 
SPLFIT  -       Cubic  spline  fitting  for  nonlinear  material 

properties  is  the  primary  function  of  this 

routine . 
COFRIT  -       Generates  coefficients  for  cubic  spline  function 
TRIDGNL  -      Solution  of  tridiagonal  equation  for  spline  fit. 
FD  and  BD  -    Foreward  and  backward  interpolation  formulae  to 

determine  slope. 
ORDINET  -      Spline-function  interpolation  formula. 
FAILURE  -      Function  defining  stress  ratio  at  failure. 
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OVERLAY  (0,0) 
READ  TITLE,  INITIALIZE 
VARIABLES,  CONTROL  PROGRAM  FLOW 


I 


OVERLAY  (1,0) 

READ  NODAL  AND  ELEMENT  DATA, 

COMPUTE  NODAL  LOADS,  CALCULATE 

BAND  WIDTH,  READ  MATERIAL 

PROPERTIES,  SPLINE  FITTING. 


OVERLAY  (2,0) 
PLOT  FINITE  ELEMENT  MESH 


OVERLAY  (3,0) 
FIND  PROPERTIES,  EVALUATE  STIFFNESS, 
ASSEMBLE  STRUCTURE  STIFFNESS  MATRIX 
AND  LOAD 


I 


OVERLAY  (4.0) 
DETERMINE  STRAIN  AND  STRESS,  ADD 
RESULTS  FOR  INCREMENTAL  ANALYSIS, 

STORE  RESULTS 


T 


OVERLAY  (5,0) 
PRINT  OUTPUT  RESULTS  OR  RECYCLE  FOR 
INCREMENTAL  ANALYSIS 


OVERLAY  (6,0) 
MODIFY  STIFFNESS  MATRIX  FOR  BOUNDARY 
CONDITION  AND  SOLVE  FOR  INCREMENTAL 

DISPLACEMENTS 


Figure  II    General  Flow  Chart  of  Program  FINLIN 
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BLOCK 


TRANGLE 
RING  - 
JOINT  - 


This  is  the  most  important  routine  in  FINLIN. 
Primary  functions  are  (1)  selection  of  elements 
in  a  given  layer,  (2)  evaluation  of  nodal  loads 
due  to  a  qiven  layer  and  number  of  increments 
of  load,  (3)  generation  of  element  stiffness 
matrix  for  all  three  basic  types  of  elements, 

(4)  modification  of  appropriate  material 
properties  for  qiven  state  of  stress  condition, 

(5)  formation  of  global  stiffness  matrix  and 
load  vector,  (6)  storing  the  matrix  on  tape  for 
future  use. 

Evaluates  element  stiffness  matrix  for  linear- 
strain-tri anqul ar  elements  (Type  III,  IV,  V). 
Creates  stiffness  matrix  for  curved-bar  or  ring 
elements  (Type  I). 

Evaluates  stiffness  matrix  for  interaction 
elements  (Type  II). 


In  this  program  several  means  have  been  used  to  make 
efficient  use  of  memory  spaces.   For  this  purpose  the  same 
locations  have  been  used  several  times  for  entirely  different 
purposes.   If  a  user  wants  to  modify  or  change  any  portion  of 
the  program,  he  should  make  sure  that  the  intended  changes  do 
not  wipe  out  or  overwrite  a  portion  of  memory. 


3.   Limitations  of  the  Program 

In  its  present  form  the  program  takes  about  112,000 
(octal)  core  memory  spaces   to  load  and  execute.   FINLIN  uses 
nine  disc  tape  units,  out  of  which  two  units  are  used  for  INPUT, 
OUTPUT,  five  for  random  access  mass  storage  and  two  for  inter- 
mediate storage  purposes.   Limitations  of  the  program  are 
listed  below  (they  can  be  modified  easily  by  changing  a  few 
cards ) : 
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(1)  Total  number  of  nodes:   (NNODES)  =  550 

(2)  Total  number  of  elements:   (NELEMNT)  =  250 

(3)  Different  types  of  materials:   (MATRIAL)  =  25 

(4)  Increments  of  stress  ratio  in  Figure  7  and  8  =  0.1 

(5)  Maximum  number  of  points  for  nonlinear  material 
properties  =  7 

(6)  Maximum  semi -bandwi dth  for  global  stiffness  matrix 
(including  diagonal)  =  103 

However,  semi -bandwi dth  capacity  can  be  modified  easily 
by  making  changes  as  follows: 

(a)  Estimate  the  maximum  semi -bandwi dth  (including  the 
diagonal  term).   If  it  is  even,  add  one  and  make 
it  odd  numbered,  which  is  the  final  semi -bandwi dth 
(N)  of  the  system.   Numerals  in  the  indicated 
places  should  at  least  be  N. 

(b)  OVERLAY  (0.0),  Line:  FLN-36,  NSIZE  =  N 

(c)  OVERLAY  (3.0),  Line:  BLK-15 
DIMENSION  A(N,N)  ,  ARRAY  (N) 

(d)  OVERLAY  (6.0),  Line:  SOL-10 
DIMENSION  A(N,2N),  B(2N),  ARRAY(N) 

(e)  SUBROUTINE  MODIFY,  Line:  3 
DIMENSION  A(N,2N),  B(2N). 
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CHAPTER  III.   INPUT  DATA  CARDS  FOR  PROGRAM  "FINLIN" 


1st  Card  Type 

Col.  2-78: 
Co.  79-80: 

2nd  Card  Type 


Col.  1-5: 

Col.  6-10: 

Col  .  11-13 

Col.  14-16 

Col.  19-20 

Col.  23-24 


3rd  Card  Type 


Col.  1-5: 
Col.  6-15: 
Col.  16-75 


TITLE,  MESH 

FORMAT  (13A  6,  12) 

One  card  identifying  the  problem. 

alphanumeric  description  of  the  problem  to  be 

printed  in  the  output. 

if  greater  than  0,  a  plot  of  finite  element 

mesh  is  generated  by  CALCOMP  plotter.   If  zero 

or  blank,  no  plot  is  generated. 

riNODES,  NELEMNT,  MATRIAL,  NPRSR,  LAYERS,  ISTOP 

FORMAT  (215,  213,  2X,  12,  2X,  12) 

One  card  defining  problem. 

NNODES  -  Total  number  of  node  points. 

NELEMNT  -  Total  number  of  elements  (all  types 

incl uded) . 

MATRIAL  -  Total  number  of  different  soil  types. 

NPRSR  -  0 

LAYERS  -  Number  of  construction  layers. 

ISTOP  -  If  zero,  linear  elastic  soil  properties, 

if  greater  than  zero,  nonlinear  soil  properties. 

NTYPE,  GAMA,  TYPE 

FORMAT  (15,  FIO.O,  10A6) 

Soil  type  cards,  one  card  for  each  soil  type, 

total  number  =  MATRIAL 

NTYPE  -  Soil-type  identification  number. 

GAMA  -  Unit  weight  of  soil. 

TYPE  -  Alphanumeric  description  of  this  soil 

type. 
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4th  Card  Type:   ANLSIS,  DELTA 

FORMAT  (A6,  F5.0) 

One  card  specifying  type  of  analysis  and  angle 

of  friction  between  soil  and  pipe. 
Col.  1-6:     PLSTRS  -  For  plane-stress  analysis. 

PLSTRN  -  For  plane-strain  analysis. 
Col.  7-11:    DELTA  -  Anqle  of  friction  between  pipe  material 

and  soil  adjacent  to  pipe,  in  degrees. 

5th  Card  Type :   Two  cards  per  soil  material  type,  total  number  = 
2X  MATRIAL,  specifying  initial  modulus  and 
Poisson's  ratio  for  each  soil  type. 

(a)  1st  Card  -  E 

FORMAT  (EIO.O) 
Col.  1-8:     E  -  Initial  modulus. 

(b)  2nd  Card  -  NUE 

FORMAT  (F5.0) 
Col.  1-5:    NUE  -  Initial  Poisson's  ratio. 

6th  Card  Type:   These  cards  are  required  only  for  nonlinear 
materials  i.e.  if  ISTOP  is  greater  than  zero 
(in  2nd  Card  Type).   This  set  of  cards  is 
repeated  for  each  nonlinear  soil  type.   If 
ISTOP  =  0  or  blank,  these  cards  are  not  required 

(a)  NP,  PSI,  PHI,  ANISO,  DELTA,  FACTOR 
FORMAT  (15,  5F10.0) 
One  card  identifying  soil  properties. 

Col.  1-5:     NP  -  Number  of  points  on  each  tangent  modulus 
and  tangent  Poisson's  ratio  vs.  o      ^    curve. 

Col.  6-15:    PSI  -  Factor  defined  as 


,0^+03 


Col.  16-25:   PHI  -  Friction  angle  for  soil,  in  degrees 
Col.  26-35:   ANISO  -  Anisotropy  ratio, 

E 

y 

If  isotropic,  ANISO  =  1.0 
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Col.  36-45:   DELTA  -  Angle  of  friction,  in  degrees  between 

pipe  and  soil. 
Col.  46-55:   FACTOR  -  Conversion  factor  for  o^^^    and  tangent 

modulus,  e.g.  a   ,.  =  FACTOR  *  a^^^  and 

E^  =  FACTOR  X  E^. 


If  no  conversion  is  required,  FACTOR  =  1.0 


(b)  XP(I) 


FORMAT  (BFin.O) 

a      .    values  of  nonlinear  property  cards,  total 

OCt  f        t  J 

number  values  =  NP,  up  to  8  values  per  card. 

Col.  1-10:    MP(1)  -  1st  value  of  a^^^. 

Col.  11-20:   NP(2)  -  2nd  value  of  a^^^. 
and  so  on. 

(c)  EP,  PSNR 

FORMAT  (2F10.0) 

These  set  of  cards  define  the  nonlinear  material 
properties.   Two  sets  of  curves  are  required 
(1)  for  tangent  modulus,  E^  vs  o      .    for  stress- 
ratio  ranging  from  0.0  to  1.0  and  (2)  same  type 
of  curves  for  tangent  Poisson's  ratio  v^,  each 
curve  of  E.  and  v.  is  defined  by  NP  number  of 
points,  so  total  number  of  cards  =  NP  x  11. 

Col.  1-10:    EP  -  Tangent  modulus  value. 

Col.  11-20:   PSNR  -  Tangent  Poisson's  ratio  value 

Note:        (1)  First  NP  cards  should  read  the  values  of 

E   and  v.  for  increasing  values  of  cr   .  starting 

from  a   .  =  0. 
oct 

(?)  The  first  set  of  NP  cards  are  for  curve  of 
stress  ratio  =  1.   The  second  set  of  NP  cards 
will  stand  for  stress  ratio  =  0.9  and  so  on. 
The  last  set  of  NP  cards  will  read  values  of 
E   and  v.  for  stress-ratio  =  0.0. 
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(3)  A  set  of  cards  which  includes  cards  from 
(a)  to  (c),  define  a  complete  set  of  Type  6 
cards.   As  no  material  number  has  been  attached 
to  it,  the  sequence  in  which  the  set  of  cards 
are  placed,  will  define  the  soil  type.   For 
example,  the  first  complete  set  of  type  6  cards 
(which  includes  (a),  (b)  and  (c)  type  cards) 
will  automatically  be  defined  for  Type  1  (NTYPE 
in  3rd  Card  Type)  soil  and  the  second  set  of 
type  2  and  so  on. 

7th  Card  Type:   NREAD,  FACTOR 

FORMAT  (15.  FlO.n) 

Col.  1-5:    NREAD  -  Total  number  of  node  point  data  cards. 
The  trianqular  finite  elements  have  six  nodes 
i.e.  three  corner  nodes  and  three  mid-side 
nodes.   Co-ordinates  of  mid-side  nodes  are 
calculated  by  the  program.   Except  for  defining 
boundary  conditions,  these  nodes  need  not  be 
defined. 

Col.  6-15:    FACTOR  -  Conversion  factor  for  node  point 
coordinates  such  as  Ft.  to  Meter.   If  no 
conversion  is  required,  FACTOR  =  1.0. 

8th  Card  Type:   NN,  NCODE,  X,  Y 

FORMAT  (15,  5X,  I?,  FIO.O,  20X,  FIO.O) 
Nodal  cards,  one  card  per  nodal  point,  total 
number  =  NREAD. 

Col.  1-5:    NN  -  Node  point  number. 

Col.  11-12:   NCODE  -  Node  point  boundary  condition.   The 
following  table  describes  values  of  NCODE  to 
represent  desired  boundary  condition. 
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NCODE    Node  Type 


Boundary  Condition 


?  degrees  of 
freedom  in  X 
and  Y 
di  recti  ons 


Free  X  and  Y  directions. 
Fixed  in  X  and  free  in  Y 

di  recti  ons  . 
Free  in  X  and  fixed  in  Y 

di  recti  ons . 
Fixed  in  both  X  and  Y 

di  recti  ons . 


3  deqrees  of 

Free  in  X ,  Y  and  9 

freedom  in  X, 

di  recti  ons . 

Y  direction 

Fixed  in  X,  free  in  Y 

and  in  rotation 

and  9  directions. 

(e) 

Fixed  in  Y,  free  in  X 

and  9  di  recti  ons . 

Fixed  in  9  ,  free  in  X 

and  Y  directions. 

Fixed  in  X  and  9,  free 

in  Y  directions. 

Fi  xed  i  n  Y  and  9  ,  free 

in  X  directions. 

Col.  13-22:   X  -  X  co-ordinate  of  node  point 

Y  -  Y  co-ordinate  of  node  point 

Nodal  data  cards  may  be  placed  in  any  order. 


9th  Card  Type:   lEL,  IX 

FORMAT  (13,  2X,  2(315,  5X),  4X,  II,  3X  ,  12) 
Element  cards,  one  card  per  element  total 
number  =  NELFMNT  Cards. 
lEL  -  Element  number. 

IX(1)  -  Node  number  1  of  element  no.  lEL. 
IX(2)  -  Node  number  2  of  element  no.  lEL. 
IX(3)  -  Node  number  3  of  element  no.  lEL. 
IX(4)  -  Node  number  4  of  element  no.  lEL. 
IX(5)  -  Node  number  5  of  element  no.  lEL. 
IX(6)  -  Node  number  6  of  element  no.  lEL. 
IX(7)  -  Element  type  identification  number. 
IX(8)  -  Material  type  identification  number. 
Element  types  1,  2,  3  and  corresponding  node 
numbers  are  shown  in  Figure  1  to  Figure  3. 
This  numbering  scheme  has  been  used  in  the 
program.   For  type  4  and  type  5  elements  see 
Figures  4  and  5. 
The  element  data  cards  can  be  placed  in  any  order. 


Col  . 

1-3: 

Col  . 

6-10: 

Col  . 

11-15 

Col  . 

16-20 

Col. 

26-30 

Col  . 

31-35 

Col  . 

36-40 

Col  . 

50: 

Col. 

54-55 
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10th  Card  Type:  ZAI(l),  ZAI(2),  ZAI{3) 

FORMAT  (3F10.0) 

One  card  defining  the  point  where  stress  of 

triangular  element  is  required.   (See  Figure 

12  for  definition  of  ZAI). 
Col .  1-10:    ZAI(1 )  - 

Col.  11-20:   ZAI(2)  -      Area  coordinates. 
Col.  21-30:   ZAI(3)  - 

nth  Card  Type:  XCFN,  YCEN,  RADIUS,  EI 
FORMAT  (3F10.0,  F10.0) 

One  card  defines  center  of  circular  pipe, 
radius  and  stiffness. 
Col.  1-10:    XCEN  -  X-coordinate  of  center  of  pipe. 
YCEM  -  Y-coordinate  of  center  of  pipe. 
RADIUS  -  Radius  of  pipe. 

EI  -  Stiffness  of  pipe.  Young's  modulus  times 
moment  of  inertia  of  pipe  cross  section  per 
unit  length. 

12th  Card  Type:  This  card  is  required  only  when  a  plot  of 

finite  element  mesh  of  the  problem  is  required, 
which  is  specified  in  1st  Card  Type.   If  Col. 
79-80  is  zero  or  blank  no  plot  is  generated. 
If  any  number  is  punched  in  Col.  79-80,  a  plot 
will  be  generated  and  in  that  case  only  this 
card  is  required. 
XMAX,  YMAX,  YCEN 
FORMAT  (3F10.0) 
Col.  1-10:    XMAX  -  Maximum  size  of  mesh  in  X-direction. 
YMAX  -  Maximum  size  of  mesh  in  Y-direction. 
YCEN  -  Y-coordinate  of  center  of  pipe. 
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Figure   I2»  Definition    of  Area    Co-ordinates 
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13th  Card  Type 


Col.  1-2 


Col .  3-4 


Col  . 

5-6: 

Col  . 

11-20 

Col  . 

21-30 

Col  . 

31-40 

Col.  41-50 


NOTENSN,  INTER,  NSTEP,  Hi,  H2,  KN ,  KS ,  NITER 

FORMAT  (312,  4X,  2Fin.0,  2E10.0,  15) 

This  card  contains  very  important  information 

regarding  type  of  analysis  required.   Total 

number  of  cards  =  LAYERS  (as  specified  in 

Col.  19-20  of  2nd  Card  Type). 

NOTENSN  =  0,  if  soil  is  allowed  to  take  tension. 

=  1,  if  'no-tension'  in  soil  analysis 
is  required. 
INTER  =  0,  if  'no-interaction'  between  soil  and 

pipe  is  required  in  type  II  elements. 

In  this  case,  pipe  and  soil  are 

rigidly  connected. 
=1,  if  interaction  between  soil  and 

pipe  is  desired  which  will  permit  slip 

depending  upon  the  state  of  stresses 

in  type  II  elements. 
NSTEP  -  Number  of  equal  increments  of  application 
of  gravity  load  in  a  particular  layer. 
HI  -  Starting  height  of  a  construction  layer 
being  analyzed. 

H2  -  Finish  height  of  a  construction  layer 
being  analyzed. 

KN  -  Normal  stiffness  of  pipe-soil  interaction 
element  (type  II)  before  failure  in  tension. 
In  case  tension  develop  in  an  interaction 
element,  the  program  will  modify  value  of  this 
stiffness . 

KS  -  Shear-stiffness  of  pipe-soil  interaction 
element  (type  II)  before  failure.   If  INTER  =  0, 
value  of  KN  and  KS  are  kept  unchanged  in  all 
interaction  elements  throughout  the  analysis, 
which  simulates  rigid  connection  if  values  of 
KN  and  KS  are  considerably  high.   If  INTER  >  0 
values  of  KN  and  KS  are  modified  based  on 
failure  conditions  specified  for  an  interaction 
element. 
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Col.  51-55:   NITER  -  Maximum  number  of  iterations  specified 
for  'no-tension'  analysis.   If  the  solution 
does  not  converge  after  specified  number 
(=  NITER)  of  iterations,  farther  execution  will 
be  stopped  and  a  message  will  be  printed. 

Input  Data  Cards  for  Program  'PROPRTY' 

In  program  FINLIN,  an  octahedral  stress-strain  system  and 
tangent  modulus  and  tangent  Poisson's  ratio  values  have  been 
used  for  representing  nonlinear  soil  properties.   Required 
data  for  nonlinear  properties  in  FINLIN  (data  card  type  6),  is 
difficult  to  get  from  conventional  triaxial  or  plane  strain 
tests . 

In  program  FINLIN,  some  routine  calculations  and  inter- 
polations have  to  be  performed  to  prepare  data  for  card  type  6. 
The  program  PROPRTY  has  been  written  to  aid  in  generation  of 
the  data  required  for  interpolation.   This  program  accepts 
actual  test  data  and  interpolates  using  spline  function,  prints 
values  of  octahedral  stress,  strain,  stress-ratio,  tangent 
modulus,  tangent  Poisson's  ratio  and  similar  informations. 

1st  Card  Type:   TITL 

FORMAT  (lOAP) 
Col.  2-80:    TITL  -  Alphanumeric  identification  of  soil  and 
test;  one  card. 


2nd  Card  Type:   TEST,  NOCURVS,  RF 

FORMAT  (A6,  13,  FIO.O) 

One  card  defines  type  of  test  performed  and 
other  information. 
Col.  1-6:     TEST  =  PLSTRS  for  plane-stress  test  and 

=  PLSTRN  for  plane-strain  test  performed 
on  the  sample. 
Col.  7-9:     NOCURVS  -  Number  of  confining  pressures  used  in 

the  test;  should  at  least  be  three. 
Col.  10-19:   RF  -  Ratio  of  t^^^/o^^^  ^^    failure. 
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3rd  Card  Type:   SIGMA3,  NP 

FORMAT  (5X,  FIO.O,  15) 

One  card  specifies  value  of  confining  pressure 

and  number  of  points  for  this  test  with  given 

confining  pressure. 
Col.  6-15:    SIGMA3  -  Confining  pressure. 
Col.  16-20:   MP  -  Number  of  points  on  a  curve  of  given  a^. 

4th  Card  Type:   XP,  YP,  VP 

FORMAT  (3F5.0) 

One  card  per  point,  total  number  =  NP;  containing 

information  about  stress-strain  and  volume 

change  data  for  a  given  a-. 
Col.  1-5:     XP  -  Axial  strain  in  percent  (e,  %) . 
Col.  6-10:    YP  -  Deviator  stress  (a^-a^)  corresponding  to 

axial  strain,  XP. 
Col.  11-15:   VP  -  Volume  strain  in  percent  (e  %) . 

Note:   3rd  and  4th  card  types  are  repeated 

NOCURVS  times. 


APPENDIX  -  I 


FINLIN 
PROGRAM  LISTING 
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aV£RL.HV(RDY«Os  ij> 

aVERLAV'TRaY!  0.  IV) 
PRQGRrM    F  I HL ,[  ri    ■::  1  tiPUT  >  DUT  PUT  -  ThPF  2  -  TAPE9  -  ThPE':.  =  I HPUT  -  TAPEt'.=DIJTPUT 
i TAPE  1 ..  ThPE 3 .  TriPEa .  TAPE? «  i  APE  1  0 >  i ■■LOT > 

THIS   FINITE   ELFMENT   PRCIGRhM   HAS   EEEH   WRITTEN  FHR   PH.D=    THESIS 

PREiilCTIMG   PERFDRpiANCi:    DF    PIPE   CULVERTS    IN  SOIL 

EY   M.    B,    ROY,    GPA.UiJATE    RESEARCH   ASSISTANT 

SCHDOL    DF   CIVIL    EHhlNEER JiTG?    GEDTECHHICAl    ENGINEERING? 

PURHUE    UHIVERSTF.  <    west    LAFAYETTE.     INBIANA    -    47907 

ri;-iTE    BEl  EMBER.     1975. 

CDHMDN  HNDDES.  HElEMHT?  NBlif"?  MBAHD?  NB5  NTS;  I  STOP  <  NCYCLE;  LAYERS?  I  STEP 
lNSTEp5NT!ajETH.HTl.NT£.NnTENSN-IFLAG!NSIZEjNCDDE'::55i:0!X<55n>,Y<55 
2  > !  JND>i  ■::  5  i  J' «  ANLS I S  s  I X  <  8  ?  £!  ■  U  >  ?  ARE  Aft  <  £5  0  >  «  I NDI'-;  C  £5  0 ;'  •  I NDEX  ■:;  £5  0  >  ?  GAMA  <  i 
35::'jZAIt:3> 

CDMMDf i   -■■  !-'■   t  <  NUE ,  RAD T  US 5 ;  CEH .  YiJEN 5  E I  ?  KN  -  KS  -  H  i  •  H£ ,  I  f  jTER 

CUPlMDH    ■■£■■'■   D 

CDriPiuri   ^3-'   P 

CDMhuH   -'4.''   PiiJliU(£6)<DELTr- 

COMMDri    .■•■B-'    Q  <  1  i  0  0  >  »  L I S-T  ■;    i  0 1  '■■ 

DIMENSIuN    TITLE'::i3>!    R<liu0>n    irfip,;:i7> 

BIMENSIDH   E(10>3    NUE(10)i    IKlOiiiO) 

REAL  NUE 

REAL  KNjKS 

THIS  IS  pIhIh  UVERLhY?  WHICH  BIRECTS  El.ECUTinN  DF  OTHER  DVERLAYS 
DEPtruiING  UPDN  T' "PE  OF  hNALYSIS  HtID  DTHER  COMMANDS 
FuLLuWInG  TWO  STATEMENTS  ARE  LThRhRY  ROUTINES  FDR  FDR  BLDCK  hND  S 

CALL  FTNBIri  (l?  0> 
CALL  SETSTAK  (U> 

MAX.  SIZE  DF  SEMI-BAND  WiDTHQF  rOTAL  STIFFNESS  MATRIX  INCLUDING  B 

NSIZE=i03 

NSTEP=i 

IFLAG=U 

REWIND  S 

READ  (5)55:5  i  iTLEnMESH 

IF  CEDFjS)  5.iO 
i  GO  TO  40 
I  CONTINUE 

CALL  ZERO  <  R  ?  1 1  0  0  ;■ 

WRITE  <£>  <R'' T  :.  ?  1  =  1  5  1  i  00  ■ 

ISTEP=i 

WRITE  (t-fcO)  riTLE 

CREhTE  mass  STORAGE  FILES 
INITIALISE  MASS  STORAGE  UNITS 

CALL  OPENMS  < J  « INBXj  £50? lO 
CALL  DPEHMS  (S?  INDEXj  £50?  O!) 
CALL  OPENMS  'M. MDDU? £6- 0> 
CALL  OPENMS  <.  7  j  JNDX .  5 1 ,  0  > 
CALL  DPENMS  k i  0 »  L I ST  j 1 1 0  0  >  0 > 
ND=NSI2E 

DVERLAYv  iv  u;'  READS  PFQELEM  GEOMETRY  AND  MhTERIhL  PROPERTIES 
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30 

FLN 

31 

FLN 

3£ 

FLN 

33 

I  FLN 

34 

FLN 

35 

FLN 

36 

FLN 

37 

FLN 

38 

FLN 

39 

FLN 

4  0 

FLN 

41 

FLN 

4£ 

FLN 

4? 

FLN 

44 

FLrf 

45 

FLN 

46 

FLN 

47 

FLH 

48 

FLN 

49 

FLN 

50 

FLN 

51 

FLN 

5£ 

FLN 

53 

FLN 

54 

FLH 

55 

FLN 

56 

FLN 

57 

FLN 

58 

FLN 

59 

31 


CALL  OVERLAY  OHRaV^  1  ?  U- f  .HRECALL) 

IF  CMESH.GT.O:.'  CALL  DVERl  AV  •■  SHrDYj  £.  0.  6HRECAL,L> 

CALL  ZERD  aEMPn  l?:* 

nu  15  I  =  i?HFLEriHT 

CALL  WR I Tn  ?  <  i  <  TEMP  < 1 7 ' I > 
15  CDMTIHUE 

FnLLGWlNG  PARAMETEf S  DESCRIBE  TYPE  ilF  ANALYSIS  DESIRED 
HOTENSH  .GT.  0  =  MD-TEM;  IDN  AtlALYSISi.  .  =  0  .  WD   CHECK  FOR  TEMSIDNFLH 
INTER  .GT=  ZERO)  =  INTERACT I DN?  lrrrER=0=  NQ  INTERACT  I GN 
NSTEP  =  NO.  DF  INCREMENTS  PER  LAYER 
HI  =  STAFTIHG  HEIGHT  DF  THIS  LAYER 
H£  =  ENLINQ  HEIGHT  DF  THIS  LAYER 
KN  =  NnRMAL  STIFFNESS  FDR  INTERACTIDN 
C  KS  =  SHEAR  STIFFHESS  FDR  INTERACTION 

C         NITER  =  HD.  DF  ITERATION  SPECIFIED  FOR  ND-TENSIDN  ANALYSIS 
C 

£0  READ  '::5;fc.5>  NOTENSH?  INTEF  >  NSTEP;^  HI  <  H£?  KN;  KSs  NITER 
NCYCLE-i 

WRITE  '.:6;70>  IHTER<NSTEPi.Hl5H£OiOTENSN 
£5  nD=NSIZE 
HCuUNT=  0 

C  FDRH  STRUCTURAL  STIFFNESS  MATRIX  AND  LOAD  VECTOR 

C 

CALL  OVERLAY  <  3HRGY !  3<  0 :■  t  HRECALL  > 

WRITE  ';6!?5> 
C 

C  MODIFY  FOR  BOUNDARY  CONDITIONS 

C  SOLVE  FOR  DISPLACEMENT 

C 

CALL  OVERLAY  ■:.  3HR0 Y  ?  6 ;  U  ;■  fcHRECALL  > 
C 

C  FIND  STRESSES  AND  SIPAIHS.  AND  PRINT  RESULTS 

C 

30  CALL  OVERLAY  (  SHROY- 4?  0?  cHRECALL-) 

CALL  OVERLAY  (3HR0Y<  5- ijj  6HRECALI  > 

IF  (IFLAG,LE.O:>  GO  TO  35 
C 

C  RECYCLE  H!R  NO-TENSION 

C 

MCnUNT=NCDUNr-.-l 

IF  CNCDUHT.GT.NlTEft;'  GO  TO  50 

CALL  OVERLAY  <  JHRQY -  6 ?  0^  6HRFCALL ;:■ 

GO  TO  Sij 
C 

C  RECYCLE    PRIjGRhM    FOP    INCRCliENTHL    OR    NON-LINEAR    ANALYSIS 

C 

■35    HC¥CLE=NCYCLE-!-l 

IFLAG=iJ 

IF  CNCYCLE.GT. NSTEP >  GO  TO  40 

GO  TO  £5 
40  ISTEP=ISTEP+1 

IF  <ISTEP.GT.  LAYERS,)  GO  TO  45 

GO  TD  £0 
45  WRITE  <:6j80::' 
C 
C  CLOSE  ALL  MASS  STORAGE  UNITS 


FLN 

60 

FLN 

61 

FLN 

6£ 

FLN 

63 

FLN 

64 

FLN 

65 

FLN 

66 

FLN 

6? 

FLN 

63 

IFLN 

69 

FLN 

7  0 

FLN 

71 

FLN 

7E 

FLN 

73 

FLN 

74 

FLN 

75 

FLN 

"7^_ 

FLN 

77 

FLN 

78 

FLN 

79 

FLN 

3  0 

FLN 

81 

FLN 

S£ 

FLtH 

P,~i 

FLN 

34 

FLN 

S5 

FLN 

86 

FLH 

87 

FLN 

88 

FLN 

89 

FLN 

90 

FLN 

91 

FLN 

9£ 

FLN 

9:": 

FLN 

94 

FLN 

95 

FLH 

96 

FLN 

97 

FLN 

98 

FLN 

99 

FLN 

1  0  0 

FLN 

101 

FLN 

10c 

FLN 

103 

FLN 

1  04 

FLN 

ins 

FLN 

1  06 

FLN 

107 

FLN 

1  08 

FLN 

1  09 

FLN 

110 

FLN 

111 

FLN 

11£ 

FLN 

113 

FLN 

114 

FLN 

115 

FLN 

116 

FLN 

117 

FLN 

118 

FLN 

119 

FLH 

137 

FLM 

13S 

FLN 

139 

I  FLM 

140 

FLN 

141 

FLfi 

14£ 

FLtl 

-143 

32 


CALL  CLOSEMS  (1>  FLH 

CALL  CLDSEfiS  <3y  FLN 

CALL  CLuSEn;  ';4;'  FLH 

CALL  CLDSEnS  <?;'  FLH 

CALL  CLDSEMS  a  HO  FLH 

FLH 

STOP  FLH 

50  WRITE  (t.<i35)  NCDUNT  FLN 

STDP  FLN 

FLN 
55  FGRMAT  <13AS5lE:>  FLH 

60  FCRMAT  <£0:-:f  13A6>  FLN 

65  FDRNAT  <3I£5  4X=  SFIO.  U;£El  0.  0,  15  :■  FLN 

7  0  FORMAT  (. -■■!  OK s  1 7HLAYER  I NFDRNA  !  i  DN ;  .•■■  1  OX  <  1  4H I HTERACT 1 DN  =  .  1 5--- 1  OX «  FLN 
i  35HNU.  DP  INCREMENTS  FDH:  THIS  LAYER  =  !  IS--!  OX:.  3£HSTARTIHG  HEIGHTFLN 
£  DP  THIS  LA'fER  =  j  PI  0.  £■•' i  OX?  S^HFIHISHING  HEIGHT  DP  THIS  LAYER  =  ?  FLH 
3FiO=£-'-luX?  i4H  ND-TENSIDii  =  .I5--10X-  30HHD,  DP  ITERATIOHS  SPECIPIEFLH 
411  =  jIS.-.'-::' 
75  FORMAT  (lOX?  ESHDyERLAY  ■ 3? U>  CnNPLETEDs 
SO  FDRMAT  (.vv. -SX?  i5H  EHD  UP  PRDLLEM) 

•35  FORMAT  ■:.■'■--- 1  OX <  45HND-TEHSIDN  ITERATIDH  liOES  NRT  CONVERGE  AFTER 
15?  i3H    rrERHTIuNS-.-nlO;;!  £0HEXECUTIDH  TERMINATEIO 

END 

SUBRQUriNE  ZERO  CAsH^  ■  ZRD 

2RD 

THIS   SUHKLiiiTIiiE   GENERATES   A   NULL  VECTOR  ZRO 

ZRD 
D I MENS  I OH    A<N> 
no   5    1=1  J N 
Aa>  =  0=0 
5   CONTINUE 
RETURN 

EHD 

SUBROUTINE  NULLMAT  (A-M'ti) 

THIS  ROUTINE  GENERATES  A  NUl.L  MATRIX  ?  A  ?  OF  SIZE  <M  X  M 

D I  MENS  I  OH  A':MhN> 
DO  5  1=1 ?M 
no  5  .J=i?N 

A<i?  •j:"'  =  o.  0 

5  CONTINUE 
RETURN 

EHD 
aVERLAY<PDY<  i  <  0> 

nVLRLRVcRDY, 1,  o> 

PROGRAM  INPUT i                                                INi  3 

COMMON  NNODES?  NELEMNTi  NDL'F-  MEAHD-  HIi-  r'T3»  I  STOP  i  HCrCLEj  LAYERS)  ISTEPj  INI  4 

lNSTEPiiNTi£<ETAjNTii'HT£«NDTENSN)  iPLAG?  HSIZE<  NCDfiEOSO!)  s  XCSSO  !  Yf  550IN1  5 

£>j  JNDX';5i>i  ANLSIS)  IX<8;  £50)  ?  AREAA(£5u>  ?  INDX<£50>)  INDEX(:e5Ct>  ,  GAMA';£IN1  6 

35>jZAI'::3>                                                     IHl  7 

COMMON  .'l--  Fi  NUEi  RADIUS- XCEH-YCEN)  EI  ^KNjKS)  HI '  HS' IHTER            IHl  o 

COMMON  '-'2--'  D                                                  IHl  '? 

COMMDH  ■4--  MuDU'iSb;;".  DELTA                                        IHl  10 

D I  MENS  I DH  L  ■;  1  0  :>  <  HUE  <  1  0  "> .  D  ( 1  0 «  1  0  >                                I H 1  1 1 


ZRD 

8 

ZRD 

9 

ZRD 

10 

ZRD 

11 

ZRD 

1£ 

ZRD 

13 

ZRD 

14 

r<UL 

£ 

NUL 

"I 

NUL 

4 

NUL 

5 

NUL 

fc 

NUL 

7 

NUL 

HUL 

9 

NUL 

10 

NUL 

11 

riUL 

1£ 

i-iUL 

13 

IHl 

1 

IHl 

£ 

33 


DIMENSinN  IVT'EfiU) 

D I MENS  KIN  NL(4> 

DATA  MAKHP  ?  M AXEL  -  MAXMAT-'  35  U  ?  S5  0 » £5-^ 

3ATA  PLSTRS?  PLSTRN.-'bHPLSTRS?  6HFLSTRH-' 

REAL  K2ERQ 

REAL  KNjkS 

REAL  HUE 

NT3=i:i 

NDuF=0 

nEANri=L! 

NT 1=0  INI 

NT£=n  INI 

HT1£=0  INI 

INI 
READ  PROBLEM  STATEMENT  AND  DTKER  PARAMETERS  INI 

INI 
READ  <:5:?5::'  NNDDESj  NELEMHT?  MATRIAL:.  NPRSRj  LAYERS?  ISTQP  INI 

IF  (NNniiES.LL.nArtNP)  Gu  TD  5  INI 

MRITE  <6?S!J>  NNDDES  INI 

5  IF  (HELEMNT.LE.MAXEL)  GD  TO  10  INI 

LiRITE  ':e,e5j    NELEMNT  IHl 

INI 

PRINT  LAYOUT  DF  THE  PRDELEM  INI 

INI 
iO  URITE  <65  9U>  NNQDES-NELE  INT-MATRIALnNPRSRsLAYERS  INI 

IF  CMATRIAL.LE.MhXMAT)  G'J  TQ  ib  INI 

WRITE  Cb!'?^)  M ATRIAL 
15  CONTINUE 

READ  MhTERIAL  TYPES j  PROPERTIES  AND  DESCRIPTION 

WRITE  <6<105;' 

CALL  ZERO  n3AMA<£5> 

DD  £0  I=l!MATRIAL 

READ  ■;  5 .  1  0  0  >  NTYPE  :■  GAi1A  <  NT YPE ;'  <  TYPE 

WRITE  <6<ilU>  NTYPE?  GnMA'.NTV'^E).  TYPE 
£  0  COMT I  NiJE 

READ  <5?li5)  ANLSIS'DELTi-i 
Du  £5  I=l!MATRlAL 

CALL  INPUTS  U> 

IF  CANLSIS.  EC!.  PLSTRS >  CALL  ^''.NSTRS  CI) 

IF  CANLSIS.EQ.PLSTRN:;'  CALL  i  ARAMTR  ( I  > 

IF  TSTDP  .GT.  Ji  RD;  Nl'N-LINEAR  PROPERTY  SPECIFIED 

INI 

IF  CISTOP.GT.  O:)  CALL  ::PLINE  ;i>  INI 

£5  CONTINUE  INI 

INI 

INI 
CALL  ZERO  CXj NNDDES)  INI 

CALL  ZERO  CYiiNNODES:)  INI 

CALL  ZERO  (NCEiDEiiSSO:)  I  Ml 

INI 

READ  HuDAL  POINT  DATA?  LlJDE!.  COORDINATES:.  LOADS >  ETC.       INI 

NREAD  =  TOThL  fiO.  OF  MODAL  DATA  CARDS  SUPPLIED  INI 


INI 

1£ 

INI 

13 

INI 

14 

TN1 

1  5 

INI 

16 

INI 

1? 

INI 

18 

INI 

13 

INI 

£0 

INI 

£1 

INI 

4  0 

INI 

41 

INI 

4£ 

INI 

43 

INI 

44 

INI 

45 

INI 

46 

INI 

47 

IHl 

4S 

INI 

49 

INI 

50 

INI 

51 

INI 

5£ 

INI 

53 

INI 

54 

INI 

55 

INI 

56 

INI 

5? 

INI 

jy 

INI 

59 

34 


FfHCTOR  =  UNIT  CaHVE-'SIDN  P^lCTDP  FOR  HODhL  CDDRDINATE  DATA     INI 

INI 

READ  <5'1.£0)  NREAIifFACTDw  INI 

DD  3u  I  =  i-NRLAri  INI 

READ  <: 5  n  1  eS  >  NN n  NCDDE  •;  HH )  j  M  v  NN  ;■  ?  V  ( NN  >  I  \\  1 

IF  (  NCDDE  ^hN>.GT  =  3>  Hr:S=NT3tl  INI 

X  C  NN  J  =H  ■:.  NN  >  "FACTOR  I N 1 

Y(:NN>=Y<Nr^>*FACT^R  INI 

30  CaHTIHiJE  INI 

:iS«j»:**«>  ;.,*<•*******«*****  ************** 

READ  ELEMENT  DATA-  TYPE-  MATERIAL  INI 

CALL  NULLMAT  < IHt 3? NELEM fT> 
CALL  2ERD  (AREAA? NELEHNT  ■ 
Du  33  I=i)NELEMNT 

READ  tS'lSCO  lELn  (Ix(  J?  IEL>?  J=i?8::' 

NT=IX<:75lEL> 

IF  aHT.EQ.l)  NTl=NTl+i 

IF  (NT.EQ.E)  NT£=NT£+1 
35  CDNTINUE 

NT1£=NT1+HT£ 

WRITE  NDDAL  PniNT  INFnRMATIDNS 

WRITE    '"Ci  135:j 
DD    40    I=^l<HHnriES 

WR I TE    ■;  6 !  1  '*  0  >    I  ?  NCDDE    I  >  ?  H  <  i  ,■■  ?  Y  <  I '? 
40   CDNTINUE 

NRI7L   ELEMCNT    INFDRHATIDNS 

URITE    ■::&j145) 

DD  45  I  =  i!HELEriNT 

WRITE  Cb'lSO::'  I  j  aX'::  J»  I )  ?  J=l  ■•  8> 
45  CDNTINUE 

READ  ( 5  •  1 6  0  >    2A I  <  1  >  ?  Z  A I  < .':  >  .  ZA I  <  3  ':> 
WRITE  <t<  155-'  ZAKi>?ZAI' £)'ZAI '^S;' 

Dl.  TERM  I  HE  BAND  WIDTH 

DETERMINE  NDDAL  FDRCIS  DUE  TD  SELF  WEIGHT.  FDR  THIS  PURPOSE 

WR I TE  ( 6  J  1 65  > 

DD  65  II=i!NELEMNT 

HT=I>^f;7rII> 

GD   TD    (ISO?  55?  60?  60'^  60    n    NT 
50  Ii  =  IxaiiII) 

I£=IJiC2»II> 

IDF-I1-I8 

NBAND=  <  I rB'< '■  I DF >  + 1  > *3 

1 F  <  NBAND :,  GT .  MBAND  >  ME AND=NE AND 

GD    TD   65 

55        ii  =  ixa.ii>    , 

I£=IX<£-II> 

I3=IX':.3?  Ii:; 

I4=IX':;4!lI> 

HLa>=3*r[-£ 

NL<£>=3*I£-£ 

NL(3>=£*I3^-1+NT3 


I  HI 

R4 

INI 

85 

INI 

86 

INI 

37 

INI 

83 

INI 

8'^ 

INI 

30 

INI 

91 

INI 

92 

INI 

oo 

INI 

34 

INI 

35 

INI 

96 

INI 

97 

INI 

98 

INI 

99 

INI 

1  0  n 

INI 

T  N4  1 

101 

1  no 

INI 

1  '.'c 
1  03 

INI 

104 

INI 

105 

INI 

1  06 

INI 

107 

INI 

1  08 

IHl 

1  n? 

INI 

no 

INI 

111 

INl 

HE 

INI 

113 

INI 

114 

INI 

115 

INI 

116 

INI 

117 

INI 

118 

INI 

119 

IHl 

120 

INI 

lei 

INI 

1££ 

INI 

1£3 

INI 

124 

INI 

1£5 

INI 

126 

INI 

127 

INI 

128 

I  HI 

129 

INI 

130 

INI 

131 

35 


NL<4>=e»iI4-l+HT3 

MAX=MHX0'::HLa:J?NL<£)!ML<:3:)jiiLC4>> 

MIN=M!Hn<HL(l>>HL<e>?:iL'::3;>.HL<4-0 

NEHnri-=MHl-v-f'iIN+£ 

IF    '::NBHND,GT,MI:AHri>    f-ISAHD^H&tiHD 

GD   TD   65 

CALL   AREhS    <TIjNT:> 
CDNTIHUE 
NrinF=S*MHuIiLS+NT3 

FRIiiT   SIZi:    OF    THE   PRDELEM 

WR I TE    ( fa . 1 7 U  >    NDDF  j  MEhMD 

IF    (NT3.LE,  0)    GiJ   TO   70 

kEhD   CENiER.    RhiiIUS,    and   EI    DF   PIPE 


KLHH    (.SjlZr:.. 
WRITE    <6!l8i 


0;:EN.YCEN.RnIiIUS<EI 

;«;CEHjYCEH,r:HDiu.s;.Ei 


riETERMIHE   ETA   FOR    TRIANGULAR    ELEMENTS 


NT1£P1=HT1£+1 

IF    (NTiePi.GT.NLLEIMNT)    GO    TD    70 

II  =  IKa<NTl£Pl;' 

ie=IK<£!NTl£Pi:) 

;-';i=X';ii> 
Ti=S'';ii::' 
j^£=x<ie> 

V£=Y<I£> 

SPAN=SQRT  C  ■::  X 1  -X£  >  «*£+  ■;  V 1  -  Y£  > **£  > 

SN=0.  5«SPAM-'RAriIiJS 

TH=ASIN<SN> 

THETA=£.'*TH 

R£=kADIUS**£ 

Ri=0.5«THETH+R£ 

A£=0.5*R£*SIN(THETA> 

A4=A1-A£ 

ETa=A4 ■ AREAA  ^  NT  i  £P  i  > 

ETA=0. 0 

70  CDNTINUE 

URITE  (b-lSS,' 
RETURri 

75  FaRMAT  (£15^  £13,  £Xi.  I£?  £Xn  IE?  IX- 1 1 1.  IXs  I£::' 

80  FDRHAT  CSX?  4£HNa=  OF   NDPES  EXCEEDS  LINIT<=900>?  NNnDES=  .  I5:i 
S5  FORMAT  CSXj  4t.HNn.  DF  ELtNENTS  EXCEEDS  LIMITOSOO;  ?  NELEMNT=  n  I5> 
90  FORMAT  a  OX.  i9HPR0BLEN  DESCRIPT IDN?  .-•.•-•SX?  19HNn,  OF  NODE  POINTS=;  I 
i5.''5X!.  ibHND.  DF  ELEMENTS=hI5.--5Xj  17HND.  OF  MATERIALS=-- IS-'SX-  EbHND 
£.  UF  BOUNDARY  PRE3SURES=<  IS-'SX-  £9HND.  DF  CONSTRUCT  I  ON  LAYERS  =  ?I 
35.'-.-'> 
FORMAT  CSX?  4?riND.  OF  MATIRIALS  EXCEEDS  LIMIT  •=  £5>?  MaTRIAL=!IS> 


i  0  0  FORMAT  <  1 5  J  F  1  0  .Oil  0 A6 ;' 

105  FORMAT  (.-'.•-•^'■iOX?  £iHMATER3  iL  DESCRIPTION:; 

IWEIGHT?  i  OX?  i  IHPESCRIPTID  H?  -';:■ 
i  i  0  FORMAT  ■;  5X  ?  I  £  ?  8X  ?  FS .  4  ?  1  OX  ■  1  0 Ab ) 
115  FORMAT  CAfe.FS.  O:' 
1£0  FORMAT  (I5?F:!0.  0> 


5HND. 


11  HUN IT 


INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 

INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 
INI 


132 
133 
134 
135 
136 
137 
138 
139 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 

1 5  0 
151 
15£ 
153 
154 

156 
157 
158. 
159 

1 6  0 
161 
162 
163 
164 
165 
166 
167 
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125  FDRMHT  (  15- Si^Me?  Fl  0.  O,  £0K- Fl  C.  n>  INI  i'. 

130  FDRMAT  <  13?  aX?  SOIS?  5X:>  :■  4X?  1 1  .SX?  I£?  £F5.  0- 5Xi.  £F5.  CO  IMl   193 

135  FORMAT  C-'.-'^'l  OX?  16HHDDHL  PDIHT  DATA!  ■■-■■■SX,  SHIiaDE  NO.  i.  £X-  9HriIRECirU  I'i 
1TIDNI.3K!  iOHFIXED-'-FREE^S  ;-  1  IHCD-DRDIHATE1.5X5  luHNDDAL  LOHDneX-  9IH1  V\ 
£HriIRECTIuH'.3XH  1  OHFIXED.-FREE.  SXj  1  IHCD-nRDINHTE, '5Xj  1  OHNDnAL  LOAD  j  INI  V. 
3-'->  I  hi 

140  FaRMAT  (7Xi.  I4?yXi.   IHX;  1  OX?  II  j  1 IX?  F8.  £1.  £4Xi.   IHY:.  ££Xi  F8.  £::•         INI 
145  FORMAT  C'.--.--i  0X>  ISHELEMEriT  BA7A  INPUT,  ■-.  IX,  IIHELEMENT  HD.'EXj  591-IINl 
i  _p_   _Q_   -.R._   .-s-   _PQ-   -QR-   _RP-   -RS-   TYPE   MATERIAL) --■-O    INI 
150  FDRMAT  CaX,  15?  6X5  3a3,  eX>  ,  SKs  SKM,  ex:> ,  6X?  £X«  I£,5Xj  I3?4X?F6.  i.c:X,F6INl 
l.£5  1XjF4.  0,  iXiiF4.  n>  IHl 

155  FORMAT  auX?  yRZAICn  =  jFS.f'SX,  9HZAI<£>  =  ?F5.£?5X,  3HZAI<:3:ilHl 
1  =  ,F5.  £.-':>  IHl 

IfcO  FORMAT  <3F10„n;:.  TNI 

1 65  FORMAT  •;  ■''  i  O;-: ,  £  UHELEMENT  1  HFDRhl  AT  I  DNS , .- 1 X ,  ?HEL.E .  HO .  ,  9X ,  6HZA I  <  1 1 N 1 
i>,9X,  6HZAI<e>,9X-  6HZAI(:3>,  lOX,  4HhREA-5X«  £4HX-CDDRDINATE- Y-CINl 
tORn  I  NATE -.■■■■>  INI 

170  FORMAT  ■:>----l  OX,  1£HPRDBLCM  SIZE, --■--■£5X,  £?HND.  OF  DEGREES  OF  FREEDDINl 
IM  =,I5.'-35X5  17HSEM I -BAND-WIDTH  =,  I5.->  INI 

175  FORMAT  (3F1  0.  0,  El  0.  O!)  INI 

ISO  FORMAT  (.'lux,  lyHPIPE  SPECIFIChTIDHS,.'1  OX,  IIHX-CENTER  =  ,Fb.£,3X,INl 
1  liHV-CENTER  =  ,  F6.  £,  •■•1  OX,  14HPIPE  RADIUS  =  ,  Fb.  £,  .--1  OX,  15HPIPE  ST  INI 
£IFFHESS  jEi£,£/>  INI 

185  FORMAT  (IX,  £3HDVERLAt  a,0>  COMPLETED)  INI 

INI   £16 
END  INI   £1? 

SUERDUTINE  INPiJTP  (ID  IN£    £ 

COMMON  ---l-'-  E,NUE,RADIUS,,-,CEH,YCEN,EI,KN,KS,Hi?H£,INTER  IN£     3 

D I  MENS  ION  Ed  0  >  ,  NLiE  <  1  0  >  I  N£    4 

REAL  KN,KS  IN£    5 

I  N£     6 

?;**:^fl):-**ft:"ft*ifl>;fl:.^f^;n':i'i':if'i:*:i^^  ' :  lift  *  *  ■*  .^P,  '1^   *!*:  fl^  ift  i*  *  ^i^  *;  i*^^fc  :*.*;*  fl>  :^  :*;  .^>;  ;♦;*  Ift  .*:  ^P  ;**;  :^^  *  *  I  Np  7 

THIS  ROUTINE  READS  DATA  AND  NtCESSARY  VARIATIONAL  PARAMETERS  FOR  IN£    8 

MATERIAL  PROPERTIES  ANO  CONSTITUTIVE  RELATIONSHIP.         INE     9 

THIS  RDUTINF  IS  DIRECTLY  CONNECTED  TO  *RARAMTR  f-    THROUGH  COMMON  INE    10 

REAL  HUE                                                      ni£  IE 

INE  13 

READ  E  ,  HUE  I'OR  FIXED  PROPERTY                 INE  14 

INE  15 

READ  '::5,5:i  ECIi::'                                              INE  16 

READ  (SjIO)  NUEdl)                                            INE  17 

tJRITE  <6,15)  lI.EaDnNULan                                   INE  IS 

RETURN                                                       INE  19 

INE  £0 

5  FORMAT  (luES.U)                                               INS  £1 

10  FORMAT  (F5. 0>                                                 IN£  ££ 

15  FORMAT  C /'--•,•',-•  i  OX,  £?HPROPrftTIES  FDR  MATERIAL  NO  ,  15,  .'--SX,  15HYDUNGSINE  £3 

1  MDDULUS=,1X,Eil.4.---5X,  15HPDISS0HS  RATIO=,  F5.  E)                  INE  £4 

IN£  £5 

END                                                             INE  £6 

SUBROUTINE  PLNSTRS  (K)                                         PES  £ 

COMMON  .'-l--  E?  HUE,  RADIUS,  ::CEN,YCi:N,  EI,  FN,  KS,  HI,  HE,  INTER            PES  3 

COMMON  --■£.■'  D                                                  PLS  4 

REAL  HUEjNUEK                                                 PLS  5 

REAL  KHjKS                                                       PLS  6 

DIMENSION  Da  0,  10  >,  ECIO:,  NUEOOj                                PLS  7 

THIS  ROUTINE  13  FOR  PLANE  STRESS  PLS    9 

PLS    1 0 
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HUEf 

•:>HIJE 

Ky 

C=E 

;k'. 

..'■  ( 

.  -hit 

ek: 

i^n 

..IE 

U'iK 

8; 

=  0 

0 

IKK 

7'. 

=ri 

K  J  s : 

IKK 

fcV 

=ri 

Kj?: 

D'::k 

3!. 

=ri 

K?6: 

DCK 

5: 

=c 

iiCK 

1> 

=ri 

Kjs:. 

TK.K 

4; 

=C*NUEK 

D'::K 

£; 

=ri 

K54:; 

IKK 

9; 

=  0. 

5*C^^ 

<  i 

1M 

D'::K 

n; 

>=HUEK 

RET! 

jpr 

PLS  1 1 

FLS  IE 

PLS  13 

PLS  1 4 

PLS  15 

PLS  16 

PLS  1 7 

PLS  18 

PLS  19 

PLS  £0 

PLS  21 

PLS  ££ 

PLS  £3 

PLS  £4 

END                                                                                                                                                                PLS  £5 

SUBRDUTINE    PhRrMTR    (K>                                                                                                                PLM  £ 

CDPiMari   .'-l.--   Ej  NUE?  RADIUS?  SCENjYCEN,  EI  .KH?KSi.  HI  <H£.  INTER                             PLfl  3 

CuriMDN   .'■•£■'■    D                                                                                                                                         PLti  4 

REAL   KtisKS                                                                                                                                      PLM  5 

PLM  6 

THIS  RDiJTIHE    IS   FDR   PLAhE    STRAIN                                            PLH  7 

PLN  S 

PLIi  9 

J?    EaU'!    hUEUO;'                                        ■                                          PLN  10 

F1.N  1 1 

J'--<(l.+liLiEf<>*<t.--e.*HIJEK-J>                                                           PLN  12 

PLN  13 

PLN  14 

-NUEK>                                                                                                              PLN  15 

PLM  16 

SNUEKX'-Cc.  *a.-NllEfO:5                  ■              .   .                            .               PLN  17 

PLN  18 

PLN  19 

PLN  £0 

PLN  £1 

PLM  ££ 

PLN  £3 

PLN  £4 

ENJ.I                                                                                                                                                                FUi  £5 

SUBRDUTINE   AREAS    ■::MN-NT>                                                                                                     ARE  £ 

COMMON    MNDHES,  NELEMHTj  MDNf?  MEAND:.  NDj  NT3.-  ISTDP-  NCYCLEj  LAYERS?  ISTEP?  ARE  3 

iNSTEP'NTlPjETAsNTl'NTE'NriTENSN?  IFLAG?  NSIZE?  NCniiEOSO  5  K(:550>  5  YC550ARE  4 

£;■.  JND;-^(:51>?ANLSIS!  IX<S!£;:iJ>!AREriA<£5u>?  INDX(:£5D>!  IMDEXCESCD  ,  GAMA<£ARE  5 

35)?ZAI<3;'                                                                                                                                                 hRE  6 

ARE  ? 

ARE  8 

>;*******¥*;*****«*  w**w^(t  **;«**:**>»■  ******!***  9 

RDUTIHE  CALCULnTES  ELEMENT  AREA-  AND  SEMI  BAND  WIDTH     ARE  10 

ARE  1 1 

ARE  1 E 

ARE  13 

ARE  14 

ARE  1 5 

ARE  1 6 

NK                                              ARE  17 

ARE  18 

:.N..  nN>+riT3-l                   ■                                                                                             ARE  19 

ARE  £  0 

'nN>-e                                                                                                                          ARE  £1 

'Mn:>-£                                                                                      are  ££ 


REAL   NUE 

HUEK 

DIMEH: 

.  I  DM    D •  1 U  <  1 

nuek:=nue 

K> 

C=fE(f 

•;>* 

I . -NUEh 

D'iKjS 

=C 

Da;:?! 

=D 

K?5:) 

D'::K!4 

=C«NUEK.-a 

DO<.£ 

=D 

K5  4> 

DCKjA 

=C*<i.-£=* 

D<Ki.8 

=  0. 

0 

L<Kj? 

rrri, 

K?  8> 

DCKjb 

=D' 

i< :.  7  > 

DCKjS 

=-^D' 

K?&> 

D(:K?ii 

i::'=HUEK 

RETURf 

\ 

nif 

1EN 

SIDN 

LM 

J              *:**i*j> 

****!*«* 

THIS 

RD 

AREA-- 

0 .  0 

1  = 

H< 

1  ?  MN 

J= 

IK< 

£?MN 

K= 

X< 

3  5  MM 

MK= 

=NT 

-£ 

GO 

TD 

<5ii 

5-, 

5   DD 

IC 

H=l 

6 

iO   LM 

N) 

=£*i:- 

:(.N 

GD 

TD 

3  0 

15   LM 

i> 

=3*1,' 

Iv  i 

LM 

£;■ 

=3*1) 

;<£ 
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LM<3>=e*IX<3.hrO-l+NT:3                                           ARE  £3 

Lri(4>=3*IX<45MH>-£                                               ARE:  £4 

L|-'i<5>=£*i:«;(5?MN::'-l+MT3                                         ARE  £5 

LM'::b::'=e*IK':;fc.i.MN>-l+NT3                                           ARE  £6 

GD  TD  30                                                     ARE  £7 

£0  DD  £5  H=l!b                                                  ARE  £3 

£5  LM';N)=£*IX<M:.MrO-l+MT3                                         ARE  39 

LM';:3>=3*IX(3i.M|-0-£                                               ARE  30 

30  AREA=AESi:.5*<K(J::'*Ya(>--T<J>**Xa')+Xa>*Or(J>---Y<K>J+Ya>*<X'::K>-XU>>ARE  31 

i  >  ::■                                                              ARE  32 

riAX=riAxo<:LM(:  i  >  ?  lm<:£>  <  lM'::3>  ?  lm>::4>  ?  lmcs::-  ?  Ltice:'  ::■               are  33 

M I H=M I N  0  <  LM  a  )  <  LM  (£::'-  LM <  3  >  j  LM ':  A),  LM  <  5  >  5  LM  C  6 !>  >                     ARE  34 

NWiriTH=MA:<-MIH+£               _                                 ARE  35 

IF  (HWIDTH.Gi  .f'ihiANIi)  r1BpiiIi=Mi.-rj  D I H                               ARE  36 

IF  ';HUniTH.QT„MD>  WRITE  Oi!40>  HWIDTHi.MN                          ARE  37 

IF  U-iREh.GT.  0.  0>  GD  Tu  35                                      ARE  3S 

URITE  <6.45>  HN                                                  ARE  39 

RETURN                                                          ARE  40 

ARE  41 

S5  xx=ZAia)*xa>+ZHi(£)*x<j>+ZAi>:3>*Xf:K:)                  are  42 

VY=ZAia)*Va>+ZAI<£)*Y(..!)+ZAI':3:i*Y'::K>                                                                    ARE  43 

AREHH<nN)=AREA                                                                                                                                    ARE  44 

iJR  ITE    ( &  J  5  0  >    MH  5  (  ZA  I(  I  :^  ?  I  =  1 !  3 )  ■■  AREA  ?  y.K  <  YY                                                            ARE  45 

RETURN                                                                                                                                               ARE  46 

ARE  47 

4  0  ruRMHT  (3X5  3£HBANr!WIDTH  EXCEG.l'S  LIMIT 5  MEArili-^  ,15.5X5  17HIH  ELEMEARE  48 

IHT  HD.  =  I.  15.:.                                                       ARE  49 

45  FORMAT  <5Xj  35HHEGATIVE  OR  ZERO  HREAj  ELEMENT  Nn„=.I5>            ARE  50 

5  0  FORMAT  C £X .15?  5X .  4E 1 4 .  £ .  .i'F  1  £ .  £ >                                  ARE  5 1 

ARE  5£ 

END                                                             ARE  53 

SUBRDUTINE  SPLINE  (I MAT)                                       SPL  £ 

CuMMDH  ■•■'4--'  MuriU(:£6>!liELT-i                                        SPL  3 

niMENSinfi  XP';7>.  tPdlj?-?  SEP<il.7::'.  PSNR(li-7>!  SPSNRai5  7>     SPL  4 

fl  I  MENS  I  Ufi  PROP'  (319  > .  SRO  ;  1  i  >                                     SPL  5 

SPL.  t 

SPLINE  FITTING  FDR  TAN,  Mnii.  TAN.  PDISSSDN  RATIO  VSi.  SIGMADCTA   SPL  7 

NP=  NO.  OF  Di4T,J|  POINTS  FDR  STRESS                SPL  8 

PHI==  FRICT'ON  ANi.i.E  FDR  SOIL  (DrGREES)              SPL  9 

PSI=  RATIO  OF  (:SIG;iA£>.--(:  :  [GMA1+SIGMA3;:'  FOR  PLANE  STRAIN      SPL  10 

AHISD=  DEGREE  OF  AH  i  ::nTRDPY  OF  SOIL  STRENGTH           SPL  11 

SPL  1  £ 

JDIL  AND  PIPE     SPL  13 

'ACTDR           SPL  14 

SPL  1 5 

SPL  1 6 

SPL  17 

SPL  1 8 

SPL  19 

XP=  DCThHI  AURAL  NORMAL  STRESS  VECTOR               SPL  £0 

(  P=  TAfHiENT  MODULUS                       SPL  £1 

PSNR=  TANGEi^T  PDISSDHS  RATIO                   SPL  ££ 

i:PL  £3 

READ    <5.65>    (XPd  >  ?  1  =  1 5  Nf  O                                                                                                      SPL  £4 

DD   5    1  =  1  ..HP                                                                                                                                           SPL  £5 

XP  < I >  =F ACTDR*XP ( I >                                               SPL  £6 

DD  15  HD=i. 11                                                 SPL  £7 

READ  (:5.?0>  (EPtNQ.  J>.PSNR(HD?.J::'.  J=i.HP>                       SPL  £8 

SRD ( NO >  =  ■:;  11. -FLOAT (NO ;■::'. ■■•10=                                    SPL  £9 

DD  10  J=1.HP                                                  SPL  30 


FDR  ISllTPDPIi 

MATERIAL  ANISD=1 

delth=  hngle  of  FRICIDN  (. 

OEGREES  >  BETWEEN 

fhctqr=  strength  Pi 

OPERTY  CONVERSION 

REhE  (S.55>  HP 5  PS 1 5 PH I ! an ISO s 01 

J.TA. FACTOR 

IF  (ANIED.Ll.  0.  0:.  ANISO=  .  0 

FRATIO=FhILUPE(PSI.PHI> 
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iO    EP(ND,J>=f"HCT!JP*EPai[;:5  J::' 

CUBIC  SPLINE  FITTIHC  FDR  T;>riGEHT  hiniiULUS  hHD  PDISSDNS  RAT  ID 

CALL  SPLFIT  CHP, ND- XP- EPjSEP> 

CALL  SPLFIT  (HP:.  ND?  XP-PSNR.  SPiHR> 
15  CDHTINUE 

WRITE  <6!?5::'  IriAT 

URITE  <&f60?    FHI<PSI'iFRftTia5r!E!.TA5ANISO 

WRITE  <b)3i:i>  (XPvi;.'i.  1  =  1,  MP) 

WRITE    CbjSS)    <SRDa>.  1  =  1' 11) 

nn  £0  1=1 jNp 

£A   WRITE  >;6j90>  <EP<J,  I  > »  J=  i  >  11  > 

WRITE  (6.95:'  <SRaa.),  1  =  1' li;:' 

DD   £5  I  =  1 5  NF' 

£5   URITE  ';6j9n>  ':.PSNR<  Jj  I  >  j  J=l  n  1 1  > 

xi'DRE   MATERIAL    PPDPEFTIFx   DH   MASS  STDRAGE   UHIT 

pRDPa:j=ANisn 

PROP<£>=riELTA 
PRnP''3>=FRATin 

riD=3 

Du   30    l=l?7 
MQ=Ma+ 1 
30   PRQP':MD>=JiP(I> 
DD   35    I  =  1 J  1  i 
DD   35    J=l^? 
hiu=MD+l 
35   PRDP<rlD>=EP<Ii..J> 
DD    40    1=1, 11 
DD    40    J=1j7 
r1D=MD+ 1 
4  0    PRDPt:MD>=SEPa?J> 
DD   45    1=1  J  11 
DD   45   J=1j7 
f'iu=nQ+ 1 
45   PPuP'::f1u>=PSHR<IjJ;' 
DD   50    1  =  1.'  11 
DD   50   .J=ij? 
Mn=nD  -I-  i 

50  pRap';Ma>=spswR(i?  J) 

PROP  (.  3 1 9  ':>  =FLDAT  <  fiP  > 

CALL  WRIThS  (4i.pRap5  319<  tMAT) 

RETUPM 

55  FDRnAT  ■:.  1 5 1  5F 1  0 .  0  > 

60  FORMAT  ■■iXf    E4HFAILLIRE  RiiTIO  FDR  PHI  =  i.F6=£.  1£H   AND  PS  I  =  ?F5. 

1'   SH   IS  =   jF6.3,--10Xj   3HIIEL1A  =  ?F6.£j  ££H   hHISDTRDPY  FACTOR 

£  jF6.£';' 
65  FORMAT  <8F10, 0> 
?!J  FORMAT  '::£F10.  0> 

75  FORMAT  fi-'i  OX?  45HHaN-LIMf:AR  SOIL  PROPERTIES  FLIP  ilATEPIAL  HO  =  -,  It 
80  FORMAT  (.IX,     1  1HSIGMA(:DCT>  =  ?  1£E5=  1  > 
85  FORMAT  OOX?  ££HTAHGENT  HDDULUS  VALUES? -EX!  15HS 

11?  10';6X!F3.  i  :>'■') 
90  FORMAT  Q5X?    ilE9.£> 
95  FORMAT  C-'-'-lOX!  £8HTANGEN1  PaiSSlIM  RATIO  VALUES?/ 

10  =  jF3.i!.10<6K!F3.1>.-> 


SPL 

31 

SPL 

32 

SPL 

oo 

SPL 

34 

SPL 

35 

SPL 

36 

SPL 

37 

S:PL 

oo 

SPL 

39 

SPL. 

4  0 

SPL 

41 

SPL 

4£ 

SPL 

43 

SPL 

44 

SPL 

45 

SPL 

46 

SPL 

47 

S:PL 

48 

SPL 

49 

SPL 

50 

SPL 

51 

SPL 

5£ 

SPL 

53 

SPL 

54 

SPL 

55 

SPL 

56 

SPL 

57 

SPL 

53 

SPL 

59 

SPL 

6  0 

SPl_ 

61 

SPL 

6£ 

SPL 

63 

SPL 

64 

SPL 

65 

SPL 

66 

SPL 

67 

SPL 

6  'i' 

SPL 

69 

SPL 

7  0 

SPL 

71 

SPL 

7£ 

SPL 

r'  I' 

SPL 

74 

i:PL 

75 

SPL 

76 

SPL 

77 

SPL 

78 

SPL 

79 

SPL 

30 

SPL 

81 

S  RATIO  = 

«F3.SPL 

SPL 

SPL 

15HSTRESS 

RAT  I SPL 

:SPL 

SPL 
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EHl!                                                             SPL  91 

SUBRDUTIHE  SPLfl r  (HPM? ^Qj XP? Vr) YDP>                              SPF  2 

niMEHsiDH  xp'-?:'!.  Ypciij?;?  vnr  ii??)?  hc?:?.  Y';?)^  ai<7;'<  ek:?)?  cispf  3 

i(7>j    DI'::7>)    SIiPC?)                                                                                                                   SPF  4 

SPF  5 

CUBIC   SPLINE    FITTING                                                           SPF  6 

SPF  ? 

HPi=HPH-i                                                                                                                                                 SPF  e 

IiQ   5    r1=l<NPl                                                                                                                                         SPF  9 

5  H(:ri>=j<;p<M+i:5-xp';M::'                                spf  in 

SLDPl=FD(:H'::i::'5H(£>i.Tp<NDa>«YP'::NDi>£>?TP<Na;3>>                   SPF  11 

SLnPN=BD <:. H < NPN-e >  5  H < NPN- 1  >  ;  YP < HD ,  NPN-S >  ,  YP < HQ .  NPH- 1  >  -  YPC MD ,  NPfH )  )   SPF  1  d 

riD  10  N=l!NPN                                                    SPF  13 

0  Y':M;'=YP<Nu?M>                                                    SPF  14 

CALL  CDFRIT  <HPN»  XP?  Vj  SLOPl ,  SLDPHi.  AI  ?  BI  ?  CI  <  TiJ  >                    spF  15 

CALL  TRIDGHL  ';NPN,  rI  ,  EI>  CI .  DI «  SDP>                                SPF  16 

DO  15  I  =  l!HPN                                                    S'PF  1? 

5  VDP  ■; ND  5  I  >  =SDP  ■:  I  >                                                 SPF  1 8 

RETURN                                                       SPF  1 9 

SPF  £0 

t-Nii                                                             SPF  £  1 

SUERDUriNE  CljrRiT  CNPNjXfj Yp5SLapi5SLGPN?Hl!EI?CI«DI>             CDF  £ 

EiIMENSIDH  -;p<?:'?  YP':;?)-  ■n(:7>.  EI<7)j  CI<7>!  DI';?)                CDF  3 

CDF  4 

GENEP iTE  SPLINE  COEFFICIENTS                   COF  5 

CDF  6 

Hl<i>  =  i:i.O                                                       CDF  7 

Bia>  =  0-;P<£>-XP<l>>.'3=                                           CDF  S 

CI  i;0=B I  (l  ■)■■'■£,                                                   CDF  9 

iii<i>=cYP<£>-Ypa>:;.-';XP(:£>-xp'::n;.-sLapi                  cdf  lo 

H I  <  NPH ')  =  (;  XP  i  NPH  >  -NIP  i  NPN-  i  ;■ )  .'■6 .                         '           CDF  1 1 

BI<NPN:i=ftI<NPN>*£.                                               CDF  IE 

CI<NPN>=U.O                                                      CDF  13 

D  I<  NPN> =-  <  YP':;  NPH  ;■  -YF'  <  NPN  ■  1  >>  ■■  <  XP  <  NPN  >  -XP  <  NPN- 1  >  >  +SLDPH            CDF  1 4 

Hi=NPN-l                                                         CDF  15 

DD  5  I=£?N1                                                   CDF  16 

Hia;'=<xp<i>-xpa-i;o--6.                          cdf  17 

Eia;'=<:xp<i+i>-xpa-i  :o---3.                         cdf  i8 

CI<:I>  =  <XP'::t  +  1>-XP'::I>  ;■  -6.  ■                                     CDF  19 

Di<i>=crp' !+i>-Ypa>;'  ■<xp<i+i>-xp a ;:■>-(: YP(!>-YP<i-i> >-'<:xpa>-xpcaF  pfi 

1    5I-n>                                                        CDF  £1 

5  CDHi I HUE                                                         COF  £E 

RETURN                                                          CDF  £3 

CDF  £4 

t-Nli                                                             CDF  £5 

SUERDUTINE  TRIDGML  (HPNmsI < BI « CI ? DI ? YDP >                        TRI  £ 

DIf'iENSIDN  Rl';?>5  BK")?  rK?),     i:il(;7>j  VDP<7::'<  UdiV)?  U<10>         TRI  3 

TRI  4 

SDL'E  TRI-iJlAGDNAL  NATRIX                    TRI  5 

TRI  6 

TRI  7 

TRI  S 

TRI  9 

TRI  10 

TRI  11 

TRI  1£ 

TRI  13 

TRI  14 

TRI  15 

TRI  16 


p=Bia> 

Q<i;'=-ci' 

:.1).'T 

ucn^DK:: 

L-:'-'P 

DO   5   K=£! 

.  NPN 

p=Aia 

(>«Q':: 

K- 

-i; 

'+b: 

kk: 

Q':.K)=- 

-C I  <  a 

: ) , 

-p 

U':;k:;=' 

::di<k 

-A] 

:(k: 

)  *U ' 

::k- 

■1 

CQNTINUE 

VDP'::  NPN  >  = 

=U'::NF 

■n: 

N1=HPN-1 

41 


DO    iij   1=1.  HI 

K=rU  +  l-L 

YDP';:K>=Q':K.'«YDP'::K+1>  +  J':;K> 
1  0   CONT I  f iUE 
RETURN 

ENS 

FUHCTIDH  FrifSi<S£?Rt<R£>R3::' 

FDREUARD  niFFEREHCE  FUNCTIDM 

IF  CSl-SE)  5? 10? 5 
5  FD=<RS-R1>.--S1 
RETURH 
i  0  FIi=  ( - 5 .  *R  1  +4 .  *k£'-R3 >  ■■  <  15 .  *S i  > 

RETiJf;:H 

END 

FUHCTION  BPvSi  .>  $■-?  Rl  ?  RiE.'?  R3) 


BHCKUiiPD   Dlf^fERENCE    FUNCTIDH 


II-    fSi-S£>    :dii10?S 

RETURfi 
L  0    BIi=<3.  *R3-4.  *R£+Ri  :).■-■<£.:*:  :l  :> 
RETURN 

END 

FUNCT I DN    DRIJ I  rCT  •:  H  -  B  -  C  ■>  D   E .  F  ^  P  - 

H J=D-C 

B1=D-P 

E£=P-C 

Hl=Bl**3 

Fi£=E£*:*3 

Ti=Bl*H.'<b.*H  i> 

T£=ft£*B.-  Cb.sHJ;' 

T3=';E-A*HJ**a.'6.  ;'*<D-Pj-'HJ 

T4=c:F-BS«HJ**3. -6.  ::'*(P-C>^-'HJ 

URD  INET==T  i  +T2-I-T3+T4 

RETURN 

END 

FUNCTIDH    FAIiJJRtt  RSI  -  PHI  J ) 


DETERN]HE   FrtfUJRE    STRESS   RATIO 


REAL    HPHI 

PHI=PHII*PI/Iy!:i. 

SINPHI=SIH(;pHi> 

NPHI  =  a.+SInPHI>--<i.-SINFHi: 

T 1  =6 .  *HF'H  I  .-■  '■  <  NPH I  + 1 .  >  **£ :.' 

T£=£.*'::PSI**£-PSI  +  1.  > 

T3=SQRT';T£-Ti> 

T4=T3'-a.+PSi::' 

FhILURE=T^ 

RETURH 

END 
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!?■ 

TRI 

IS 

TRI 

19 

TRI 

£0 

TRI 
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TRI 
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FD 

o 
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3 

FD 

4 

FD 

5 

FD 

6 

FD 

7 

FD 
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FD 
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11 

FD 

1£ 

BD 
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BD 

3 

BD 

4 

BD 
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ED 
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BD 
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BD 

3 

ED 

5 

BD 

1  0 

ED 

11 

ED 
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DRD 

£ 
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DRD 

4 
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DRD 

3 

gPD 
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10 

DRD 
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DRD 

1£ 

DRD 
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DRD 
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15 

FR 
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3 
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FR 

5 

FR 

& 
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FR 

R 

FR 
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FR 
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■j  1 
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1£ 
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FR 
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OVERLAY (RDY.; 

j  (.0 

DS'ERLi^ 

PRDGRAM  LAVrii 

T 

CDnr-iuH  hHDDEr 

<HELEMHTmi:ii 

)F?MEAi 

INSTEP.  HTi£,  ETHj  NTl  ■■  NTgi.  N[iTENSN< 

S).  JHDKfSn-ANLSIS.IXcSjPf 

;ij;'i.ARL 

35>7ZhK3::' 

DIMENSIDN  XPaLY<t.>?  YPOL' 

'■  6 ') 

CnLL  PLDTS 

REttB  >::5?£5::'  ;■ 

NAXs YHAX? YC! 

:n 

SF=ii:i./YMAX 

CALL  FACTOR  ■■ 

SF) 

CALL  PLDT  (0. 

n.  u.  Of  -3> 

YPOLy';4;'  =  0.  0 

XPDLyc:4::'=YPDI 

Y(4> 

YpaLY<5;'=i.o 

MPDLY'::5:>=YPui 

Y':  5> 

CALL  AXIS  '  lY 

U5  0=  U'EHM-j- 

-£.XMA>- 

CALL  AMIS  kQ. 

0. 0. Oj£HY-iL 

^jYMAX< 

CALL  PLDT  OJ. 

OjYMAXfS:' 

CALL  PLDT  i.K 

AXsYMAX^e:; 

CALL  PLDT  ':;^f 

AX?  O.  ii!"£> 

DD  £0  I=i!NELrMNT 

riT=IX'::7?i: 

GD  TD  <5,1 

U'lS)?  NT 

ll=lxi£! i; 

i£=iX';i?  i: 

XPDLYa>=; 

'.:  1 1  > 

VPDLY':.  i>=' 

'!  I  i  > 

XPDLY  (?::■  =  : 

<I£> 

YPDLY  (£:■=■ 

'  I£) 

YPDLY' 3)=i 

■  0 

XPDLYi:;3>=-- 

PDLYO;:' 

YPDLY  (:4>  =  ; 

.  Ij 

;«;PDLY<4>=' 

PDLY(4> 

CALL  LINE 

(XPDLY'YPDL' 

.e.i-i; 

GD  TD  £0 

CDHT I NuE 

Il  =  l>i<A,  1 

I£=IK(£.i: 

I3=IM<3?i: 

XPDLY';;4)=:^- 

<  1 1  > 

ftPDLYa)=? 

PDLYc;4> 

YPDLY<:4>='i 

<  I  i ; 

YPDLY<i>=Y 

PDLY<:4> 

XPDLY  <£;■=; 

•■.l£':> 

YPDLY  <:i:::'='i 

'..IQ-? 

XPDLY(:S)=?' 

aS';' 

YPDLY<3;:'='i 

■as;' 

ypdlv<:5>  =  i 

.  Ij 

J:;PDLY(5;'  = 

PDLYCS' 

YPDLY.::6;:'  =  1 

=  Ij 

xpcly(:6>='t 

PDLY<6> 

CALL  LINE 

(XPDLY.YPDL- 

-4;  i-,i. 

y.c=':y.<:.ii-)-i 

x<i£>+xa3>: 

••'■3.  -0. 

VC=iVtIlM 

Ya£;'+vci3;.: 

'■3. 

CALL  NUMEf 

R  ';xc.ycjCi.^ 

.  I. 0. r 

CONTINUE 

i>^C=knAx-'£.-£. 

'i'C=Yr'iAx+£. 

iRDY' 


U.J 


DjNE<NT3i.  I  STOP.  HCYCLE.  LAYERS.  I  STEP. 
lFLAG>NSIZEjNCDDE<550;>i.K';550>.Y<55i: 
AA<£5u:).INIiX':£50>.  INDEX (:S5n::.,GHNA<c 


ij=  O.XPDLY<4>.XPaLY(:5; 
JU.iJjYPDLYC4>.YP0LY(:5: 
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PLT 

17 

PLT 

18 

PLT 

IS 
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£3 
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££ 
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PLT 
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PLT 

£7 

PLT 

£3 

PLT 

c!H 

PLT 

30 

PLT 

31 

PLT 

3£ 

F'LT 

33 

F'LT 

34 

PLT 

35 

PLT 

3-6 

PLT 

37 

PLT 

■jp 

PLT 

39 

PLT 

40 

PLT 

41 

PLT 

4£ 

PLT 

43 

PLT 

44 

PL.T 

45 

PLT 

46 

PIT 

47 

PLT 

48 

PLT 

49 

PLT 

50 

PLT 

51 

PLT 

5£ 

PLT 

53 

PLT 

54 

PLT 

55 

PLT 

56 

PLT 

57 

PLT 

58 

PLT 

59 

PLT 

60 
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CALL  SYMBOL  ■;  i  i ,  0  ?  VCEH  ?  0 .  : .  1 1 «  0 ,  0 ?  - 1  > 

CALL  SYMBOL  OiC?  YC?  ij.  t- I'^HFIHIi  E  ELEMENT  MESH.  0.  0,  19;:. 

CALL  PLOT  ';0=  0.  0.  U<99y> 

RETURN 

55   FORMhT    <3Fin. CO 

END 
DVERLAY(RDY?35  U> 

nVERLnYc:RnY!35  0> 
p-pnj~pAM  rij_  or;f>' 

COMMON  NNDHES?  NELEMHTs  HDiir?  HBAriO?  ND?  NT3s  IS  fGP?  NCYCLE?  LAYERS?  I 
INSTEp!  HT1£.  ETHj  HTl  ?  NT£..  NNTENSNi^  IFLAG-  HSIZE-  HCDriE';55ri:> ,  XCSSO;. , 

e>  <  JNDXOi  >  5  hnlsis?  ik<6k  e'>o;'  <  ARLf^Acesco  ?  ir-ffii-icssivj  ?  iNDEKceso:'  ■  g 

35>!ZAI(3) 
COMMON  --- 1  .-■  E  =  NUE  ?  RAD  I  US ,  yCEfi  ?  Yi. T N  ?  E I .  KN  -  KS  -HI?  Ha  -  I HTER 
COMMON  -'3''   R 

COMMON  ■•■■4.'-  MniiU';:26>?DELTf; 
COMMON  .•■'5.--  Q  <  1 1  0  0  >  f  L I  ST  <  !  1  0 1  > 

DIMENSION  RaiOiJ).  TEMPe/)-  LM<6::'j  SKClS?ie> 

DIMENSION  KP(.7y,    EPC 1 1  -  7> ,  SEP':!!,?^,  PSNRai.rJj  SPSNRai,?;) 
DIMENSION  EaO>!  NUE<10)f  PRDP0319> 
DIMENSION  Aa03jl03?j  ARRAY  a  03) 
DIMENSION  SR(6!6:;'j  SKJTCS.S) 
DIMENSION  TENPa7> 
REAL  NUEKjHuLY 
REAL  KHjKS 
REAL  NUE 
PI=e£..-V. 
ANISD=1. 0 
NDINC=i 

IF  <NT3.LE.  i)>  NDIHC=0 
IF  <:NT3,LE,  CD  NB=ND-1 
NEND=3*NT3 
NIi£=£*ND 

NLOF=£*NNDDE  :;.-' NT3 
CALL  HLILLMAT  ■:  A  :■  NS I ZE  •>  NS 1  2E  > 
CALL  NULLMAT  <SK? 1£- 1£> 
ChLL  zero  '::G'.'ND0F;' 

FJfiD  NUj'iihL  loads  FCK  ELEMlNTS  IN  THIS  PARTICULAR  LAYER 

DO  40  n  =  i!NLLEMNT 
NT=IX(?. n  J 

GO  TO  C40)4Cii.5.5.5>.  NT 
5    Il=IX<l! II ' 
I£=IX(£?II > 

i3=i:«;';3i.ii> 
Yc=';:rai>-i-Y<i£:)+Y(i3:):'  '-s. 

IF  (YC.GT.HE.OR.YC.LE.Hi;-  GD  TO  40 
nT=IK'::3j  II  ' 

wt=gamh<:mt>*areaa<  i  i  ;■ 

MK=NT-S 

GO   TO    CIOwJOjSS:"    MK 
.U  DO    15    1=1-6 

,5  LM(I>=£*IX':I?  IH+NTS 

GD   TO   35 
?ci  LMa>=3*l;v  1-  II)-1 

LM<£>=3*IX'  £j  I  D-l 

LM'::3>=£*IX' 3-  ID+NTS 


PLT 
PLT 
PLT 
PLT 
PLT 
PLT 
PLT 
PLT 
BLK 
BLK 
ELK 
STEP? BLK 
Y(550ELK 
AMA'::EBLK 
BLK 
BLK 
BLK 
ELK 
ELK 
BLK 
ELK 
BLK 
BLK 
ELK 
BLK 
BLK 
l:LK 
ELK 
ELK 
BLK 
BLK 
BLK 

blk: 

BLK 
BLK 
ELK 

:e;lk 

BLK 
ELK 
BLK 
ELK 
BLK 
ELK 
BLK 
ELK 
I:LK 
BLK 
BLK 
ELK 
BLK 
BLK 
BLK 

BLK 
BLK 
BLK 
BLK 
ELK 
ELK 
ELK 


10 
11 
1£ 
13 
14 
15 
16 
17 
18 
19 


40 
41 
4£ 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 


44 


LMi::4>=3*IX'.:4i.  ID-l 

LMi::5:5=£*i:v  'ii?  iI)+HT3 

L\-u.e>=a*iH'-.ef  ii>+nt3 

GD  TD,35 

£5 

DO  30  l^i,h 

3  0 

LMa>=£*IX(Ii.  ID+HTS 
LM':;3>=3*lK(3j  IH-i 

35 

Il=LMa> 

I£=LMi::£) 

I3=Lr'K3> 

I4=LM(:4> 

I5=LM(5> 

I6=Lrii::6> 

Q'::ii;;i=Qai>-WT--l£  = 

QO£>=Q';ie>-WT'--l£. 

Qa3>=Q(:i3>-WT.^'i2. 

Qi;:i4>=Qa4:;i-iJT.'-4. 

Qa5>=0<I5)-WT.'-4. 

Qi;ifa>=Qa6::i-UT.-'4. 

40 

CDNTINUE 

DO 

45  i=ii.ririar 

45 

QC 

n=Q(I>.'FLaiHTi:MSTEP::i 

REWIMIi   9 

FDRM   S1  IFFHEri    MATRIX    IM   BLIIICKS 

IMAT=0 

K3HirT=U 

HLit1ELH::=l 

NM=Nli 

HL=1 

THETA=IiELTAWpI.-'18  0. 

3LIP=TANaHETA;:' 

SLIP=1.1*SLIP 

MRITE    ■■6?310>    KSHIFT?NMMT.L.K,Hf1.nLiHD 

DO  £90  k;i  =  l.HELEMHT 

MT=IXi;7?KT> 

MT=IX'::8!KT;. 

EriDRnL=F:N 

ES=KS 

IF    '::NT.LE.£>    go    TO    85 

IF    aSTDP.LE.O.:'    GO    TO    35 

IF    (MT.NE. IMAT)    GO   TO   55 

GO   TO   85 

IMAT=rlT 

CALL   READrif;    i:;45  PROP?  319?  MT> 

ANISa=PKUP<t> 

FRRTIU=pR'npi::3::i 

r1u=3 

no  60  1=1 5  / 

MD=MO+i 

xpa;'=pRi:ipi:Ma) 

DD   65    1  =  1  J  11 
no  65   J=1..7 

M0=MD+1 
EPCIj  J)=PRaP<MD;:i 
DD  70  1  =  1'.  11 
DD  70  J=1j7 

Mn=MD+l 


ELK. 
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ELK 
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BLK 
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BLK 
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ELK 

94 

ELK 
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ELK 
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BLK 
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ELK 
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ELK 
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elk; 
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ELK 

10£ 

I:LK 

103 

ELK 
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ELf( 
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BLK 
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BLK 
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BLK 
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BLK 

1  09 

ELK 

lie 

BLK 
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11£ 
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70        SEPa?  j>=pRnp<:MD> 

DD  75  I  =  i.:i  1 
DD  75  J=l-7 
MD=MD+i 

75   PSNRa?j:i=pRap<:Mn:) 

DD  80  I=i5 11 

DD  80  J=iH7 
MD=M0+1 
80    SPSNR(IjJ>=PRnp';Mn> 

NP=INTCPRDP':319>) 
85    CDHTIMUE 

SEhRCH  FaR  f  IXMEMTS  EELDHiS  TO  THIS  LAYER 

GD  TQ  <3  0  ■  1  £  0 !  £  0  0  n  £  0  0  <  £  0  0  >  ?  f  ^T 

T'.'PE  I  ELEMEHTS 

90    DD  95  1=1 ?£ 

95   LMa:)=3*i;:::a«Kn-3 

DD  100  1  =  1 5  £ 

IF  '::Lna>  +  l.LT.NL>  GD  TO  100 
IF  <LM<n  +  l.LE.rHM)  GD  TD  105 
100    CDNTIriUE 

GD  TD  £90 
105     Ii  =  IK'::i?Kf> 
l£=TH<:£.KT> 


ChLL  R  I  r  iG  <.  K I ,  RHD I  US .  E I  <  VCEM  .  K  <  11  >  j  X  <  I  £  >  ? 

Y  <:   1 1 

JDIH  TYPF  I  ELEMEhTS 

DD  115  1=1 ?£ 

DD  115  K=i?3 

~ 

II=LM<i::'+K-KSHIFT 

IF  ai  =  LE,O.DR>II  =  Cr.MD>  GD  TD  115 

KK=3*I-3+K 

DD  110  J=li£ 

DD  ilu  L=1.3 

JJ=LMCJ::'+L-KSHIF  T-T  I  +  i 

IF  V  JJ.LE, 0)  GD  TD  110 

LL=3*.J-3+L 

A  <  J. .!  V  1 1  >  =A  <  J.J  <  I  I'J  +:::R  <  f-  r!  <  LL  > 

110 

COHTIHilE 

115 

CDfiTifiUF 
GD  TD  £90 

TYPE  II  ELEMEHTi 


i£0 

DD  i£5  1=1 

Lna;-: 

«IX<I»KI>-3 

1£5 

CuNTINUt 
DD  130  1=:: 

«4 

130 

Lria::'=£*i: 

DD  135  1=1 

'.:ii.Ki>-£+nT3 

.4 

IF  '::LM' 

t::'  +  l.LT.HL;) 

]D  TD 

1  3 

IF  .:lM' 

I>+1.LE.HM> 

jD  TQ 

;4 

135 

CDNTIriUE 
GD  TD  £90 

140 

ii=iMa-K; 

I£=I;^;£i.KT 

;' 

BLK 

113 

BLK 

114 

ELK 

115 

BLK 

116 

BLK 

117 

ELf;; 

118 

ELK 

119 

ELK 

1£0 

ELK 

1£1 

ELK 

!££ 

blk: 

123 

ELK 

124 

ELK 

125 

ELK 

126 

ELK 

127 

ELK 

128 

ELK 

129 

BLK 

13  0 

ELK 

131 

Et-K 

132 

ELK 

133 

ELK 

134 

BLK 

135 

ELK 

136 

ELK 

137 

BLK 

138 

BLK 

139 

BLK 

140 

BLK 

141 

BLK 

142 

ELK 

143 

BLK 

144 

BLK 

145 

BLK 

146 

ELK 

147 

ELK 

148 

ELK 

149 

ELK 

150 

ELK 

151 

BLK 

152 

ELK 

153 

BLK 

154 

BLK 

155 

BLK 

156 

BLK 

157 

BLK 

158 

ELK 

159 

ELK 

160 

BLK 

161 

BLK 

162 

ELK 

163 

elk: 

164 

ELK 

165 

BLK 

166 

ELK 

167 

ELK 

163 

ELK 

169 

BLK 

170 

BLK 

171 

ELK 

172 

46 


EH=EfiDRML 

YC=i;i.5*(:V':  i:i>+v<:i2)> 

CALL   READIi:"    (  1  ?  TEMP '  1 7i  l<  i  > 

CALL   READil;;    ■:;  1  ?  TEhP <  1 ;  .  I< I  + 1  ) 

IF    (IMTER.,LE„n::.    QQ   TD    160 

IF  <VC.GI„He>  GO  TD  If^O 

I F  < YC . GT . H  t . AND . NCYCL  E . EQ . i  >  GD  TD  1 6  0 

SIGMA='::TEMP(0+TENPa:  :'.■£. 

TAU='::TEPlP'.:e::'+TENP<e::'  fyd. 

IF  '::ABSCi:iGMA>.LT.GAf'ifi(r1T:J>  GD  TD  160 

IF  ^SIGnA+GAMAsiMT^J  IPS?  145- 145 

Ru=TAU'-SIGMA 

ratid=aes(rd;! 

if  (ratio.  gt.slip;i  gd  td  150 

GD  TD  160 

ES=  i  0  0  u . 0 

WRITE  (:6.320>  KI 

GD  TD  165 

ES= 1 0  0  0 . 0 

EN=ES 

WRITE  (b?  3i5;'  KI 

GD  TD  165 

ES= 1 0  0 . 0 

EN=ES 

GD  TD  165 

EH=ENQR('iL 

CALL  JQIi  n  ifKI  <  VrEhf  Et:.  ES<  X':  1 1  >  «  X'!  I£>  n  Y<  1 1  >  -.yf  Igi 

lK'::5,Ki::.  =  iriT<THl> 

I;<;<6jKI>  =  IHT(:TH£> 

TEnP'::3>=EN 

TEMP<4:>=ES 

CALL  WRITM5  a < TEMP. 17 ? KI > 

JDIH  TYPE  11  ELENEHT 


DO  175  K=i  '3 

n=LM':  I  J+K-KSHIFT 

IF  a  I  LE .  U .  DR .  1 1 .  GT .  HD >  '.'.Q   TD  175 

DD  170  J=li.4 
DD  170  L=lse 

JJ=Li1(J)+L-KSHIFT-II  +  l 
LL=£^J-2+L 

IF  CJJ.LE. 0'  GD  TD  170 
Ac:J.J.  II>=A(  JJ-  II,:'+SKJT-  KKjLL> 
li'O       CDNTINuh 
175     CDNTIHUE 

DD  1 3  U  I  =  1  ■  £' 
180    LM(I>=3*i;v  I-KI+l>-3 

DD  185  1=3-4 
1 S5    Lt'K  I  ■:•  =£*  I H  k  I ,  K I  + 1  >  -£+N  f  3 

Il=IX<l!KI+i> 

I£=IX(:£?KI-i-i  J 

CALL    JDIHT    (KI  +  l.  ,'CEH-  TH?  ESv  X  1 1  >  ?  Xa£>  ,  Ya  1  >  -  YOc 

1       t:.' 

IX<5.KI  +  l>-=INTaHl.J 
IX':;6i.k:i  +  i:;-IHT<TH£> 
TENP':3:'=^EH 


BLK 
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ELK 

174 

ELK 
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BLK 
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BLK 
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BLK 
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ELK 
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BLK 
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BLK 
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ELK 
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ELK 

184 

ELK 
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ELK 

186 

BLK 
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BLK 
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ELK 
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BLK 
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ELK 
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BLK 
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BLK 
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ELK 
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ELK 
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ELK 
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ELK 
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ELK 
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THl-THS.SKJT'BLK 
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ELK 
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BLK 

£01 

BLK 

£0£ 

BLK 
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BLK 
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ELK 
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ELK 
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ELK 
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ELK 
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BLK 
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BLK 

£10 

BLK 

£11 

ELK 
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ELK 

£13 

ELK 

£14 

BLK 

£15 

ELK 

£16 

ELK 

£17 

BLK 

£18 

BLK 

£19 

ELK 

220 

ELK 

£21 

BLK 

£££ 

ELK 
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BLK 

££4 

BL.K 
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ELK 

££6 

ELK 

££7 

::'.TH1.TH£.SKJBLK 

££8 

ELK 

££9 

ELK 

£30 

ELK 

£31 

ELK 

£32 
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TENP<4>=ES 

CALL   MPITMS    Q - TENP. ITjKI+I ) 

DO   1^5   1=1-4 

DO    195  K=l!c 

II=LMa>-t-K-KSHIFT 

IF    ai.LE.O.DR.II.GT.NriJ    GO   TO    195 

KK=£*I-e+K 

DO   iyO  J=l?4 

DO   190  L=l>£ 

JJ=Lr1'::  J  )+L-KSHIF  T-I  I+l 

IF    '  Jj.LE.O)    GD    TD    190 
ACjJ.  II>=A<.JJ.II>+SKJT'KK?LL> 
CDhT I nUE 
CDriTIMUE 


TRIAhGULAR  lLEMEHT  TYPE  III 

200    MK=HT-£ 

GO   TD    <£05=c!15!££0>?    hK 
£05  DO  £10   1=1 ? 6 

£10  LMa>=£*Irt<l5KI>-£+NT5 

GD   TD   £30 
£13  LM'::i>=3SIX<:i.Kn-3 

Ln<£>=3*IX<£jKn-3 

Lri<3>=£«I«<3?  KI  >-£+HT3 

LM<:4>=3*IX'::45KI>-3 

LM<;5>=£*IX<55KI>-2+hT3 

LrnS>=£*IX<e.5KI>-£+!HT3 

GD   TO  £30 
££0  DO   ££5    1  =  1:- 6 

££5  LrKI>=£*IXt:ijia>-£+MT3 

Lr-K3:>=3*iK':3jKn-3 

£30  DO   £35    1=1-6 

IF  'CLh':  i:>  +  i.LT.HL>  30  TD  535 
IF  'Xr-U,  I  ::•  +  !=  LE.rirO  -30  TO  £40 
£35    CONTINUE 

GD  TD  £90 
£40     I1=IX<1?KI> 
I£=IX<£jKI> 
I3=IX';35KI> 
'fC=iY<Il>+V':T£>+V'::i3>>  -3= 


VPE  lil  ELEMENT 


IF  (ISTuP-LE.  0>  GO  TD  I'rO 

IF  ';YC.GT,H£:j  go  TD  £6  5 

IF  <YC. GT, HI . AND. NCYCLE. EO, 1  •  L 

CALL  REAIifiS  <  1.  TEMPj  17.KI> 

STRESSX= ABS  <  TEMP  <?'>'} 

TAUaCT=AB3  C  TEMP  <  S  >  > 

TAUF=STRESSK*FRhT I 0 

RrT  I  D=THUaLT.'-T  AUF 

IF  <RATID=GT=  1  .  0>  WRITE  (6~3P5'> 

DO  £45  LI=i! 10 

R1=FLDAT'::LI-1>.'10  = 

r£=float':li>.-io. 


blk 

£33 

ELK 

£34 

BLK 

£35 

BLK 

£36 

BLK 

£37 

BLK 

233 

BLK 

239 

ELK 

240 

BLK 

£41 

ELK 

£42 

BLK 

243 

BLK 

244 

ELK 

245 

BLK 

246 

ELK 

247 

ELK 

£48 

ELK 

£49 

ELK 

250 

ELK 

251 

BLK 

252 

BLK 

£53 

BLK 

£54 

BLK 

£55 

ELK 

£56 

BLK 

£57 

ELK 

£58 

ELK 

259 

BLK 

260 

BLK 

£61 

ELK 

£6£ 

ELK 

263 

BLK 

£64 

ELK 

£65 

BLK 

£66 

BLK 

267 

BLK 

268 

ELK 

269 

ELK 

270 

BLK 

271 

BLK 

272 

ELK 

£73 

ELK 

274 

BLK 

275 

BLK 

276 

BLK 

£77 

ELK 

278 

BLK 

£79 

ELK 

230 

ELK 

281 

ELK 

2S2 

BLK 

ooo 

ELK 

234 

ELK 

£85 

BLK 

286 

BLK 

£87 

BLK 

2SS 

BLK 

289 

ELK 

290 

BLK 

291 

ELK 

292 
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LJ=11-LI-M  ELK  £93 

IF  (RKT[n.bT.Rl.ANr:  =  RATIIl-LE.R£>  bD  TD  £50  BLK  £94 

CDHTIHUE  ELK  £95 

LJ=S  ELK  £96 

ELK  £97 

LJ   *    IS    THL    ET   VS.    SIGhH3   CURVE    HD.    SFLECTEn   EASEC   ON   STRESS  RATELK  £98 

iiDW   SELECT    INTERVAL   OF   C0PFriIHQ7RRESSURE  ELK  £99 

IF    (STRESSX.  LT.  XP<  1 ;' .  DR.  STRESSX.  GT.  ftp(  HP^  >    WRITE 

XjKI 

HU   £55   LK=£5HP 

LI=LK 

IF    fiSTRESSX.  GE.  XP<LK-1  :> .  AHD.  STRESSX.  LE.  XP<LK>  ) 
CDNTIHUl 
LT=HP 

YriPi=SEP<LJ5LI-i:' 
YDFH 1  =SF  ;  f  IR  ■;  LJ «  L I  - 1  > 
VriP£=SEP<LJiLi:) 
TriPH£=SRSHRa-._bLI> 

¥l=EP<L.J«LI-i;' 

VHi=PSNR':Lj!.LI-l> 

Y£=EP(L-_NLI> 

¥ti£=PSriRa  J?LI> 

PP=STRlSSX 

ET  1  =DRD  T  MET  <  YDP 1 ?  YDP£  ?  X  i  5  X£ ?  Y 1 5  V£  ?  PP ) 

PSHi=ORriIHET';YDPHi  J  YDFNS,  XI  ?  ;<Sj  YNi  ?  YN£i.  PP> 

VDPi=SEP':LJ-i.LI-i::' 

VDPH1=SPSNR<L J-1 ?  LI-1 > 

YBP£=SEPa-J-l5LI> 

YDPH£=SP3nR<Lj-i . LI > 

Yl=EP<LJ-l!LI-l> 

YHi=PSNR(LJ-l)LI-l> 

Y£=EP(LJ-li.LI> 

YN£=PSNR'(,  J-ijLI) 

ETS---DRDINE  TC  YDPl .  YDP2!  Xl  j  XE;  Yl  j  t£j  PP;j 

PSN£=DRDINET'::YDPNi  !  YDPHfi.  Xl !  XE<  YNI  ?  YN£h  PP> 

EY=ET£+';ETl-ET£::'*<RATia-Rl) 

HUEY=PSN£  +  ■;  PSN 1-PSNE  :>  *  <  RaT  I D--R1  > 

EX=ANISnsEV 

NUEX=NIJEY*EX.-EY 

GD  TD  £75 

EY=E'::MT::i-iiJi:n;i.  0 

EX=ANISD»EY 

NUEY=NLiE(MT> 

NUEX=N!JEY«lX.-'EY 

GD  TD  £75 

EY=Ea'1T> 

EX=ANISG*EY 

NUEY=NUE(nT:J 

NUEX=NUEY*EX^EY 

CALL  TRANGLE  (KI 5  NTs  EX'  EY=  NUEX?  NUEYj  SK!) 

no  £85  I  =  l.=  & 

DD  £85  K=l5£ 

II=LN(r:^+K-KSHIFT 

IF  ai.LE.i:i.DR.II,Gr.ND>  GD  TD  £85 

KK=I  +  ';K-1>*6 

DD  £80  J=1j6 

DO  £80  L=l?£ 


ELK 

300 

6)330)  STRESSELK 

301 

ELK 

30£ 

BLK 

303 

ELK 

304 

GO  TD  £6  0    ELK 

305 

BLK 

306 

ELK 

3  07 

ELK 

308 

ELK 

309 

BLK 

310 

ELK 

311 

ELK 

31  £ 

BLK 

313 

BLK 

314 

BLK 

315 

ELK 

316 

ELK 

317 

BLK 

318 

ELK 

319 

ELK 

3£0 

ELK 

3£1 

ELK 

000 

ELK 

000 

ELK 

3£4 

ELK 

3£Pi 

BLK 

3£6 

BLK 

3£7 

BLK 

3£8 

ELK 

3£9 

ELK 

33  0 

ELK 

331 

BLK 

33£ 

ELK 

333 

ELK 

334 

ELK 

335 

ELK 

336 

ELK 

337 

ELK 

338 

ELK 

339 

ELK 

34  0 

ELK 

341 

BLK 

34£ 

BLK 

343 

ELK 

344 

ELK 

345 

ELK 

346 

BLK 

347 

ELK 

348 

ELK 

349 

ELK 

350 

BLK 

351 

BLK 

■Pi^,^ 
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.JJ=iJ1':;j>+L-KSHIF  r+l-H 
IF  '.  JJ=LE,  O)  GO  TG  £S0 
LL=J+'::L-i>*t. 

H  <  J  J  ■  n  >  =A  ■;  J  J  -in  +sk  <  kk  .  ll  > 

PyO  ClIHTIHUE' 

£S5  CUHTINUE 

S90   CDHTINUE 

IJPITE   STIFFNESS   EQUftTUIGMS   DH   TAPE    1 1: 

ijRITEaiVi    (<:H<M<!1>?ri=i5MBAHl!>?M= 

CALL  SERD  'J-iRRAY;  HSIZE> 

NX=HL- 1 

DiJ  300  M=l5h;iZE 
Du  £95  H=i!NSIZE 
£95    ARRAY  (;N:J=Aai!N:> 

CALL  WKITMS  <i  OjRRRAYjiiSIZE!  f';«:+M> 
SOU  CDhTIHuL 

IF  CNM.GE.HEDF)  GO  TO  305 

CALL   HULLMAT    <A? HSIZE. MSI lEJ 

KSHIFT=KSHIFT+Nri 

IF    <HUMBLK.EQ=i>    riD=HIi-HEIHC 

HLir'iBLK=HUMBLf:  +  i 

HM=HM+ND 

nL=HM-nIi+i 

QD   ID   50 
305    IF    '::NT3.GT,0>    riD-riD+NDIHC 

RETURN 

310  FGRMAT  i5H^       9HK,SHIFT  =  -IS^SX?   9HrHUMELK  =  <  I'; 

iSH    NL  =  5  15-   7H   HD  =  ^  15  :■ 
315  FDRMAT  (SX?  £41-!  INTERACT  I  DH  ELEMENT  Ha=  ^15-'3>i^ 

i  fHS:  I  UN  > 
3£0  FDRMAT  CSK?  26HINTERACTIDN  ELEMENT  ND.  =  -ISsS; 

1  EXCESSIVE  SHEAR;' 
3£5  FuRMAT  <10H^    4iHSTRESS  RHfin  IS  GREATER  THAN  1 

i  ITHIN  ELEMEnT  ND.  =  j  15;' 

330  FDRNAT  (AQX,    56HCDNFINING  PRESSURE  IS  DUT  DF  RANGE  DF  SIGNAS.  STREELK  391 

iSSX  =  jEIS.S-SX!  17HIN  ELEMENT  HQ.  =  ?  15;'                         ELK  392 

ELK  393 

EHIi                                                             ELK  394 

SUERGUTIME  TRhNGLE  ';K?  NT?  EK?  EYs  riUE;^  NUEY;  SK>                      TRG  £ 

CDMMuN  NNDHEC  ,  HELEMNT  5  NDDi-"  5  ME  AND ,  ND  ?  NT3  5  I  STOP  -  HCYCLE  -  LAY'ERS  j  I  STEP  ?  TRG  3 

INSTEP?  NT  l£.ETA«NTl!NT£!.NDrENSN?  li^LAG?  NSIZE- NCnDE'iSSOj  ?  X<:55n>  ?  Y'^SSOTRG  4 

£;'?JNDy'::51>?ANLSIS)  IX '35?  £5  L?  ?  ARE' iA';£5  0>  ?  INDX';£50:'?  INDEX<£50>  j  GAMA'^eTRG  5 

35::'!ZAI'::3>                                                     TRG  6 

CDMMDN  ■•  1/  E'.  HUE?  RhBIUS?  X  ;EN?  YCLN^  EI  ?  KN?  KS?  HI  ?  H£?  INTER            TRG  7 

CQMMDN  .■-■£.--  n                                                     TRG  S 

DATA  PLSTRSsPLSTRN.--6HPLST^^S.6HPLSTRN.-'                           TRG  9 

REAL  NUEXjNUEVjHUENS                                             TRG  10 

DIMENSION  SK<l£s  !£>?  EaO'?  NUE'riO;'?  D<10i»KO                     TRG  11 

DIMENSION  DD'::3?3;'?  DMATC90;'                                      TRG  IE 

REAL  KH?KS                                                      TRG  13 

TRG  14 

****a!«*sw*?ihft**s«*******K?!>*w?*;i*;****«?i(w;*w  15 

GENEkhTE  STIFFNESS  MATRIX  FDR  TRIANGULhR  ELEMENTS        TRG  16 

TRG  17 

AA=AREAACK;'  TRG  IS 

MT=IX<85K>  TRG  19 


BLK 

353 

ELK 

354 

ELK 

355 

ELK 

356 

BLK 

357 

BLK 

358 

BLK 

359 

ELK 

360 

E  10 

BLK 

361 

ELK 

36£ 

?hd:) 

ELK 

363 

ELK 

364 

ELK 

365 

ELK 

366 

ELK 

36? 

ELK 

368 

ELK 

369 

ELK 

370 

BLK 

371 

BLK 

37£ 

BLK 

373 

ELK 

374 

ELK 

375 

ELK 

376 

ELK 

377 

ELK 

378 

ELK 

379 

BLK 

380 

ELK 

OCj| 

BLK 

382 

8H   NM  = 

1 5  5   BLK 

383 

ELK 

004 

IHHAS  FAILED 

IN  TEBLK 

385 

BLK 

386 

31H  HAS  Ffili 

.ED  INBLK 

387 

BLK 

388 

RATIO  =  ?F6 

P.5X.ELK 

389 

ELK 

390 
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=  1 .  'AA 


TYPE    III    ELEMENL 


L=IX 
M=Ir-=: 

H=i;«; 


-XCL 


;n) 


d  J=  r 

r' 


Y31  = 
IF  < 
GO  T 
Iiii  = 
ril£= 
Iii3= 
ri£:l  = 
Ii£S= 
Ii£3= 

ri3  i = 

ri3£= 
ri33= 

GD    T 


:r1;'-Y<N> 
:h::'-v<l) 

^TDP.LE.n;'    go 
10 

iriTs  !;■ 
:mt!.2> 

J'iT 

:riT 
:mt 
:riT 
:mt 
:mt 

15 


Al 


'■  HflL 

(AHL 

£=DE!': 
3=DD  ( 
1=DD':: 
£=DIi':: 
3=DIn: 
1=DDC 
£=DrK 
3=Dri':: 
=Iili* 
=ri3i* 
=riii* 
=ri3 1  * 

=D11* 

=ri3 1  '*■ 
=rii£*: 

=D3£* 

=11  le*; 
ij=ri3£^ 

1=D1£ 


SIS.EQ.PLSTRS>    CALL   PLNSTRS    CEX? EYj MUE>^  NUEYj DD) 

SIS.EQ.PLSTRH::'    call   values    < Eft jEYhHUEX?  nicy jDD> 

1>  1> 

i 

1 


23+D33*X3£ 
31+D13*:«;i3 
31+D33«K13 
ri£+Di3*<X£l 
T'i£+D33*;«;£i 
K3£+ril3*Y£3 
«:S£+ri33*T£3 
«;i3+l!l3*V3i 
*X13+ri33WY31 
■KJsEi+niS^YlE 
*k£1+D33*t1£ 


B 1  =D£  i  *  Y£3+DS3*;>v3£ 

Be=r!£  1  *  y3  i  +rie3«x  i  3 

B3=D£  1*Y  1  £+D£3*:X£  1 
B4=D££*>i3£-Hi£3*Y£3 
E5=D££*N  1 3+IiS3*Y3 1 
B6=D££*H£l+r!£3*Yl£ 


::ii  = 


Y£3*Ai  t-y;3e 
Y£3*A3t;»;3£ 


TRG 

2C! 

trg 

£1 

TRG 

££ 

TRQ 

23 

TRG 

£4 

TRG 

£5 

TRG 

£6 

TRG 

df 

TRG 

'-"P, 

TRG 

£9 

TRG 

30 

TRG 

31 

TRG 

32 

TRG 

P.P. 

TRG 

34 

TRG 

35 

TRG 

36 

TRG 

37 

TRG 

38 

TRG 

39 

TRG 

40 

TRG 

41 

TRG 

4£ 

TRG 

43 

TRG 

44 

TRG 

45 

TRG 

46 

TRG 

4? 

TRG 

48 

TRG 

49 

TRG 

50 

TRG 

51 

TRG 

5£ 

TRQ 

53 

TRG 

54 

TRG 

55 

TRQ 

56 

TRQ 

5? 

TRG 

58 

TRG 

59 

TRG 

60 

TRG 

61 

TRG 

6£ 

TRQ 

63 

TRG 

64 

TRG 

65 

TRG 

66 

TRG 

67 

TRG 

68 

TRG 

69 

TRG 

70 

TRG 

71 

TRG 

7£ 

TRG 

73 

TRG 

74 

TRG 

75 

TRG 

76 

TRG 

77 

TRG 

78 

TRG 

79 
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Ci4=V£3*A7+> 
C15=Y£3*h9+;' 
Cit.=Y£3*iiilH 


Sl=¥31*HirXi3*;H£ 

&£=Y3  i  *A3+K 1 3*h4 

£3=Y3 1 *h5+X 1 3SA6 

C£4=Y3  i*H7+X 13*A8 

5=Y3i*A9+Xi3*HlO 

6=Y3i*Aii+;a3*Ait 


C3i= 

i'i£*Hi+X£i*A£ 

i::3£= 

r'i£*A3+J':£l*A4 

C33= 

( i  £* A5+X£ 1 «A6 

C34= 

rl£*A7+X£i*A8 

C35= 

t'1£*A9+K£1*A10 

C36= 

T'i£*Rii+;:£i*Al 

C4i  = 

<3£«El+Y£3«Ae 

C4£= 

-:;3£*B£+Y£3*A4 

C43= 

«£SB3+Y£3*A6 

C44= 

<;3£*B4+Y£3*AS 

C;45= 

■::3£*E5+Y£3*h10 

C4fa= 

«£*B6+Y£3*Ai£ 

5  i  =F^  1 3*'B  1  +  Y3 1  •*«£ 
5£=Xi  3*E£-i-Y3  1  *A4 
53=K  i  3*B3+ Y3 1 *A6 
54=K 1 3*B4+ Y3  i  *A8 
C55=X13«B5+Y31«A10 
56=X  i  3*E6+*v'3 1*A  1  £ 

63=X£l*E3+Yi£*Ab 
bi=H£l*Bl+Yl£:*A£ 
6£=K£1 *B£+ Y 1 £* A4 
64=k£  i  *B4+Y 1 £*A3 
65=X£1*E5+Y1£«h10 
6t.=X£  i  *B6+ Y  i  £*  A 1  £ 


IF  (K.GT.NTiS.HNri.K.LE.  <MTi£+MTl.'£>>  GD  TD 

Qu  TD  £5 
£0  ETTA=lTA-'£. 

ETATMn=l.-£=*t;TTA 

ETAFnUP=l,-4,*ETrA 

ETA£=ETATWU*ETATUD 

ETA4=ETHFDUR*tTH!-XllJR 

GD  TO  30 
£5  ETTA=0. 0 

£TAFOUK=i. 0 

ETATWD=ETAFDUK 

ETA4=1.U 

ETA£=ETA4 


SK 

i5i>=CilKT 

SK 

£»i;'=-Cl£* 

SK 

l5£)=SK<£= 

SK 

■1j3:)=C13*T 

SK 

:3.i>=sKa:. 

:i  =  +4.*etia.--3, 

•Cl.E.^ETAFaUR; 


■4/*ETTA*ETTA:'-'-<ETATUD*ETATWC 


'^ETTAH3.*ETTA*ETTA>-'(4.*ETATWa) 


trg 

80 

TRG 

81 

TRa 

8£ 

TRG 

83 

TRG 

84 

TRG 

85 

TRG 

86 

TRG 

37 

TRG 

88 

TRG 

89 

TRG 

90 

TRG 

91 

TRG 

9£ 

TRG 

93 

TRG 

94 

TRG 

95 

TRG 

96 

TRG 

97 

TRG 

98 

TRG 

99 

TRG 

100 

TRG 

101 

TRG 

10£ 

TRG 

1  03 

TRG 

104 

TRG 

105 

TRG 

106 

TRG 

107 

TRG 

108 

TRG 

109 

TRG 

110 

TRG 

111 

TRG 

112 

TRG 

113 

TRG 

114 

TRi3 

115 

TRG 

116 

TRG 

117 

TRG 

118 

TRG 

119 

TRG 

1 8  0 

TRG 

1£1 

TRG 

1££ 

TRG 

1£3 

TRG 

184 

TRG 

1£5 

TRG 

1£6 

TRG 

1£7 

TRG 

1£8 

TRG 

1£9 

TRG 

130 

TRG 

131 

TRG 

13£ 

TRG 

133 

TRG 

134 

TRG 

135 

TRG 

136 

TRG 

137 

TRG 

138 

TRG 

139 
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SKf  £?3>=-C?3«TH.-'i;E:.  TRG  140 

SKc:3!.e)=SK<£i.3;'  TRG  141 

SK < £ !  £  >  =C££STh. ••  (. 4 .  *ETA4 ;■  TRG  1 4£ 

SK  i  3  J  3  >  =C33*T  A*  ■;  i .  +3 .  *ETT  A.--3 »  +  i  6 .  *ETTA*-:ETTA  )  -A .  TRG  1 43 

TRG  144 

SK<i.4>=TA*<Cii«£.*ETTA-H:i£«<4=  ■3.+8=*ETrA.'3.::'>--<4.*ETft£«ETATMn>   TRG  145 

SK'::4!l>=SK<l-4>  TRG  146 

SK<£!4>=C£i*TH.-'::3.*ETA£*LTAFaLiR)  TRG  14? 

SK';4j£>=SK'::£?4>  TRG  148 

SK<33  4>=4=s(C3i+C3£>STA«ETTH.'<3.*ETA£>  TRG  149 

SK<:4?3>=SKC3;4>  TRG  150 

SK(:43  4:;  =  '::£.  *C1  l+Cl£+CSl+£  ,  ^CSe^s^TA--  <3=  *ETA£*ETfte>  TRG  151 

TRG  15E 

-■C3.*ETHTM.a«ETAFniJR:)  TRG  153 

TRG  154 

TRG  155 

TRG  156 

5*4,«ETTA>WTA.'-(3.*ETAFnUR>  TRG  157 

TRG  153 

E:3-'^'<3.*ETAFDUR;«ETh£)  TRG  159 

TRG  160 

.*r  ^3 >.■■'< 3.  *ETA4>  TRG  161 

TRG  16£ 

^Cl3::0^•'(:3,*FTATWD*ETA^DUR>  TRG  163 

TRG  164 

TRG  165 

TRG  166 

i-C33>>.--<3,*ETAFaiJR>  TRG  167 

TRG  163 

i3>'<3=*ETAFDUR*ETA£>  TRG  169 

TRG  170 

53 >.'■':; 3,  *ETh4>  TRG  171 

TRG  17£ 

.*C33:5.-<3.*ETA4:i  TRG  173 

TRG  174 
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QRriIHET=Tl+T£+T3+T4 

RETURN 

END 

SUERDUT 1 NE  K I I'iG  (R I  ?  R «  E I  ?  r'CEH •  X  I  j  H£ -  Y 1  s  Y£ ?  SK ) 

DIMENSIDH  SRf6!6>!  SK<6?6> 

DIMENSIDN  TLCSr-fc.:)!  TCCfaje) 

DIMENSIDN  ThP<36> 

REAL  NUE 


STIFFHESi  MATRi:-:  FUR  RING  ELEMENT; 


SP AN=SQRT  < < X 1  -XE  >  **£+  <  Y 1  - Y£  >  «*£ ) 

BETA=SPAN-'R 

BETAS=£.*BETA    : 

SN=SIH':.BETA:^ 

CS=CnS':;BETA) 

SN£=SIH( BETAS > 

A=BETA-SN 

B=CS+SN*SN^--e.  -i, 

C=l .  5*EETA-£,  «SN+SN£.--4. 

D=  ij .  5*EETA-SNS>-'4 . 

E=CS- 1 , 

AA=E*E-'EETA-ri 

BE=B-AWE.-'EETA 

CC=A«D-E*E 

Dri=H*R.-EETA-r: 

EE=C*E-ASB 

rr=E*E-C*D 

G=B*'::B-£.  ^A*t. •'BETA) +C*^ (:£*£■■'■?£ Th   ri>+H*f: 

ft=ei.--<<rw*3>*g::' 

S:R';is  1>=FT^Ah 
SR(1?£)==FT«BB 

SR(£f.i>=SR(:i,e::' 

SRa !  3>=FT*CC:*R''BETA 

SR';3!  i>=SR(i!3::' 

SR<£.-£>=FT*riri 

SR'::£!  3>=FT*EE:«R.-'-EETh 

SR'::3j£::'=SR<:S?3> 

SR<3j  3>=FT*FF  *R*R--BETA 

SR(4!  i  >=-FT*':;Aft*CS+BE*SN::' 

SRa!4>=SR(:4i.  1> 

SR<4!E>=-FT*':.  bE*CS+DD*SN  ^ 

SR';£?4>=SR-;4j£> 

SR  <  4 r  3  ;■  =-FT* ■::  CC*CS+EE*SN  >  WR.-'BETA 

SR':3!4>=SR(4!3> 

SR  >;  1 !  5  >  =FT*  •'.  AA*S!-i-EE*CS  > 

sRkSi.  i)=SRaj5;' 


drd 

o 

rjRD 

4 

□RD 

5 

□RD 

6 

ORB 

7 

DRD 

R 

DRD 

9 

DRD 

10 

ORD 

U 

DRC 

IE 

DRD 

13 

DRD 

14 

DRD 

15 

QRD 

16 

DRD 

17 

DRD 

18 

RNG 

£ 

RNG 

3 

RNG 

4 

RNG 

5 

RNG 

6 

RNG 

10 

RNG 

11 

RNG 

IE 

RNG 

13 

RNG 

14 

RfiG 

15 

RNG 

16 

RNG 

1? 

RriG 

18 

RNG 

19 

RHG 

SO 

RNG 

£1 

RNG 

EE 

RNG 

£3 

RNG 

£4 

RNG 

£5 

RNG 

E6 

RHG 

£7 

RNG 

£S 

RNG 

pC( 

RHG 

3  0 

RHG 

31 

r-,1  1'J 

RNG 

33 

Rr-JG 

34 

RHG 

"Dcr: 

RNG 

36 

RNG 

37 

RNG 

oo 

RHG 

39 

RNG 

40 

RHG 

41 

RNG 

4E 

RNG 

43 

PN|- 

J.  A 

57 


SRC55£)=FTS(BB*SH-DriSCS::' 

SR<:£55>=SR(55e:' 

SR  ( 5  J  3  >  =rT*  ■;  CCSSH-EE*CS )  ^R.-'BETA 

SRi;3j5:'=SRi;5.«3> 

SR<6«  i:J=FT*(:ftii|*';:CS-l.  >+EI;*SN-Cr.'BETA>*R 

SR<ljb>=SR'::6j  i> 

SR'::6?S>=FT*(:BE*';CS-1.  :;'+Dl*Sfi-EE-BETA>*R 

SR<£?6:)=SR':'6?2;' 

SR(6!.3>=FT*'::CC*(CS-1=  J'+Et  *SN-FF  ;'*R*R---EETA 

SR<3!6>=SR'::6j3> 

SR<4j4::'=SR<l5l> 

SR<4!5::'=-SRaj£> 

SRC5!4>=SR'::4.5::' 

SR<4j6>=SR';iii3;' 

SR  <  6  ?  4  >  =SR  i!  4  J  6  > 

SR';555>=SRi::e'£;' 

SR<5j6:)=-SRi::Sj3;i 

SR<6j5:)=SR<5j6> 


iK':.fc.jb.: 


TRAHSFaPI-lfiTinN  GF  CQQRDTMATE  SYSTEM 


CALL 

NULLMAT    CTL 

!  b  !  6 

,;) 

CALL 

NULLMAT    <TC 

■  636 

,> 

T1  =  T 

1-YCEH 

t£=y; 

E:-YCEN 

ih    0 

■•%  i .  Eu  ■  u .  0  ;■  X 

1=0. 

U  0  0  0 1 

IF    <; 

Kd.EQ.  u .  Ij >    Xi 

£-0. 

Ij  u  Ij  Ij  i 

TNi=Ti.-Ki 

TH£=" 

r£.--H£ 

PHI1  = 

=rtaH';tni> 

PHI£= 

^ATAN'CTNS:' 

CS1=( 

;as(PKii> 

SH1=; 

iIH(PHIl> 

CS£=i; 

JuSCPHIc') 

SH£=; 

i:IN';PHl£":' 

TLU; 

.  0=-SHl 

TL;!: 

.£>=+CSi 

TL(£. 

. 1>=-CS1 

TLcS. 

.£;)=-SHi 

TLC4; 

.4:j=-SN£ 

TL<4; 

i5>=+CS£ 

TL<5; 

.4>=-CS£ 

TL<5; 

.5;'=-SN£ 

TL<6: 

1 6  > = i .  0 

JU.3: 

.3;'=TL(6!6> 

CALL 

MhTMULT    <SK: 

-  TL- 

TC1.6'' 

CALL 

MATRMUL    (TL; 

=  TC? 

SKi.r.|.i 

riQ=u 

DD   5 

1  =  15  6 

DO  5 

J=l5  6 

Mn=Mn+i 

TMP'.Mu>=Sk<Ij., 

J  ':> 

cdht: 

[NuE 

CALL 

UPITrtS    <.?5Tf 

".P?!-! 

D^Kn 

RETURH 

END 

SUERDUTIHE   jniHT 

(KI 

J  VCEf-i 

DIMENSION    SKjT03: 

•  8;') 

TL<t 

RHG 

46 

RNG 

47 

RHG 

4S 

RNG 

49 

RHG 

50 

RNG 

51 

RHG 

5E 

RNG 

53 

f:NG 

54 

RHG 

55 

RHG 

56 

RHG 

57 

RNG 

58 

RNG 

59 

RNG 

60 

RHG 

61 

RtHG 

62 

RHG 

63 

RHG 

64 

RriG 

65 

RHG 

6-6 

RtiG 

67 

RHG 

68 

RHG 

63 

RHG 

70 

RNG 

71 

RNG 

7£ 

RNG 

73 

RNG 

74 

RNG 

75 

RNG 

76 

RNG 

77 

RHG 

70 

RtHG 

7''^ 

RHG 

8  0 

RHG 

81 

RHG 

S£ 

RHG 

S3 

RriG 

34 

RHG 

85 

RriG 

86 

RHG 

87 

RHG 

38 

RHG 

83 

RNG 

90 

RNG   91 
RNG    9E 


. VIhYE'THIsTHSkSKJT) 
.  TMP*:64;:' 


RNG 

94 

RNG 

95 

RNG 

96 

RNG 

97 

RNG 

'58 

RNG 

99 

RNG 

1  0  0 

RrUG 

101 

RHG 

1  0£ 

RHG 

1  03 

JNT 

£ 

JNT 

3 

58 


EQU I  VALENCE  <  TC  ■;:  1  -  1 ':'  -  TMP '  1  >  > 

STIFFNESS  MATRIX  i-DR  IMTERACTinN  ELEMEN 

PI=£i2..--7. 

SPAN=SQRT';  (:;■■■; i-H£ >**£+(  Vl-  V£ >**£;' 

Ti=Yl-YCEH 

T£=Y£-'rCEN 

Tl=-Ti 

T£=-Te 

IF  '.! T 1 .  EQ .  0 .  0 >  T 1  =  0 ,  0 0 0 0 ^ 

IF  afi.ER.O.O)  Tl^O.OOOO; 

TNi=ABS<:>a.--Ti:) 

TN£=HBS'::H£'-T£> 

THETRl=ATAfKTNi;' 

THETh£= AT AH  <  TH£  > 

IF  (  T i  .  LT .  U .  0 >  THETA 1  =££ .  .'■? .  -THETA 1 

IF  C  I  £ .  LT  =  iJ  =  0 ;■  THETA£=££ ,  .•-? .  -THETAS 

THi  =  i  0.  *u  HETHl*i80.  .■■■PI  ::■ 

THe= 1  0 .  * '■: THET A£-*  1 8 0 .  .-p I  .:■ 

CALL   NULLf'iAT    <  SK  JT s  8 ?  S ::■ 

CS=SPi4N*ES.-6 . 

CN=SPAfl*EN.-fc., 


SKJT 

<1 

?  1 

>=£. *CS 

SKJT 

(.1 

>3 

;'=CS 

SKJT 

a 

J  5 

>---CS  ^ 

SKJT 

'::i 

)=-£=«CS 

SKJT 

':!£ 

,  £ 

>=£.*CH 

SKJT 

■::£: 

=  4 

>=CN 

SKJT 

';£i 

.6 

>=-CN 

SKJT 

'•'  £  ■ 

i8 

)=-£.*CH 

SKJT' 

t:3: 

.  i: 

>=CS 

SKJT' 

'  3'. 

:'=£=*CS 

SKJT' 

( :3 1 

1  5! 

J=-£.*CS 

SKJT' 

C  3  ? 

,  7; 

:'=-CS 

SKJT' 

;4! 

'£: 

:'=CN 

SKJT' 

■Mi 

i4; 

:'=£,*CN 

SKJT' 

::4, 

i  6; 

>=-e„*CN 

SKJT' 

'As 

8 '. 

:'=-CN 

SKJT' 

■.5n 

i: 

:'=-C:S 

SKJT- 

-.5, 

3; 

:'=-£.*CS 

SKJT': 

:5, 

5; 

»=£,*CS 

SKJT- 

:5, 

7: 

'=CS 

SKJT< 

'.6} 

£; 

"=-CN 

SKJT': 

'€■} 

4; 

■=--£.  *CN 

SKJT': 

'.  6 ) 

6!: 

'=£.»CH 

SKJT< 

!6! 

8; 

'=CH 

SKJT': 

!7? 

1:: 

'=-£.*CS 

SKJT': 

:7? 

3; 

'=-CS 

SKJTt 

:?) 

5) 

i=CS 

SKJT': 

!?? 

?) 

'=£.*CS 

SKJT': 

!&! 

£) 

'--£.«CH 

SKJT': 

;8j 

4) 

'=-Cn 

SKJT< 

!Sj 

6; 

■=CN 

SKJT  ■; 

i-ii 

8) 

=£.*CH 

CALL 

HULLMAT  <TL5  8<8> 

CALL 

NULL 

.MAT  ■::TC5  8j8) 

CSi=r 

:;ds 

:'::  THETA  1;' 

SNi=S 

■■!!■' 

l<THETAi> 

CS£=C 

:D3 

:!:THETA£:' 

JNT 

4 

JNT 

5 

JNT 

6 

JNT 

7 

JNT 

S 

JNT 

9 

JNT 

10 

JNT 

11 

JNT 

IE 

JNT 

1:3 

JNT 

14 

JNT 

15 

JNT 

16 

JNT 

17 

JNT 

18 

JNT 

19 

JNT 

£0 

JNT 

£1 

JNT 

££ 

JNT 

£3 

-..INT 

£4 

■JNT 

£5 

JNT 

£6 

JNT 

£7 

JNT 

£8 

■JNT 

£9 

■JNT 

:jO 

JNT 

31 

-JNT 

3£ 

JNT 

33 

JNT 

34 

JNT 

35 

JNT 

3^3 

JNT 

37 

JNT 

•3y 

■JNT 

39 

JNT 

4Q 

.JNT 

41 

JNT 

4£ 

JNT 

43 

.JNT 

44 

..JNT 

45 

JNT 

46 

JNT 

47 

JNT 

48 

JNT 

49 

JNT 

.5  0 

JriT 

51 

JNT 

S£ 

.JtiT 

5:3 

■JNT 

54 

JNT 

55 

.JNT 

56 

JNT 

57 

JNT 

58 

JNT 

59 

JNT 

60 

JNT 

61 

.JNT 

6£ 

JNT 

6':3 

59 


:iM£=SIH(THETAP> 
TL'::7j?>=C:Si 
TLi:!"  1>=TL(;7j7;' 

TL';i!t£>=TL'::7!iy;' 
TL(:3.7>=-SH1 
TL'iE)  i>=TL';8:i7;' 

TL'::£!.e:j==TL(;8'.8;' 

TL':5!5>=CS£ 

TLC3!3>=TL(5'i5> 

TL(;5j6::'=SH£ 

JL(.3,4':>=lL':5^bj 

Tl<:6?5>  =  -SH£ 

TL<4j3::'=TL'::iE.i.5> 

TL(:6-b>=CS£ 

TL'::4«4>=TL'::t.»6;' 

CALL  MATnULT  (SKJTj TLj TC? yjS;' 

CALL  HATRPiUL  'CTL!  rC?  SKJT5  Ss  8> 

DD  5  I  =--  i  J  8 
DD  5  J=l!8 

nn=MD+i 

TMP<:ND  ;•=<»:  JT'::  I  ?J> 
5  CDNTIHUE 

CALL  URITMS  i.?)  TMP.  MDs  KI ) 
RETURN 

lND 

aVEPLftY<RaY<4jO> 

□yLRLAV<RaYii4.0;' 

PRDGRAM  STRESS 

CDMMDH   NNaiiES?  HELEriHTf  HiiDF!  MEAHHf  MD;  NTS-  ISTDp.  HCYCLE?  LAYERS?  ISTEP 
iHSTEP<HTi£!ETA!NTi5NT2sNDrEHSM!  IFLAG- hSIZE?  HCDDE^SSCDs  X<:55u>  j  Y(55Cr 
£;■  ?  JNDXcSl  >  !  AHLSISf  IX(8?  £50>  ?  AREf  lACESO  -  lNDX'::£5u)  ?  mW'i'iS^O':'  >  GAMA^S. 
35>5ZAI(3::' 

CDMMDM  -'l-'  E.HUE<RArUUS!:«'CEH»Yi;LN!EI.KH!KS|,Hl.H.='.  IMTER 

CuMfiDft    ■•■-£■■     Ii 

CDMMDti   -'3^   R 

CahlMuH   ••■'4.'-    hULiU<£6>?  DELTA 


REAL  NUEl5NLiL£-HLiE 
niMENSIDH  TEMP a 7. J 
III  MENS  I DN  DMhTOO) 
I:  I  MENS  I  ON    In'50?4> 


LISTdlUl  :' 
NUEKjKSjKN 

Raico;'!  E'lo: 

lja£>.  SKJTClt 
FR<£> 


HUE'.  10..' 5  !!■  10?  ]  O;:' 
l£>j  PPa£>!  TENP(i; 


THIS  RDUTIhE  COMPUTES  STRESS-  STRAIN.  QCTAHEADRAL  STRESSES  AND 
STRAINS  FRDM  COMPUTED  NDDAL  DI  "PLACEMENTS  STORED  AS  Ra>  DN  TAPE 


?;S:ft**S'i 


K**;s«ft«*at5^*«s*i**is*-ii«***«sw«s***i*?K*i*s;**K**i« 


CALL  ZERD  (Q^NDDF::' 

IFLAG=0 

THETA=IiEL  rA*P I.-' 1 8 0 . 

SLIP=TAHaHETA:' 

SLIP=1.1*SLIP 

DD  100  KEL=1 iMELEnNT 


JMT 

64 

jrrr 

65 

JNT 

66 

JNT 

67 

JNT 

68 

JNT 

6'? 

JNT 

70 

JNT 

7! 

JNT 

7£ 

JNT 

73 

JNT 

74 

JNT 

75 

JNT 

76 

JNT 

77 

JNT 

78 

JNT 

79 

JNT 

8  0 

JNT 

81 

JNT 

Re 

JNT 

S3 

JNT 

84 

JNT 

85 

JNT 

86 

JNT' 

87 

JNT 

P.p. 

JNT 

k9 

JNT 

90 

JNT 

91 

JNT 

9£ 

STR 

1 

STR 

£ 

STR 

3 

STR 

4 

STR 

5 

STR 

6 

STR 

7 

STR 

8 

STR 

9 

STR 

iO 

STR 

11 

STR 

1£ 

STR 

13 

STR 

14 

STR 

15 

STR 

16 

STR 

17 

STR 

J.  0 

STR 

19 

STR 

£0 

STR 

£1 

STR 

££ 

STR 

£3 

STR 

£4 

STR 

£5 

STR 

'-'f. 

STR 

£7 

STR 

£8 

STR 

£9 

STR 

30 

STR 

31 

60 


iTVP=i;-;';7:.KEL> 

MT=lK'::8=KEL> 

GD   TO    cIOOsSjSS.SS-SS/i     ITYP 


TYPE  II  INTCRACTIDN  ELEMENT 


CALL  ZERO  (U:!  IS;' 

DD  iO  I=ij2 

Du  1  u  -1=  1 J  £ 
K=a-i:>«£+J 
M=3«I;■^::I.KEL>-3+J 
Ir(K:EL;I>=N 

U';K:'=Rao 

CDHTINUE 
DD  15  1=3-4 
DD  15  -J=  1  ?  2 
K=  a- !;■*£+ J 
H=£*lK';  I?  KEL>-£+HT3+J 
IA<KEL5l>=H 

U'::K>=R<H> 
CDNTIMUE 
I  =  IKa;KEL::' 
J=IX'::ei.|<EL> 

SPAH=SQRT<  (.Hi  1>-)«J':>  )^^*£+< 
THl=FLnAT<IX'::5i.KEL>> 
TH£=FLuFiT''  Il'-Xb!  KEL!'  > 
THETp  1  =TH  1  *:p  I  '■■■  1 8  ij  0  = 
THETA£--TH£*P  I  .•■'  i  S  0  d . 
CSi=CaS':.THETAl> 
SHi=SirKTHETAl> 
CS£=COSv  rKETA£:' 
Sri£=£IN<THETR£) 
CALL   ZERO    ipPs 1£) 
CALL    NULl.fiAT    CSKJTj  1£?  1£> 
SK.JTc:7<7>=CSl_ 
SKJl  (I?  i;'==SK.JT(7!7> 
SKJT<:73  8>"SH1 
SK  JT  ■;  1  •  £  >  =SK. JT  <  7  J  8  ;■ 
SKJT'::8j7>"-SM1 
SKjT'::£- i>=SKJT<8?7> 
SKJT<8!8::'==CS1 
SKJT(:£5£:>=Sk:JT(8j3::i 
SKJT(55  5>=CS£ 
SKJT<3!3>=SKJT<55  5::' 
SKJT(:5!6>=SH£_ 
SKJT(35  4;'-=SKJT'::5j  6; 
SKJT';6!5)=-SN£ 
SK  JT  <:  4 !  3  :>  =SK  JT  <:  fc  ?  5  > 
SKJT<:6.6:>=CS£ 
SKJT'::4!4)==SKJT<:6i.6> 
DD   £5    1  =  1  ?S 

SUM=0. 0 

DD   £0   K=1?S 

SLIM=SUM+SKJT-::  I  j  K>*iL(K> 

ppcn^suM 

CDHTINUE 

I  i=-PP(:i  )-PPC3;+PP<;5::'H  PP';? 

T£=-PP<£>--PP';4>+PP':;6:>  +  PP<3 

CALL    REAIiNS    CI «  TEMP?  17.  KEL 

FIXH=TEMPa> 


STR 

32 

STR 

o -^ 

STR 

34 

STR 

35 

STR 

36 

STR 

37 

STR 

38 

STR 

39 

STR 

4i:i 

STR 

41 

STR 

4£ 

STR 

43 

STR 

44 

STR 

45 

STR 

46 

STR 

47 

STR 

48 

STR 

49 

STR 

50 

;:TR 

54 

iTR 

55 

JTR 

56 

JTR 

5? 

;tr 

58 

::TR 

59 

:;TR 

60 

sTR 

61 

;tr 

62 

;TR 

63 

:TR 

64 

;TR 

65 

;tr 

6  if' 

:TR 

6-7 

:TR 

68 

:TR 

69 

:TR 

70 

:TR 

71 

:TR 

7£ 

:TR 

73 

:TR 

74 

:TR 

75 

:TR 

76 

:TR 

77 

;TR 

78 

:TR 

79 

■:TR 

80 

;TR 

81 

TR 

82 

-JR 

83 

TR 

84 

TR 

85 

TR 

86 

TR 

87 

TR 

88 

TR 

89 

TR 

90 

TR 

91 

61 


FIXS=TEriP<£> 
EH=TEMPC-!; 

ES=iEfip(:4:J_ 

PN=  u .  5*  T  £*  EN  ■■  Si-'HN 

PS=  0 .5*71  SFS.-  SPAN 

TEMPa>=TEnPa>+PH 

TEMP(:S)=TEMPC£>+PS 

DD  30  K=5j  i£ 

TEriP(:K>=!J(K-4> 

TEMPa3>=PS 

TtrtHp  /  ;^  4  :,  =:pfJ 

TEhP<15>=FIXS 

TEMP<:i6:>=f"i>;H 

CALL   MRITMS    a  s  I  Eht-'»  li-'.  KEL> 
GD   Tu    iOu 


TYPE    Til    ELEMENT 


I  =  IX(l<KEL:i 
J=IW':£?KEL.J 
KK=IM<:35KEl:' 

H=iyx6,K.u.y 

t1K=ITYP-£ 

GO  Tn   <A(!^45f50')3    MK 

Ii=£*I-i+NT3 

I3=£*i<H:-l+HT3 

I4=:£SL.-i+NT3 

I5=£*M-i+Ni3 

I6=£*N-1+NT3 

GD   TD   55 

Ii=3*I-£ 

I£=-3*J-S 

13=£*KK-1+HT3 

I4=5*L-S 

I5=£*M-1+NT3 

I6=£*:N-1+NT3 

GD    TD    55 

I1=£*I  -l-tNT3 

r£=£*.J-i-tNT3 

I3=3*KH::-£ 

I4=£*L-i+HT3 

I5=£«H-i+NT3 

I6=e*H-i+NT3 

CDNTINUE 

AA=iHREAr'-!<KEL> 

yc='::y(:i::'+t(;j)+y(kk>>-^:" 

IF    vYC=G7=H£>    GD  TD    ii 

CALL.   REHTinS    <  3 ?  DMAT  ?  8' 

riii=riMAT<i> 

ril£=DMAT':£> 

Di3=riMAT';3> 

D£1=LMAT(:4;' 

D££=Dr'tAT':;5> 

ri£3=Iif1AT  <  6  ':> 

D3i=DHAT(7> 

D3£=DnAT(8> 


:TR 

9£ 

;TR 

93 

?TR 

94 

;tr 

95 

::TR 

96 

3:TR 

9? 

;tr 

98 

STR 

99 

3TR 

100 

-■T  p 

101 

3TR 

10£ 

STR 

1  03 

STR 

104 

::TR 

1  05 

STR 

1  Ob 

iTR 

10? 

3TR 

108 

STR 

109 

STR 

1 1  0 

STR 

111 

STR 

11£ 

STR 

113 

STR 

114 

STR 

115 

STR 

116 

STR 

117 

STR 

118 

STR 

119 

STR 

1£0 

STR 

1£1 

STR 

1££ 

STR 

1£3 

STR 

1£4 

STR 

1£5 

STR 

1£6 

STR 

1£7 

STR 

1£S 

STR 

1£9 

STR 

130 

SIR 

131 

STR 

13£ 

STR 

133 

STR 

134 

STR 

135 

STR 

136 

STR 

137 

STR 

138 

STR 

139 

STR 

140 

STR 

141 

STR 

14£ 

STR 

143 

STR 

144 

STR 

145 

STR 

146 

STR 

147 

STR 

148 

STR 

149 

STR 

150 

STR 

151 

62 


D33=DMRT 
MUEl=r:Mtt' 


NUL£=DMBT'::ii::' 

STR 

154 

;K£i=xcj>->;a:) 

STR 

155 

Ki3=Ma:i-X(;KKj 

STR 

156 

;>i3e=X(:KK>-K>::j> 

STR 

15? 

Y£3=VU::'-Y(:KK> 

STR 

158 

Y3i=Y<KK>-V<n 

STR 

159 

Vi£=Ya>-Y';j> 

STR 

160 

U1=R<I1> 

STR 

161 

U£=R(:i£::' 

STR 

16£ 

Li3=Ra3> 

STR 

163 

U4=Ra4) 

STR 

164 

U5=K'::i5;' 

STR 

165 

Ub=Ra6:; 

STR 

166 

Vl=R(Ii+l> 

STR 

167 

V£=Ra£+1> 

STR 

168 

v3=R';i3-i-i::' 

STR 

169 

V4=R<  14-1-1  > 

STR 

17  0 

V5=R<I5-H> 

STR 

171 

V6=Ra6-i-n 

STR; 

17£ 

J 

STR 

173 

KDDE^O 

STR 

174 

IF  (KEL.  GT.  MTiE.  hiND,  KEL.  LE.  •; 

HTl£-HNTl.-£>>  GD  TD  60             STR 

175 

GO  TO  65 

STR 

176 

60 

t-TfYl  WU=i  .  -i;"TH 

STR 

177 

E:ThFlIUR=1„--£.*ETH 

STR 

178 

ETHS=ETHTlja*ETHTWD 

STR 

179 

ETP4=--ErttFnUK*ETftFDUR 

STR 

180 

KufiE=i 

STR 

ISl 

GD  in  70 

STR 

ISE 

65 

MT^IMOS-KFL) 

STR 

183 

ETr"-iFniJR=i .  0 

STR 

134 

ETHTWn=ETFiFDUR 

STR 

185 

ETh4==1.i:i 

STR 

186 

EThS=ETr4 

STR 

IS? 

2 

STR 

1S8 

I: 

geherhTE  ahd  store  stpfss 

iHfl  STRAIHS  ELEMENT  BY  ELEMENTC    STR 

189 

I; 

STR 

1 9  0 

?0 

K=KlL 

STR 

191 

F 1  =  ( 4  ■  *ZFi  I  ■:  1  '.>  -EThTWD  >  .-•  ETATWD 

STR 

1 92 

F£=<4.  *ZHK£>  -i  .  >.'-ETftFaUR 

STR 

193 

F3=4.*ZHl^3i-ETAFDUR 

STR 

194 

I: 

STR 

195 

': 

S  PRAIMG 

AT  CEMTRCIII!                       STR 

196 

J 

STR 

19? 

STRA I  HX=Y£3«F  1  *U1  -^Y3 1  h»  r£i*U£-t 

l£WF3«U3-!-(Y£3*4.*ZAKe>-i-Y3i«4.*2AISTR 

198 

I 

a>>*U4.■■EfA£+(:Y3i*4.*ZAK3>-^ 

i'l£»!4.«ZAI<e>> *U5.'-ET AFDUR-i- C Y23*4 .  i«ZSTR 

199 

a 

AI  ■::3>+Yl£*4.  *S:AI  C  1  ::■  >*Li6.-ETAF 

iUR                               STR 

200 

STRAINX=n.  5*STRAIN;k.-A(-. 

STR 

£  0 1 

STRA  I  NT=S-:3£*F  i  *V  1  +■/-.  1 3*  F£*v£-i- 

:.£1*F3*V3h-4.  *';K3£*ZAI  '::£>-HX13*ZAI  <:  1  >STR 

£0£ 

1 

>*V4--ETA£-!-4,  *<X13*ZAI  (3::•-^X£5 

♦ZAI':;2::0*V5.-'-ETAFaiJR-i-4.*<X3£*ZAI';:3>-*-STR 

£  03 

E 

XE1*ZAI  C  1  )  ;.'*V6.-ETAFDUR 

STR 

£04 

STRA r NY=  0 .  5*S  1  RA I MY-'RA 

STR 

£05 

[^ 

STR 

£06 

Al=,»:3S*Fi*iJi 

STR 

£07 

AS  :.i  ':'«T  .  "iJ£ 

STR 

£08 

H  :    ■  'i  ,-lI3 

STR 

£09 

i-vi     ! .  *'..■  j;i*ZAI  .::£>^-X1  3*;AI  (  1  > 

'*U4-ETA£                        STR 

£10 

A5=4.*^Xi3«ZAI<3>-^X£i*ZHl•::£> 

'^US-'ETAFDUR                       STR 

£11 
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Ab=4.*(:X;:a*ZAI<3)+X£lsZAia;' .'S^US-'EThiFQUR  STR  SliE! 

Bl=Y23*Fi«Vi  STR  213 

Be=Y3i*F£*v£  STR  £14 

B3=Y1£«F3*V3  STR  215 

E4=4.*(:Y£3«ZAH£>+Y3i*ZAI>::i  >  :-*V4--'ETft£  STR  £16 

B5=4.  *■:  V31*ZHI  c:3:)+Vl£*7AI  (£)  i^VS-'ETAFDUR  STR  £17 

Efc,=:4.*(:V£3*ZAK3::'+Vl£*ZHia):'*V6-'ETAFDLIR  STR  £18 

STRAHKY=Ai+H£+A3+A4+Ai.+A6+Bi i EE+B3+E4+B5+E6  STR  £19 

STRAHXY=Ci„5»«STRAMXY.-Ah  STR  ££0 

STR  ££1 

SIR  £££ 

ELEMENT  STRESSES  AT  CENTk'DID  IN  X  Y  ANT  XY  DIRECT I DNS      STR  ££3 

STR  ££4 

STRESSX=riJ  i*STRAINX+DlS*STRATNY+Iil3«STRANXr  STR  ££5 

STRESSY=ri£l?KSTRAINM+Dt£*STRAINY+D£3*STRANXY  STR  ££6 

STRES;-:Y=D3  l  *STRA  I  NX+D3£*STRA  1  NY+D33*STRANX'i'  STR  ££? 

STR  ££8 

DETERMINE  PRINCIPLE  STRESSES  aGMAln  SIGMAE; SIGMA3  ...  ...  ...   STR  ££9 

STR  £30 

STRESS  l=-<  ST  RESSX+STRESSY>. •■■•£,  STR  £31 

STRESS£=SG!RT<'::c:STRESS;:'-STRESi:Y ;'.-■£.  >**£+STRESXY**£>  STR  £3£ 

SIGMA1=STRESS1+STRESS£  STR  £33 

S I  GMA3='::  STRESS  1 -STRESS  £:■  STR  £34 

S I GMh£=  ■; S I GNA 1 +S I GMA3 >  * ( riUE  1  <  HUE£ > .■-£ .  STR  £35 

STR  £36 

DETERMINE  DCTAHEhDRAL  STRESSES  SGNAnCT.  TAUQCT  ...  ...  ...     STR  £3? 

STR  £38 

SGMAUCT=  •:  S I GMA 1 +S I  GMA£+S I  GM(-:S .)  '3 .  STR  £39 

TAUuCT=(SCRT'::<SIQMAl-iIQMA£>  **£+■:  SIGMA£-SIGMA3>**E+<SIQMA3-SIGriSTR  £40 

Al >**£)>. -3.  STR  £41 

STR  £4£ 

DETERMINE  PRINCIPLE  STRAINS  EPLISNl 5  EPLISNa?  EPLISN3  ...  ...   STR  £43 

STR  £44 

Ei  =  <STRAINX+3TRAINY>--'£.  STR  £45 

Ee=SQRT«(STRAINX-STRr;iNVV/£,  ::.*;«£+STRAHXY**£>  STR  £46 

EPLISNl=Ei+E£  STR  £4? 

EPLISH3=E1 -E£  STR  £4S 

EPL I  SH£=  ■;  FPL  I  S:N  1 +EPL I S '  iJ  >  *  <■.  IN  C 1  +riUE£ !'  ■£ .        .  STR  £49 

STR  £50 

DETERMINE  uCTAHEhDRAL  ::TRAI^^;;  -  EPSNDCT)  GAMADCT  .„.  ...  ...    STR  £51 

EPSNaCT-(EPLISNl+EPLISN£+EPLiSN3>.'-3.  STR  £53 

G  AMAC]CT=£ .  *  <  SQRT  (  <  EPL  I  SH 1  -EPL  I  SN£  :■  **£+  <  EPL  I SN2-EPL I SN3  >  **£+  (EPLSTR  £54 

ISN3-EPLISN1  >**£::0.--3.  STR  £55 

STR  £56 

STORE  STRESSES  AHD  STRAINS  DN  ThPE  1  STR  £57 

STR  £58 

CALL  READNS  <  1  <  TEMP- 17?  KEl!)  STR  £59 

DELSIG1=SIGMA1  STR  £60 

SIGi=TEMP(4;'  STR  £61 

STRESSX=STRESSX+TEMP(1>  STR  £6£ 

STRESS Y=STRESSY+TEMP<£>  STR  £63 

STRESXY=SrRESXY+TEMP';;3>  STR  £64 

SIGMA1=SIGMA1+TEMP'::4::'  STR  £65 

SIGMA£=SIGMA£+TEMPf:5>  STR  £66 

?IGMA3=SIGMA3+TEMP<6>  STR  £67 

SGnAnCT=:SGMAuCT+TEMP  ■■;  7  >  STR  £68 

TAUDCT=TAUuCT+TEMP'::S:>  STR  £69 

STRAIN;":=STkAINX+TEMP':9,'  STR  £70 

STRAiNY=STPAINY+TEMP(10>  STR  £71 


64 


STRHrirtY-STRHNXY+TEMpJ:  1 1 ) 
EPL I SN 1 =EPL I SN£+1 EMP  a,  S  > 
EPLISN£=EPLISH£+TEMP<]3> 
EPL I SM3=EPLI SHS+TEMP <  j  4 > 
EPSNai::T=EPSHDCT+TEMP';  15;' 
GAfiHDCT=GftMHnCT+TEMPi::  i  6) 

TEMPa>=STRESSH 

TErlP<t:>=S!KESSY 

TEMP<3>=STRESXY 

TEMP'::4>=SIGMA1 

TEMP<5>=SIGMtt£ 

TEr1P(:6:)=SIGMA3 

TEr'ip(;7;'=sGNHncT 

TEriP<8;'=7AUuCT 

TEriP<9:J=STRAIMM 

TEnPt:iu>=STRHlNY 

TE:MP(:il)=STRHHXY 

TEMP  ■;  1 S  > =EPL  I  SN  1 

TEMP'::i3>=EPLISHe 

TEMPa4;'=LPLISN3 

TEMP<15>=EPSHnCT 

TEMPa6>=Gi4rlADL:T 

ChLL    WRITMS    ajTEMp5l7jKEL> 

IF    ':NuTENSfi.LE.  ij)    GD    ID    luO 

hu- TEN-T  1  OH  HNrtLYS  I S 

S I GTDTL^S  [  Gl  +I1EILS IQ i 

IF    '-SIGTlI TL+GrpIACMT;' >    FSslUU-lUU 

WR I TE    (:  6 !  i  65  >    KEL  5  S I G  i  ;  DELS  I G  i 

CALL  ZERD    (PP.-  liE:> 

CALL   HULLMAT    <:SKJT5  1£;  1  £) 

RATin=l .  -AESfSIGl  >--AE?  CCELSIGl  > 

riD=ll 

DD   85    IM=lf 1£ 

DD    80   .jM=IMj1£ 
Mn=ha+i 

Sh  :JT<IM!.Jn)=SKJT<Jf' 5  IM>=jii  iATCMD;) 
CDNTINUE 

ija>=ui 

u(£;'=Li2 
U<3>=U3 
U'::4>=U4 
LK5:;'=IJ5 
Ui;6>=U6 
Ij,r7>=yi 
U'::S>=VS 
U'::9>=V3 

U'::io;'=v4 

U'::ii>=V5 
Ua£>=V6 
DD   95   JI-1.1£ 

SuM=0=  U 

DD    90   KI=ij 1£ 

SL1M=SUM+SKJT<  JI  ?  KI  >*LKKI ) 

PP'::Ji;i  =  -£UM*RATin 
CDHTIHUE 


STR 

£?£ 

STR 

£7^ 

STR 

£74 

STR 

£75 

STR 

£76 

STR 

£77 

STR 

£78 

STR 

£79 

STR 

£S  0 

STR 

£81 

STR 

£8E 

STR 

£83 

STR 

£84 

STR 

£85 

STR 

£86 

STR 

£87 

STR 

£38 

STR 

£89 

STR 

£90 

STR 

£91 

STR 

£9£ 

STR 

£93 

STR 

£94 

STR 

£95 

STR 

£96 

STR 

£97 

STR 

£98 

STR 

£99 

STR 

300 

STR 

301 

STR 

Sfja 

STR 

303 

STR 

304 

STR 

305 

STR 

306 

STR 

307 

STR 

3  OS 

STR 

309 

STR 

310 

STR 

311 

STR 

31£ 

STR 

313 

STR 

314 

STR 

315 

STR 

316 

STR 

317 

STR 

318 

STR 

3 1 9 

STR 

3£0 

STR 

3£1 

STR 

3££ 

STR 

3£3 

STR 

3£4 

STR 

3£5 

JTR 

3£6 

JTR 

3£7 

STR 

3£8 

STR 

3£9 

STR 

330 

STR 

33 1 
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;                              iiD-TEr 

^^I^i^  At 

,ALYSI 

HSTHRT=NTl+i 

Du    155   KEL=NSTbRT 

NTlE^e 

riT=i;s'.:8.KEL;' 

FR<£::'  =  0=  0 

FRa>=|-R'::2;' 

CALL   READMS    <i 

TEfiP<  i; 

.KEL> 

CALL  READnS   a 

TENP! IT 

?  KEL+ 

PMDRML=<TrplPa 

+TEr-iP'::; 

>>*i:i. 

PSHEhR='::TlMP':;8 

+TEtiP<f 

::'>«0. 

STR 
Q<Il>=QaO  +  PPa>  STR 

Q<is>=Qae>H-pp'::s>  str 

y.(.I3 >  =Q c:  I S > +PP <: 3 >  STR 

Q  (14 ;,  =Q  t;  1 4  ;>  +pp  /  4  >  STR 

Q  <  1 5 ;'  =Q  CI  5  >  +PP  i  5  >  STR 

Q  '■.  I  b  J  =Q  <  1 6  > +PP  <  6  >  STR 

Q<Il  +  l>=Qai  +  l>+PP'::7>  STR 

G!';is+i::i=Q';ie+i::'+pp<8::'  str 

QC 13+1 :5=Q<  13+1  >+PP'::9;:'  STR 

Q <  1 4+ 1  >  =Q •;  1 4+ 1  y  +PP <:.  1  LI  >  STR 

Q<I5+l>=Qa5+n+PP<ll>  STR 

Q  < 1 6+ 1 ?  =Q  a  6+ 1 >  +PP  a  £  >  STR 

IFLAG=i  STR 

STR 
CDNTIhUE  STR 

IF    aiTS.LE.O)    GQ   TD    160  STR 

IF    ariTER,LE,iJ>    GD   TD    Ibii  STR 

STR 
?  IMTERACTIQH  ELEMEHT  STR 

STR 
STR 
STR 
STR 
STR 
S:TR 
STR 
STR 
STR 
STR 
IF    CHBSCPHDRMD.LT.GAMACMTir:'    GD    TD    155  STR 

IF  <pndrml+gama(:mt;o  iiojIossIOS  str 

IF  aiCYCLF,EQ,l>  GD  TD  155  STR 

RAT  I  u= AES  <-  PSHEAR.-pNDRriL )  STR 

IF  '::RATin.LE.£LIP>  GD  Tu  155  STR 

Rl=HES'::TEr'1Pa5).-TEMPa6::0  STR 

QQ=TEPiP(i3>^'TEriP';i4>  STR 

FIXH=TEMP<16>  STR 

PN=TEMP<i4>  STR 

DR=SLIP-Rt  STR 

T3=QQ_R1  _r|j^  STR 

DSR=DR*FIXH.'-T3  STR 

T4=<PH-DSR:'.-PH  STR 

FRa>=T4+FRa>  STR 

FR(l>=-FR';i>  STR 

WRITE    i::6?i7u>    KELj  RAT]  CJi.  SLIP;  Rl  -  QQ- DP- DSR;  FIKMjPMjFRCD  STR 

Rl=ABS'::TENP<i5>.''TEMP<16;'>  STR 

QQ=TEHP  ( 1 3  >  ---TEMP  a  4  >  STR 

FIXH=TENPa6>  STR 

PH=TENP'::i4>                                                                '  STR 

DR=SLIP-Ri  STR 

T3=Qi:!-Rl-IiR  STR 

DSR=riR*FINrVT3  STR 

T4='::PH-riSR>.- PH  STR 

FR<e>-T4+rR':£::'  STR 

FR'::£;'=-FP',2:;'  STR 

KEL1=KEL+1  STR 

WRITE    Cfa-lrO    KELI-RAIIDjSLJP;.  Rl.QQ7DRjriSR.riXN«PH!FR<£>  STR 

GD    TD    115  STR 

STR 
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110 


115 
ic'u 


135 

14!] 


^;lLjKEL J  jFRU  .'I'FR':: 


iHU 


Tl  =  (TE:riP(16>+TEhP(16); 
T£:=<TEMP'..14)+TE;HP';:14):: 

FR<i>=F 

T3-1  =  -T1--I£ 

FR(S>=T3.--£. 

FRa>=FR(£> 

KEL1=KEL+1 

WRITE    <b>ir' 

NUM=KEL 

ND=1 

CALL   REFUifiS    (7 J  DMA! 

Mu=U 

Du  i£5  IM=ls8 

DD  i£5  JM=i5  8 

Mu=f'iD+i 

SKJT  <;  I  f1  ?  JM  >  =DMftT  < MCi  > 
CDNTIHUE 
DD    14u_IM-li.S 

GD    !D    'MSiij  135>!    riu 

iJ<In>  =  TEnp(:if'!+4;) 

GD    TD    i40 

Li(IM>=Tt:MP(iri+4> 
CnMTIHUE 

nn  150  J I -1.-8 

SLiM=i:i=  0 

Du    i  45   i<  I  =  i  J  S 

SUn=Sufi+SKJT<  JI !  KI  >'*Ui:'KI ! 


-ABSCm'-'ABSc;" 


FR'::Ma; 


LiM*FR';Ha> 

o>  (ppcji: 


PPUI 
CDHTUiUE 
URITE  (6«  iS 
Ii  =  lH(:NUh1j  i 
I£=IH(:NUM3£ 
I3=IA(NUri'.3 
I4=IAa-!UMn4 


Q<ii>=Qai>+ppa:) 

G!';il  +  l>=Q';Il  +  i)+PP'-£> 

Qa£>=Q<I3>-ipP<3> 

Q  f  I  £+  i  ;■  =Q  <  I  £+ 1  >  +PP  ( 4  > 

Qa3;'=Q<I3)+PP(5> 

Q<  13+1  ::'=Q<  13+1  :j+PP'::6> 

G!i:;i4::i=Q(:i4>+PP(7) 

Q  <  1 4+ 1  >  =Q  <  !  4+ 1  >  +PF'  <  S  > 

IFLAG=1 


Nn=HO+ 1 
IF  CND.GT. 
HUn=KEL+l 
QD  TD  1£0 

155  CONTINUE 

160  WRITE  <6j185: 
RETURN 

165  FDRMAT  (lOKn 

IH  DELSIGMAl  = 
IrO  FDRMAT  <1Xj 


jQ  TD  is; 


£4HTENSIDN  IN  ELEMENT  NO  =  > I5t 
-    1  "^  1  P  3  i 
SHI NT  SLIP? 


3H   SIGMAl=)Elc 


STR 

392 

STR 

393 

STR 

394 

STR 

395 

STR 

396 

STR 

397 

STR 

398 

STR 

399 

STR 

400 

STR 

401 

STR 

40E 

STR 

4  03 

STR 

4  04 

STR 

405 

STR 

406 

STR 

4  07 

STR 

4  08 

STR 

4  09 

STR 

410 

STR 

411 

STR 

41£ 

STR 

413 

STR 

414 

STR 

415 

STR 

416 

STR 

417 

STR 

418 

STR 

419 

STR 

4£0 

STR 

4£i 

STR 

4£c 

STR 

4  £3 

STR 

4£4 

STR 

4£5 

STR 

426 

STR 

4E7 

STR 

4£8 

STR 

4£9 

STR 

430 

STR 

431 

STR 

43£ 

STR 

433 

STR 

434 

STR 

435 

STR 

436 

STR 

437 

STR 

438 

STR 

439 

STR 

44  0 

STF? 

441 

STR 

44£ 

STR 

443 

STR 

444 

STR 

445 

STR 

446 

STR 

447 

STR 

448 

STR 

449 

STR 

450 

STR 

451 
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i?5  FDRHAT  ax?  SOHJUTERHCTinh   TEM;iaH,£I5=£El£.3)                  STR  45£ 

ISiJ  FDRMAT  aX;   9H I  NT.  FORCE  SEia.J)                                                                                       STR  453 

iS5  FORMAT  (IX?  PSHaVERLAY  M?0>  CarlPLETED)                         STR  454 

STR  455 

END      _                                                      STR  456 

D'v'ERLht-::RDY?  J?  0>                                                 RES  1 

D VERLft'i'  C  RD Y  5  5 .  0  >                       RES  £ 

PROGRAM  RESULTS                                               RES  3 

CDMMDM  NMDDES5  HELEMMT; HUSIF; MEAHS?  MB?  MT3= ISTOP5  MCYCLEs  LAYERS? ISTEP?  RES  4 

1 NSTEP ,  MT  i  £  >  ET  A !  MT 1  s  MT£  5  MLl  TEMSM :.  I  FLAG  ?  MS  I ZE  s  MCDDE  <  55  0  >  »  X  C  55  0  >  j  Y  C  55  ORES  5 

2:)!JMriX(:51>!  AMLSIS?  IX(:S?£!:50>5ARFAA'::£50:)!  IMIiXCESO:;?  IMDEXC£50>»GAMAi:£RES  6 

35jjZAI(3;'                                                    res  7 

CDMMDM  .-■l.--  Ei.NME.pHriIUS..::CEM.Yi:i:M>EIjKH.KS.Hl'H£,  INTER             RES  8 

CDMMOfi  -S'  R                                                     RES  9 

DIMENSIDN   tempi:!?:'!    RUlnOj    RPfllOu)                                                                          RES  1  Ci 

DIMEMSIDN    EaCO?    MUEaO)                                                                                                           RES  11 

DIMENSIuN    TMP'::3b>?    U<6>-    SRG<6-6'}f    PXYCSO;'!    QXYOOh    MDMOOj    POORES  1£ 

i>!    Gk5u)!    PP'ib;'                                                                                                                                  RES  13 

DIMEMSIDN  SQ<3u>j  Sl<9u>                                         RES  14 

REAL  HUE                                                     RES  15 

REAL  KN.KS                                                      RES  16 

REAL  MOM                                                        RES  1? 

INTEGER  UliU,":5U3?U43U5!ll:.5U7?US':U9                               RES  18 

RES  19 

THIS  ROUTINE  RECORBS  RESULTS  FROM  DIFFERENT  ROUTINES  AND  TAPES  ANDRES  £0 

STORES?  ADDS?  FDR  SUCCESSIVE  TERMS  IF  INCREMENTAL  ANALYSIS  AND 

F I  HALLS'  DUTPUTTED 

URITTIHG  HuDAL' FORCES?  lOSPLACLnEHTS?  ELEMEHT  STRESSES?  STRAINS  » 

REWIND  £ 
HDRIHG=3«HT3 
DO  5  1=1? NCR IMG 
5  Si'::D=R'::i) 

READ  •;£;'  CRRd;'!  1=1  ?MDDF- 
DO  10  I=l?NDnF 

R<i)=Ra>+pRa  J 

10  COMTIHUE 
REMIND  £ 
WRITE  <£>  (R'- li'?  I  =  1?NDDF'- 

WR I TE  < 6  ?  9 0  >  I STEP  ?  NCYCLE 
IF  a^T3.LE.0;.'  GO  TO  £0 
WRITE  '::fc.!95> 
DO  15  1=1? NTS 

Kl=3*I-£ 

WRITE  Ceplim,'  I.R'::K1>"RCK1  +  1j?R<K1+£> 
15  CONTINUE 

£0  IF  ^HCYCLE.  ME,  NSTEP)  GO  "1 Q  35 
IF  aPLAG.GT.  Ci;'  GD  TD  35 
WRITE  <&!  i  05-' 
NT3Pl=MT3+i 
DO  £5  I=NT3P1--NNDDES?3 

Ki=£*I-i+MT3 

K£=K1+1 

K3=£*a  +  i>-l+NT3 

K4=K3+1 

K5=£«';;i+£;'-i+HT3 

K6=K5+i 

I£=I+£ 


RES 

£1 

PES 

££ 

RES 

£3 

RES 

£4 

RES 

£5 

RES 

o^_ 

RES 

£7 

RES 

£S 

RES 

£9 

RES 

30 

RES 

31 

RES 

3P 

RES 

33 

RES 

34 

RES 

35 

RES 

36 

RES 

37 

RES 

38 

RES 

39 

RES 

40 

RES 

41 

RES 

4£ 

RES 

43 

RES 

44 

RES 

45 

RES 

46 

RES 

47 

RES 

4S 

RES 

49 

RES 

50 

RES 

51 

RES 

5£ 

RES 

53 

RES 

54 

RES 

56 

68 


WRITE  Cfc-llO)  I?R(:K1>-Ra(£)-  U  ?  R<K3>  <  F:'::K4::' ■>  IS,  RCKS)  -  R'::K6> 
£5  CDHTINUE 

PR:  NT  ELLflEHT  STRESSES 

URITE  tfaniis:;' 

NTl£Pl=HTie+t 

DD    Su    I=nTi£Pl!tiELEMNT  . 

CALL  REftLMS  (IHEMP?!?;!) 

WRITE  CbjlSu)  (I?  (TE(1P':;j)?J=l53>? 
3  0  CDHTINUE 

CALCULhTICIN   of    NDRMhL;.    shear    FDRCE    AND   MDMEHTS   FDR    PIPE   HDIIES 

S5  IF  ';nt5.le.o:>  gd  to  ss 

WRITE    ';t.!i£5;' 

HEL=NT3-i 

Du   55    I=1?NEL 

IF    <IK'::7?  O.NE.l)    GO    TD   55 

ii=iKa!i:) 

I£=IX<£! I' 
Ul=3sai-l>  +  i 

Uc:=ui+i 

U3=U 1 +£ 

U4=3*(:ie-i>+i 

U5=iJ4+l 
Ut.=U4+£ 
U'::i>=Sl(Ul> 
U<£:i=SI<U£> 

u>;3>=siai3::i 

1J';4::'=S1'::U4> 
U';5>=SiaJ5> 
U';6)=Si<U6> 
CALL  REhUMS  '::7jTMPj36.  I> 

MD-0 

CALL   REAIiMS    <.  1  -  TEMP  ?  1  ?  -  I  :■ 

nn  40  K=i?fc. 

DD   40   L=i^6 
MD^MD+l 

SRG':K!L>=TriP(MD) 
40  CDHTINUE 

nn    50    K=l5  6 
SUM=0. 0 
DD   45   L-1.6 
45  SUM=SLiM+SRG  <  K  ■■  L  >  «U  C  L  > 

pp<K>=SnM 
50  CDHTINUE 

PXY< 1 1 >=FP<4>+TEHPC3> 
QXV';ii:J=PP(5)  t-TEMP'::4> 
MDMai>=pR<6>+TEMP<5:i 
PHYt;  I£>=PP(  1  >+TEMP<6> 

exvae>=pp'::£>+TEMPC7::' 

MDM<I£>=PP(3::'+TEMP'::8> 

tempo:)  =PftY<ii> 

TEriPC4>=QXYai> 
TEMP'::5:)=MDM(Ii::' 

temP'::6:j=p:^y<I£) 

TEMP':;7>=Q:<Y<ie:;' 
TEMP(8>=Mf]lia£) 


RES 

57 

RES 

58 

RES 

59 

RES 

60 

RES 

61 

RES 

6£ 

RES 

63 

RES 

64 

RES 

65 

RES 

66 

RES 

67 

RES 

68 

RES 

69 

RES 

70 

RES 

71 

RES 

7£ 

RES 

73 

RES 

74 

RES 

75 

RES 

76 

RES 

77 

RES 

78 

RES 

79 

RES 

80 

RES 

81 

RES 

82 

RES 

S3 

RES 

84 

RES 

35 

RES 

86 

RES 

87 

RES 

88 

RES 

89 

RES 

90 

RES 

91 

RES 

9£ 

RES 

93 

RES 

94 

RES 

95 

RES 

96 

RES 

97 

RES 

98 

RES 

99 

RES 

1  0  0 

RES 

101 

RES 

loe 

PES 

1  03 

RES 

104 

RES 

105 

RES 

106 

RES 

107 

RES 

1  OS 

RES 

109 

RES 

110 

RES 

111 

RES 

HE 

RES 

113 

RES 

114 

RES 

115 

RES 

116 

69 


CALL   WRITi'Vi    OfTEMpjl,  jI> 
55   CDNTIMUE 

DD   70    1  =  1  >hT::! 

IF    <I.EQ.HT3>    GO   TO  60 

GD   TD  65 
60  PXY(I>=-PH¥a> 

«;■■■;¥<  I  )=-Q«Y>;  I  > 

MDMa>=-MaM<I> 
65  YC=V<I>-'t'CEH 

IF    CXC.EQ..  0.  0>    Xi::=0„  Ui  001 

TH='i'i::.--Xi:: 

THETA=AThiH';;TH> 

SM=NlN'::THErft> 

CS=CDS(THETh> 

PC  I  >  =-PXV <  I  > *SH+QXV ■;  I  >  ^cs 

Q<  I  >=-PXVa>*r:S+QXV(I  >*SH 

WRITE  <6?i30>  IjP'::i>!Q(I)?MnMi:  D 
70  COHTIHUE 

URITE  ';:65l40'^ 

NSTftRT=HTl-H 

m   75    I=NSTART=NTie5£ 

CALL    REHliMS    <  1  <  TEMP;  17?  I> 

Ti=TEMP(r^ 

T£=TEMP'::£> 

CALL   REALHS    a?TEi1Pjr,M  +  0 

TEMP'::i>  =  u=5«';Ti+TEMP(i  >> 

TEMP  ■::£;'  =  0.5#<T£+TEMP  <£■:■) 

RAT  I  D=HES  f  TEMP  (  £  >  .--TEHP  (  1  >  > 

IP1=I+1 

WRITE  Cb.lBS)  I.IPljTEMpf  1)!TEMP';£>!RATID 
75  CDHTIHUE 

PPIHT   Fl,  rriEMT   STRAINS 

WRITE    Ce^MS,' 

DD    80    I=NTl,E:Pl.hELEHHT 

CALL    REAIiPiS    <  i  j  TEMP  s  1  ?  ?  I  > 
WRITE    <6-l£0>    <l!  <TEMF(j;'!J=9)  i6>> 
80   CDHTIMUE 
85   WRITE    '.h,  150::' 
RETURfH 


14HRE 


15?  1.7H    LAYERS  DF  FILL? 


lo;-:?  SHMnnE 
shrdthtidn;  ■•-: 


SHHDDE  HD. 


90  FORMAT  (IHi^ 

IHINCREMEHT  NG  =  ?  I5.-v> 
95  FDRMhT  aOX^  £SH3  -  ruO.F^  liDIiAL  SEFLECTinr 

li>  i4HX  -  IiEFLECTiaN?5W5  i4HY  -  OEFLECTIDMj  !: 
100  FORMAT  ■:;7XjI5.7>;.3<E14.4«5X>> 
105  FORMAT  C .-■.•'■  1  OX .-  SlHMDnE  PDIHT  DE;  LECT I DN?  •■■'.'■:: 

1  HX-DEFLlCT  I DH .  :3X .  1  eHV-DEi-  LECT I  OH ;.  .•'.■■ ;. 
iiO  FORMAT  r:3i::5Xn  I4i5X!Ei£.4!  3XjE1£,.  4::0 
115  FORMAT  (IHl-'lOX!  I6HELEME1T  STRESSES?  .■■'/-SX?  IIHELEMEHT  ND.sSX 

1TRESSX.7XJ   7HSTRESSYj7Xj   7HSTPESXY. 7X.   eHSIGMAl-SX-   6HSI 

£Xj   6HSIGMA35  7X5   SHSIGMaJCTj  6X1.   7HTftUnCTA5  .■-■> 
i  £  0  FORMAT  < 5X  j  1 5  ?  gX;  SE 1 4 .  £ ;' 
1£5  FORMAT  ':>'■•-■  1  OX <  47HN0RMAL!.  AND  S  EAR  FORCES)  MOMENTS  AT  PIPE 

l-'lOXi.   SHNOIJE  NO.slOX!.  ISHHDRMAL  FDRCE1.5X.  UHSHEAR  FORCE?  5X 

£OriAL  MOMENT-.  > 
130  FORMAT  f  13X?  I  3?  lEXs  El£.  £<  a'::4Xj  E  i  £.  £;)  > 
135  FORMAT  (3X? I??   3H  +  j IE? 1£X? El u. £, l£X? El  0= £? 1£X) Fl 0. 4> 


RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 

RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 
RES 

5X.  15RES 
RES 
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140  FDRMhT  i:>'.--v-iuK.  seHHDRMRl  AND  SHEAR  STRESS  DM  SDIL  IN  IHTERhCTIDN  RES  177 

ILhVERj.-vSX?   PHElE.ND.?  liiXj  ISHNnRMAL  STRS.jIOX?  ISHSHEAR  STRESS?  IRES  178 

£UM.  ISH  RATIU  T-N   j-O                                           RES  179 

145  FUKflAT  aHl-'iijyj  i5HELEMi:NT  3TkiiINS?.-v5Xj  1 IHELEMENT  NQ,  5  3:K:.   7HSTRES  ISu 

1RAINX.7K!   7hSTRAINYj7Xj   7HSTR)^NXV? 7Xj   7HEPSILN1 j 7Xj   7HEPSILN2j RES  181 

S7Xi.   7HEPSILN3.7X5   7HEP::;H0CT!  7:^   7HGAMAnCT.  •>                   RES  1h£ 

150  FDRMhT  <iX!  S3HDVERLAY  ■-■■;=  0;^  CDNPLETED)                           RES  183 

RES  184 

EMU                                                             RES  185 

U VERLA'f  ■:.  RDY  ?  6  <  U  >                                              SDL  1 

avERLK'.'<Ravj6!i;!>  .  sdl        a 

PRDGRhM    RESDI  V  SDL  3 

CDiIMDH    NHDDES?  HELEMHT!.  NIHjFj  HEAIiil.  NDi.  NT3?  ISTDPs  HCYCLE-  LAYERS.  ISTEP?  SDL  4 

iHSTEp.HTi£jETA.HTljNTSjHriTEMSN>  '  FLAG!  MSIZE?  MCDDECRSiV)  5  X';550::' ?  YCSSOSDL  5 

SiJjJNDJ-XSi)  jAHLSISj  D<';S?£:nij>:i  ARFhH<£5u>j  IHriX^eSO::' j  INDEXCeSfOjGAMACESnL  6 

35>7ZAIC3>  SDL  7 

CDnnuN  ,'3-"  R 

CuMnON    ■■-S--   Q  ■:  i  1  0  0 )  ?  L I  ST  <  ]  1  u  1 ) 

D I  MENS  I DN   A '.  1  03  ?  2  06  'i ,    E  ( c '  06  >  ■>    RpR  AY  <  1  03  ::> 

LI  MENS  I DN    R'  :l  100> 

REWIND    9 

SOLUTIUN    DF    -^^TIFFNESS    EQiiATIDNi    FDR    GIVEN    EDUHDAR 

NUNELF-l 

HL=i 

NM=ND 

KSHIFT=n 

ND£=£*Nri 

HEND=3*NT3 

NriINC=l 

IF    (:NT3.LE.  U)    NDINC=0 

IF  cnts.le.o)  Nri=Nri-i 

CALL  NULLMAT  CAs  ND- ND2:> 

CALL  ZERD  '::E-NLi3> 

NX=NL-1 

CALL    ZERD    < ARRAY. NSIZE; 

DD    iO   N-inNSIZE 

CALL   READciS    <  1  0?  ARRAY.  NSIZE?  uX+H) 
DD   5   M=ljriSIZE 
5  A<N.M::'=ARRAY(M> 

10   uDHTINUE 
GD   TD   3u 
15   HX=NL-1 

CALL  ZERD  CARRAY. NSlZE> 
DD  £5  N=1.HS12E 

CALL  READMS  CI  0? ARRAY. NSIZE. HX+H) 
DD  £0  M=i:.HSIZE 
£0   _  A'::N.M;'=A(N.M::'  +  HRRAY<ri> 
£5  CDNTINiJE 

MfiUiriCATION  FGR  LDAD  AND  BOUNDARY  CDNJ 

3  0  I=NL 

35  IF  (KGT.NEND)  GO  TD  70 
h  I ND=|-  LD AT <  I  +  £ .'■  .'3 . 

ii=inT':;fihd:) 
pX=Q';I;' 
PY=Qa  +  i> 

ICDliE=NCuDE<  [i::'-3 
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GO   TO    ';4Cij45:.  5n?55!t.Uj65   J    ICOIiE 
40   .J=I-KSHIFT 

B<;j+l>=E<J+i:J+PY 
1  =  1+3 
GD  TD  35 
45  J=I-KSHIFT 
P=  0 .  0 

ohll  mdhify  < jjHrHp.-HjE:) 

E<J+l>=E<Jfl  ■+PV 

QD  TD  35 
50    J=I-k;SHirT 

E';j)=B';j>+PM 

P=u.  u 

CALL   MDIUFV    ■:  J+i  ?  HDs  P- ft?  I  ) 

1  =  1+3 

QD  TD  35 
55  J=I-kSHIFT 

B(J>=E<  J)  +  F';; 

E<J+l>=BvJ+i:>i-PY 

P=n.  0 

CALL    MODIFY    ■' J-J-3?  HH?  P?  ft?  T> 

1  =  1+3 

GD  TD  35 
60  J=I-KSHIFT 

P=0.  u 

CALL   MDCIFY    '.J;  rili?  PjA?  B> 

B'::j+i>=B(:j+i>+py 

CALL    tiDDIFY    C  J+E?  HD?  P?  ft?  F> 
1  =  1+3 
GD    TD    35 

65  j=i-ks.hift 
b<j::'=E(j>+p:'' 

P=0.  0 

CALL    MODIFY    ■■  J+1 5  HD- P?  A?  i"' 
CALL   MOD  I F',-'    (  j<d ,  ND  ?  P  j  A  - 1  ) 
1  =  1+0 
GD   TD   35 
7  0    IF    <I.GT.NM.DP=  I.GT.HIiDF'    GO   TF   95 
FINIi=FLDAT':  i-MT3+l  J'-'S. 
II=IHT<FIHD> 

px=Qa:i 

PY=Q.:  i  +  i;j 

ICDIiE=HCGriE  <  il  >  + 1 

GD  TD  <?5<80;.  85?9U>-  ICDDE 

?5  j=i-kshift 
b<:j>=e<.j>+pk 

BU+1>=B(J+1>i-PY 

I  =  I+£ 

GD  TD  70 
SO   J=I-KSHIFT 

P=0.  0 

CALL    MDIUFY    ■:  J.  HD- P?  A?  E  J 

B< J+i>-B<J+r^+PY 

1  =  1+8 

GD  TD  70 
85   J=I-KSHIFT 

B(J>=ECJ>+Pr-; 

F'=0.  0 
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58 
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JDL 
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::nL 
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::DL 
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::DL 

67 

:;DL 

68 

;:dl 

69 
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70 

;:dl 

71 

mL 

7£ 

;:dl 

73 

;dl 
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■QL 

75 

;;dl 

76 

::DL 

77 
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::DL 

79 

;;dl 

SO 

::DL 

81 

;:dl 

8E 

::DL 

83 

iBL 
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;dl 
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::DL 
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JDL 
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;:dl 
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::DL 
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;dl 
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;DL 
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:dl 
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::dl 
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;:DL 
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;:nL 
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;dl 

105 

;DL 

1  06 

;aL 

107 

:DL 
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:DL 
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;dl 
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;DL 
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CALL   MnrilFY    (J+ljNIivPsA.i:' 

I  =  I+£ 

GD  TD  ?u 
90  J=I-kSHIFT 

P=0.  0 

CALL   ilDDIFY    (  J.  NIn  P?  A?  10 

CALL  MDDIFV  .;  J+1  <  HE- P- A.  T:> 

I=I+S 

Gu  TD  70 
95  CDHTIhuF 

REBiiCE  BLOCK  GF  EQUATIGMS 

Du  i  1  LI  ri=  1 J  HD 

IF  (ACi  jN)  =  lG!.  0.  iO  GD  TD  ilu 

DD  105  L=£i.MEAND 

IF  ';h(LjN).EQ.  0=  0>  GiJ  TD  105 

c>;=A<L.N:'.-H<ijH:i 

I=N+L-1 

.J=  0 

DD  100  f;=L'MEANL 

luO       ACJ-  I>-^A<:.Ji  I>-CX'i»A':K!H> 

Ea;5==B<n-ACL.M::'*B0i> 

A'::L)H>=CX 
105    CDHTIHUE 
iiO  CDNTIHUE 

IF  (NM.GE=NDnF>  QD  TD  iEi^ 

m-'lTE    BlDCK   Df     CQUhTiDNS    (:PEliUCED>    OH    TAPE    9 

UR I TE    <  9  )    (  E .:  H  ;■ .  t:  A  ( M  5  M >  3 1  =- 1 5  MBi- ;!  iD  ':■  5  N==  1  ?  m  > 

SHIFT  BLQiK  DF  E'lH.  UP  FDR  HEXT  BLOCK 

Du  115  H=l.tili 

MM=HrH-H 

B';N>=B(Mn:) 

B(MH::'  =  0.0 
DO  115  M=l!MBAND 

A<MiH)=Aa'!?MM:) 

A>::n3nl'i::'  =  0,.  0 
115  CDNTIhUE 

KSHIFT=KSHIFT+HD 

IF  (HUMBLK^EQ.l)  MD=HD-HTIHC 

HUMELK=HUMELK+1 

NM=HM+MD 

NL=HM-Nri+i  I 

GD   TD    15 

BACK  SUBiTITUTJDN  IN  :RUSS  ELIMINATIDN  PROCESS 

leo  CALL  ZERO  '^Rj  HDDF) 

IF  (HUMELK.EQ.  I''    HEIHC=0 

NU=ND*HUHELK+ 1  +Hri  I HC 
HE=HUMBLK 
i£5   DD    135   fl=1jHD 
H=HD+1-M 
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DD  130  K=£!MBAND 
L=H+K-1 
i  3  0    B  '■.  N  y  =E  '■.  N )  -A  ■-'!  K  ?  N  >  *B  ( L ) 
HM=N+Nri 
IF  (:NB.EQ„£>  NM=NM+NDINr: 

B<Nfi:>=B':H;' 

NiJ=riU-l 
135  R';MLf)=Ba<) 

NE=ME-1 

IF    aHE.EQ.l)    HIi=Mri+riDIMC 

IF  CNB.LE. 0>  Gn  TD  140 

BACKSPf^CE  9 

READ  ■;  9  >  ■;  E  ( tl  >  .  ( A  ( M .  N  >  J  M- 1 .  HE  AM  Ii  >  ,  H=  1  ?  HE ) 

BACKSPACE  y 

GO  TD  1.25 
140  CDHTIHUE 

RETURN 


Erin 

SUERDUT I  HE   fit  Jl)  I  FY    ( N  :■  MD  ?  f ' .  A  -  E : 

D I  MENS  I  Dri   A  ( 1  03  <  £  06  >  J    E  <  c'  06  > 
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18£ 
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1S3 
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184 

DL 

185 
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DL 

189 

DL 
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DL. 

191 

•■:; 

4 

5 

MDEIFICATIQN  FDR  BDUriEARV  CnNIflTIDNS 

NIi£=£*ND  7 

DD  10  M-£i.HD  8 

K>H-fi+ 1  9 

IF  '.'K.LE.  U>  GD  TD  5  10 

Bc:t<>-E(i<>-A(MjK>*P  11 

A'::m?k;'=o=o  ie 

5    K=H+M-i  13 

IF  ( K .  GT  =  hri£ >  Gu  TD  1  i                                            14 

B':;K;'=E(iO-A'::M?rO*P  15 

A':M5N)  =  0,0  16 

10  CuMTIHUE  17 

A<l!N;'  =  0=0  18 

E<H>=0. 0  19 

RETURN  20 

£1 

END  pp 
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PROGRAM  PRDPRTY  a HPUT, OUTPUT > 

PRO 

£ 

DIMENSIDN  J^P'::eO>5  YPCSOj  VPCeO'.  YDPlCeCi).  YDPVCSO).  TITLCIO) 

PRO 

3 

REAL  NUEnNUEIN 

PRD 

4 

DhTh  PLSTRS. PLSTRM/6HPLSTRS5  6HPLSTRM.- 

PRO 

5 

PRD 

6 

THIS  PRDGRhM  GEMERhTES  DATh  REQUIRED  FDR  EVALUATIQH  DE  TAHGENT 

PRO 

7 

MODULUS  AND  PDISSDMS  RATIO  VS.  SIGriA<DCT>  FOR  VARIOUS  STRESS  RATIPRD 

3 

PRD 

9 

READ  70)  aiTLa>5l  =  l5lO> 

PRD 

10 

READ  75?  TESTjNDCURVSjRF 

PRD 

11 

IF  <TEST.EQ= PLSTRS)  MG0=1 

PRD 

IE 

il-  <:  1  ES  1  =  lQ=  PLS  1  KhL'  HbD=£ 

PRO 

13 

DO  65  IK=1.HDCURVS 

PRD 

14 

READ  80?  SIGMA35HP 

PRD 

15 

READ  35-  <XPa::'5YP<n?VPa>j  I  =  1jNP> 

PRD 

16 

PRINT  70.  '::TITL(I>5l  =  l5lO> 

PRO 

1? 

PRINT  90?  iy?SIGMA3!NPi.RF 

PRD 

18 

PRINT  95?  <I.XP(I>jYPa>?VPa>,I  =  li.HP> 

PRO 

19 

no  5  1=1. .HP 

PRD 

eo 

vP'::i>=vpa>-:>^pa> 

PRD 

£1 

S1IN=SIGMA3 

PRD 

££ 

XIH=0. 0 

PRD 

£3 

E3IN=0.0 

PRD 

£4 

CALL  SPLINE  CNP.  XPj YP? YDP1> 

PRD 

£5 

CALL  SPLINE  CNP.  XP?  VP.  YDPV:) 

PRD 

£6 

PRINT  100 

PRO 

£7 

X=0. 00001 

PRO 

£S 

DX=0. 000001 

PRD 

£9 

HUE=0.5 

PRD 

30 

DO  15  L=£..NP 

PRD 

31 

IF  <x,lt.xP'X>.and.x.ge.xp<:l-i>>  gd  to  so 

PRD 

3£ 

CDHTIHUE 

PRD 

33 

GO  TO  65 

PRD 

34 

A=YDPl<L-n 

PRD 

35 

B=YDP1(L::' 

PRD 

36 

c=xp(:l-i;' 

PRD 

37 

D=XPt:L> 

PRD 

3S 

e=yp(;l-i> 

PRD 

39 

F=YP<L> 

PRD 

40 

PP=X 

PRD 

41 

si=drdinet<a.b.c.d-e.fpp> 

PRD 

4£ 

PRD 

43 

DO  £5  L=£:.NP 

PRD 

44 

IF  c:x.lt.xp(:l>.and.::.ge.xp<l-i;'>  go  td  so 

PRD 

45 

CONTINUE 

PRD 

46 

GD  TD  65 

PRD 

47 

A=YDPV<L-i:) 

PRD 

48 

B=YDPV':;L) 

PRD 

49 

c=xp<:l-i> 

PRD 

50 

D=XPf;L) 

PRD 

51 

E=VP(L-1> 

PRD 

5£ 

f=vp(:l> 

PRD 

53 

PP=X 

PRD 

54 

E3=DRD  I  NET  <  A  ?  B .  C  <  D .  E  J  F  i.  PR  > 

PRD 

55 

E1=S1+SIGMA3 

PRD 

56 

B3=SIGMA3 

PRO 

57 

GD  TD  ':;35.40>j  NGO 

PRO 

53 

B£=E3 

PRD 

59 

bO  iD  45 

PRD 

60 

75 


40    B£=HUE*<Bi+E3) 

PRO 

61 

45    SGr1ADCT=  <  B 1 +E8+E3  >  -^i . 

PRO 

6£ 

A1=B1-BS 

PRD 

63 

AS=B£-B3 

PRO 

64 

h3=B3-B1 

PRD 

65 

T  AUDCT=  0  =  333*S0RT  c:  A  i  **2+ A£**e+ A3**2  > 

PRO 

66 

Ti=X 

PRD 

67 

T3=E3 

PRD 

68 

GD  TD  <:50?55>!  NGD 

PRD 

69 

50    T£=T3 

PRD 

70 

GD  TO  60 

PRD 

71 

55    T£=0-.  0 

PRD 

7£ 

60     A1=T1-T£ 

PRD 

73 

A£=T£-T3 

PRD 

74 

A3=T3-Ti 

PRD 

75 

GAMAnCT=£,  WCSQRK  Al**£+A£**e+A3*»£>  >--3. 

PRO 

76 

TAUF=RF*i.SGMADCT 

PRD 

77 

RATID=TAiJaC:T-'TALlF 

PRD 

78 

EP3=E3 

PRD 

79 

T 1  =S  I  GriA3*X-E  1  *EP3 

PRD 

SO 

T£=  •;  S I  Gri  A3+B  i )  *  (:  X-EP3  :> 

PRD 

81 

HlJE=Ti.'T£ 

PRD 

82 

EMi=Ei*'::i.-NUE**£::'^X 

PRD 

33 

Er1£=NUE*  ■;  1  =  +riUE  >  *S  I  GMA3/'X 

PRD 

84 

EM=EM1-EM£ 

PRD 

85 

EM=iOO.«Er1 

PRD 

86 

DSi=Sl-SiIN 

PRD 

87 

riEPi=;<;-xiH 

PRD 

88 

riEP3=EP3-E3IN 

PRD 

89 

ti=dsi--';depi-dep3> 

PRD 

90 

T£=DEP3.'-  <  DEP 1  -DEPS  > 

PRD 

91 

TnnD=T1*a  =  -T£> 

PRD 

9£ 

ThDD=iOO,*TnDri 

PRD 

93 

TANH1JE=-DEP3.'-  ■;  DEP  1  -DEF  3  ) 

PRD 

94 

PR  I  NT  1  05 ,  X  J  E3  J  S 1 ;  TAUaCT  ?  GAMADCT  j  HUE  ■>  EM :.  TAtiNUE ;  SGMADCT .  TMDD 

! RAPRD 

95 

1   TID 

PRD 

96 

XIH=X 

PRD 

97 

SiIH=Sl 

PRO 

98 

E3IN=EP3 

PRD 

99 

rix=i.i*rix 

PRD 

100 

X=X+DX 

PRD 

101 

IF  vX.GT.XPCNP;')  GO  TD  65 

PRD 

102 

GD  TD  10 

PRD 

103 

65  CDMTINUE 

PRD 

104 

STOP 

PRD 

105 

PRD 

106 

70  FORMAT  rlOAS) 

PRO 

107 

75  FORMAT  CAb- 13. Fl 0. n> 

PRD 

108 

SO  FORMAT  (5X,F10.0,I5> 

PRD 

109 

85  FORMAT  <3F5= 0> 

PRO 

110 

90  FORMAT  CiOX,  13HCIJRVE  HO.  =  ,  15?  IIH   SIQMA3  =  5F8.£:.  16HHD.  OF 

PDPRD 

111 

IINTS  =  1.15?  18H  FAILURE  RATIO  =  jFS.a-'lOX-  lOHINPUT  DATA^.-EXs 

6HPRD 

11£ 

£SL.HD.  )4K!.   7HSTRAIN1j8Xj   5HS1-S3)6X>   7HSTRAIN3- 3Xj  1 1  H C  PERCEMTPRD 

113 

3  >  > 

PRD 

114 

95  FORMAT  <■  4X .  r£ .  4X  -  3E 1  £ .  3  > 

PRD 

115 

00  FORMAT  (:3X?   7HSTRAIN1 .  SX-   7HSTRAIH3>  5X?   6HSIGMA1  <  7X.   6HTAUDCT.PRD 

116 

16Xj   7HGAMADCT5  6X!   3HMUE?7X.   7HMaDULUS!  4Xi.   6HTANHUEj4Xj   8HS 

IGMPRD 

117 

£ ADCT  J  6X  J   6HT AHMOIi «  6X !   5HRAT I D  s  •-•  > 

PRD 

lis 

05  FORMAT  aiEl£.4> 

PRD 

119 

PRD 

120 

76 


EKSCO..  CI(;SO>! 


END 

SUERQUTINE  SPLINE  CNPN. XPj YP- YDP) 

DIMENSIQN  KP'::SO>!  YP<£0>.  YDPCSO::'!  H<8n::.,  AKi 

iniceu) 

NP1=NPN-1 
DD  5  n=l!NPl 
5  H':M>=HP<M+i>-;«;P<M> 

SLaPl=FD(H(i::'!H';:£::'jYPa>.YP';2>jVP'::3>> 

SLDPN=BD  c:  H  ■;  NPN-S  ;■  s  H (.  NPN- 1  >  .  YP  < NPN-S  >  ?  YP  (:  NPN- 1  >  j  YP < NPH >  > 

CALL  CDFR I T  <  NPN  j  XP  j  YP ?  SL.DP  1 ,  SLDPN «  H I  ?  E I .  C I «  D I  > 

CALL  TR I DGNL  ■-.  NPH .  A I  •  B I  j  C I  •  D I  j  YDP  > 

RETURN 

END 

SUERuUTINE  CGFRIT  <NPN>  XPj  YP?  SLOPl  i.  SLnPN?  AI  -  BI  -  CI  j  DI  > 

DIMENSIDH  yp':.2IJ)>    YP<;£i:i>.  AKeiVJ.-  EI<:eO>5  CICSO?  DI<£0> 

Al'::i>  =  ij.  ij 

BI<i:)  =  <MP'::£>-XPa>>.--3. 

CI '::i>=EI  CI  >.•■■£:. 

Dia::'=(:YP(:s;'-Ypa:)>-'<;xp';£>-xP';i;:o-SLnpi 

AI<NPN>  =  OiPa<PN>-XP'::NPN-l>>.--6. 

bi(:npn::'=hI';npn>*£. 

cI':;npn>=i:i.  u 

n  i  a-ipn  >  =-  <  yp  <  npn  >  -  yp  ■::  npn- 1  ;0  .••'  c  xp  (  npn  >  -kp  ■;  npn- 1  ;0  +sldpn 

Ni=NPN-l 
DD  5  I=£-Hi 

Aia>=';xp<i>-xpa-i>>-'-6, 

EKi>=(:xp<i+i>-xP';i-i>>.''3, 

CI  (I  >  =  <:;«;?<  I  +  i;-XP<  I  >:j.--6. 

Dia:>=(Yp<i+i>-Ypa>::'.'(Xpa+i>-xpa>>-CYP<i>-YP<i-i>v<xpa>-;«;pi 

1   a-l);' 
5  CDNTINUE 
RETURN 


END 

SUBRDUT I NE  TR  I DGNL  (.  NPN  ?  A I  -  E 1 5  CI  s  D I  ?  YDP ) 

DIMENSION    hI<£0>i.    EKeO)?    CKain.    DICED) 

P=EK1> 

q<i:j=-ci<;i>.'P 
u(i>=Dia>.-p 

DD  5  K=£i.NPH 

P=AI'::K>*Q'::K-1>+EI<K> 

q<k>=-ci<;k>/p 

U'::k)=c:di(k>-aick>*u<k-i>)'-'P 
5  cdntinue 

YDP'::NPN)=U(NPN> 
Nl=NPN-i 
DD    iO   L=1jN1 
K=N1+1-L 

ydP'::k>=q<:k;*ydp<k+i>+u(k> 

10  CDNTINUE 
RETURN 

END 

FUNCTION  ET(Rl!R£t.21-Z£?Z::'!)SljS£> 

H.J=Z£-Z1 


YDP ■::  £ n >  ,  Q< 3  0  :j  J  U C 3 0 .j 


I  1=-U.! 


il*( 


T£=0.5*<;S£*';Z3-Z 

T3=CR£-R1>.--HJ 

T4=-';S£-Sl>*HJ.-b 


o!**ej..--Hj 


PRO 

1£1 

SPL 

£ 

SPL 

3 

SPL 

4 

SPL 

5 

SPL 

6 

SPL 

7 

SPL 

8 

SPL 

9 

SPL 

10 

SPL 

11 

SPL 

IE 

SPL 

13 

SPL 

14 

CDF 

e 

CDF 

3 

CDF 

4 

CDF 

5 

CDF 

6 

CDF 

7 

CDF 

s 

CDF 

9 

CDF 

10 

CDF 

11 

CDF 

IE 

CDF 

13 

CDF 

14 

CDF 

15 

CDF 

16 

-■CDF 

1? 

CDF 

18 

CDF 

19 

CDF 

£0 

CDF 

£1 

CDF 

££ 

TRI 

£ 

TRI 

3 

TRI 

4 

TRI 

5 

TRI 

6 

TRI 

7 

TRI 

a 

TRI 

9 

TRI 

10 

TRI 

11 

TRI 

1£ 

TRI 

13 

TRI 

14 

TRI 

15 

TRI 

16 

TRI 

1? 

TRI 

18 

TRI 

19 

TRI 

£0 

ET 

£ 

ET 

3 

ET 

4 

ET 

5 

ET 

& 

ET 

7 

77 


ET=T1+T£+T3+T4 
RETURN 

END 

FUNCT I DH  FD  ( S  i !  SS »  R 1 J  R£ )  F'3  > 
IF  <S1-SS>  5? 10.5 
5  FD=<Re-Rl>.-Sl 
RETURN 
iO  Fri=<-3.*Ri+4.*RS-R3>.-'':;2.*Sl> 
RETURri 

END 

FUNCTION    EriCSl.Se>RljR2)R3;' 
IF    <:S1-S8>    5.10.5 
5   BD=';R3-R£>'-'S£ 
RETURN 
10   ED=<3,*R3-4.>i*R£+Rl>.-CS.*Sn 
RETURN 

END 

FUNCT  I  ON   DRIi  I  NET  <  A ;  E «  C .  D  <  E .  F ..  P  > 

HJ=D-C 

El=ri~P 

B£=P-C 

A1=B1*W3 

A£=B£**3 

Tl=Hl*H.'(:6,*HJ> 

T£=H£*B--<6.*HJ> 

T3= <E-A*H J*:-*£''b .  > * < D-P > •■■'H J 

T4=<F-E*HJ**£-'6.  >*<P-C>''HJ 

DRD I NET=T 1 +T£+T3+T4 

RETURN 

END 


ET 

8 

ET 

9 

ET 

10 

ET 

11 

FD 

£ 

FD 

3 

FD 

4 

FS 

5 

FD 

6 

FD 

7 

FD 

e 

FD 

9 

ED 

£ 

ED 

3 

ED 

4 

ED 

5 

ED 

6 

ED 

7 

BD 

,q 

BD 

9 

DRD 

E 

DRD 

o 

DRD 

4 

DRD 

5 

DRD 

6 

DRD 

7 

DRD 

8 

DRD 

9 

DRD 

10 

DRD 

11 

DRD 

1£ 

DRD 

13 

DRD 

14 

DRD 

15 

